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ABSTRACT 


This  thesis  proposes  a  framework  by  which  optimal  system  renewal  decisions 
may  be  made  in  a  consistent  and  timely  manner  within  the  military,  under  realistic 
conditions  of  changing  environments.  A  unique  network  representation  of  the  system 
renewal  process  was  used  to  develop  a  prototype  version  of  the  analytical  core  model. 
Its  plausibility  and  usefulness  was  demonstrated  by  a  series  of  case  studies.  The  case 
studies  also  show  how  pertinent  staff  forecasts  can  be  organised  and  integrated  to 
provide  decision  makers  with  a  broad,  consistent,  long-term  perspective  of  the  issues 
relevant  to  system  renewal  planning.  They  are  presented  with  a  graphic  picture  of 
the  entire  solution  space  as  structured  by  the  scenarios  considered.  Various  solutions 
are  suggested  for  each  scenario  and  their  robustness  may  be  tested  by  thoroughly 
exercising  the  model  for  a  wide  range  of  scenarios.  Prediction  of  what  and  when 
renewals  are  likely  to  be  and  estimation  of  the  associated  cost  and  effectiveness 
allow  anticipatory  long-term  plans  to  be  formulated.  The  thesis  also  suggests  how 
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I.  INTRODUCTION 


A.  INTRODUCTION 

Military  forces  continually  undergo  change  as  obsolete  or  uneconomical  assets 
are  retired,  new  weapons  are  introduced  and  existing  ones  modified  or  upgraded. 
This  process  of  system  renewal  is  affected  by  a  myriad  of  economic,  technological 
and  politico-military  considerations.  Changes  in  enemy  capabilities  or  doctrine,  for 
instance,  may  render  certain  systems  less  effective  or  even  obsolete.  New  technology 
may  offer  opportunities  to  improve  present  capabilities.  Most  of  the  available  open 
literature,  analyses  and  models  for  system  renewal  are  either  one-dimensional,  too 
general,  or  too  commercially  oriented  to  be  useful  to  military  planners.  There  is  a 
need  to  synthesize  a  coherent,  balanced  view  relevant  to  military  planners  to  provide 
comprehensive  policy  guidance  on  upgrading,  replacement  and  related  matters. 

B.  PROBLEM 

The  system  renewal  planning  problem  is  encountered  prior  to  the  last  milestone 
of  the  defence  acquisition  cycle.  It  is  basically  one  of  being  able  to  evaluate  and 
decide  in  a  timely  and  rational  manner,  under  conditions  of  changing  technological, 
economic  and  politico-military  environments,  both  the  following: 

1.  When  w'ill  a  defence  asset,  in  its  present  form,  become  obsolete  or  no  longer 
cost-efTective  to  retain  (relative  to  existing  or  emerging  alternatives)? 

2.  What  should  then  be  done7  What  are  the  alternatives  and  which  one  should 
be  chosen? 


1 


C.  DEFINITION 


System  renewal  encompasses  all  changes  made  to  maintain  or  improve  the 
cost-effectiveness  of  the  present  defence  system.  It  includes  the  acquisition  of  new 
capabilities,  retirement  from  service  of  obsolete  ones,  relegation  to  reserves  of  sur¬ 
plus  capabilities  and  conversion  of  systems  to  new  roles.  This  view  streamlines  the 
treatment  of  changes  to  the  status  quo  and  involves  a  higher  level  of  system  repre¬ 
sentation.  It  allows  us  to  consider  complex  changes,  such  as  when  a  new  support 
ship  displaces  a  couple  of  smaller  ones,  one  of  which  is  preserved  for  the  boneyard 
and  the  other  converted  to  a  hospital  ship. 

D.  PURPOSE 

The  purpose  of  this  thesis  is  to  develop  a  rational  and  consistent  decision 
model  for  system  renewal  planning  that  is  relevant  and  useful  to  military  planners. 
Various  issues  pertinent  to  system  renewal  are  integrated  into  a  common  analytical 
framework  to  make  for  a  more  comprehensive  and  better  formalised  concept  of 
obsolescence  in  military  systems. 

This  work  is  intended  as  a  background  study  for  later  development  of  a  com¬ 
prehensive  policy  to  complement  life-cycle  management  initiatives  in  the  Singapore 
Ministry  of  Defence.  It  serves  as  a  rudimentary  framework  that  can  be  expanded 
later  for  actual  use. 

E.  IMPETUS 

The  choice  of  this  research  topic  was  motivated  by  the  following: 

•  Singapore's  defence  buildup  has  been  underway  for  two  decades  and  already 
system  renewal  decisions  are  due  for  some  of  the  aging  systems  that  were  ac¬ 
quired  initially.  A  framework  that  facilitates  taking  the  long-term  perspective 
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in  system  renewal  planning  will  help  ensure  that  such  optimal  decisions  are 
made  in  a  systematic  and  timely  manner.  As  the  assets  involved  are  consid¬ 
erable,  improvements  in  system  renewal  planning  is  of  significant  importance. 

•  Milestone  5  (upgrading/replacement)  was  introduced  in  U.S.  DoD  Instruc¬ 
tion  5000.2  in  September  1987.  No  clear  guidance  has  been  provided  on  the 
approaches  to  be  adopted  for  this  new  milestone. 

•  Available  literature  generally  either  deals  with  replacements  aimed  at  pre¬ 
empting  failure  (preventive  maintenance)  or  reduces  the  problem  to  a  purely 
economic  one  (commercial  approach).  There  is  a  dearth  of  literature  of  rele¬ 
vance  to  military  system  renewal  planning. 

F.  SCOPE 

Replacement  policies  concerning  routine  preventive  maintenance  are  not  ex¬ 
amined  in  this  thesis.  The  scope  is  limited  to  major  systems  or  forces  for  which 
the  considerable  effort  involved  in  identifying,  tracking  and  analysing  capability  re¬ 
quirements,  costs  and  measures  of  effectiveness  (MOEs)  are  justifiable.  At  the  heart 
of  the  framework  to  be  presented  is  an  analytical  model,  the  workings  of  which 
will  be  illustrated  under  various  assumptions  based  on  a  hypothetical  weapon  sys¬ 
tem.  Attention  is  focused  on  single  mission  systems  in  changing  environments.  The 
treatment  of  threat  evolution,  force  level  planning  and  higher  level  optimisation 
are  briefly  discussed.  Simple  sensitivity  analysis  is  illustrated  but  full  and  proper 
stochastic  modelling  and  risk  analysis  lie  outside  the  scope  of  this  thesis. 
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II.  ISSUES  AND  LITERATURE  SURVEY 

Systems  are  renewed  as  they  become  obsolete.  Why?  When  is  a  system 
considered  obsolete?  What  are  the  issues  involved  in  planning  for  system  renewal? 
These  are  the  larger  questions  addressed  in  this  chapter.  We  begin  by  explaining 
obsolescence  and  its  consequences.  Next,  for  a  broad  overview  of  the  system  renewal 
problem,  various  issues  pertaining  to  the  selection  and  scheduling  of  alternatives  are 
discussed.  The  chapter  ends  with  a  survey  of  some  of  the  research  that  have  been 
carried  out  on  system  renewal. 

A.  OBSOLESCENCE 

Obsolescence  may  be  categorised  as  follows: 

•  Economic:  A  system  is  considered  economically  obsolete  when  it  becomes  less 
cost-efTective  than  available  alternatives  because  of  cost  escalation  attributed 
to  physical  deterioration,  scarcity  of  spares,  inflation,  etc.  An  old  truck  that 
breaks  down  often,  guzzles  fuel  that  has  become  expensive,  or  has  little  spares 
support  is  economically  obsolete. 

•  Technological:  This  form  of  obsolescence  occurs  when  innovation  makes  avail¬ 
able  a  new,  more  cost-effective  alternative  to  an  incumbent  system.  Fighter 
planes  with  propellers,  for  example,  became  technologically  obsolete  when  jet- 
propelled  ones  became  available. 

•  Operational:  Operational  obsolescence  arises  when  a  system  fails  to  meet  its  re¬ 
quirements  because  of  deterioration  in  performance,  enemy  counter-measures 
or  shifts  in  own  strategic/tactical  doctrine.  Immobile  defence  systems,  for 
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instance,  may  become  obsolete  if  it  is  known  that  they  can  be  effectively  en¬ 
gaged  or  bypassed  by  the  enemy,  or  if  strategic  emphasis  shifts  from  attrition 
to  maneuver  warfare. 

Obsolescence  may  be  a  compound  effect  of  the  above.  A  gun  that  is  no  longer 
accurate  may  also  be  uneconomical  to  maintain. 

The  continued  operation  of  an  obsolete  system  is  always  a  waste  of  money, 
but  technological  and  operational  obsolescence  may  have  more  serious  consequences. 
Failure  to  exploit  technological  opportunities  or  remedy  operational  deficiencies  may 
lead  to  military  disasters.  Technology  is  not  exploited  merely  by  keeping  pace  with 
research  as  was  shown  by  the  French  experience  in  1939.  The  French  Army  then 
had  for  some  time  developed  prototypes  of  aircraft,  armor  and  antitank  weapons 
superior  to  those  of  the  German  Army.  Their  wait  for  still  better  ones  to  come  along 
before  introducing  them  into  the  inventory  was  a  reason  for  their  unpreparedness 
[ Possony  and  Pournelle,  1970],  The  Polish  Army’s  fatal  cling  to  cavalry  is  a  clear 
lesson  in  the  consequence  of  operational  obsolescence. 

Not  all  obsolescence  is  bad.  The  term  “built-in  obsolescence”  is  sometimes 
used  for  the  practice  of  deliberately  designing  machines  that  deteriorate  rapidly 
after  a  given  life;  the  user  being  compensated  for  this  by  a  lower  purchase  price. 
This  can  be  seen,  for  example,  in  commercial  batteries  with  different  lives. 

B.  STRATEGIC  ISSUES 

Defence  forces  exist  for  a  purpose,  namely  to  deter  aggression  and,  if  deterence 
fails,  to  successfully  wage  war.  The  strategy  to  achieve  this  purpose  should  guide 
the  formulation  of  system  renewal  plans.  Potential  adversaries  should  be  identified 
and  their  intentions,  war-proneness,  capabilities,  weaknesses  and  strategies  assessed 


to  arrive  at  a  strategy  that  is  congruent  with  foreign  policy.  Competitive  strat¬ 
egy  [Englund,  1987],  for  example,  seeks  a  favorable  cost  exchange  ratio  by  making 
obsolete  major  components  of  the  enemy's  defence  system.  Stealth  bombers  and 
SDI  are  manifestations  of  this  potentially  destabilising  strategy.  Some  references  on 
strategic  analysis  and  forecasting  are  provided  in  the  bibliography. 

There  are  many  facets  of  strategy  that  affects  system  renewal  plans.  Resources 
for  modernisation  may  be  channelled  to  boost  readiness  if  conflict  seems  imminent. 
The  danger  of  provoking  or  worsening  an  arms  race  must  be  weighed  carefully  if  plans 
call  for  a  significantly  increased  capability.  The  impact  of  a  plan  on  the  economy  may 
also  be  significant.  The  B-2  program  recommended  by  the  USAF,  for  example,  asks 
for  a  whopping  seventy-five  billion  dollars.  The  impact  of  modest  levels  of  defence 
spending  on  the  economy  is  usually  seen  to  be  beneficial.  Fredericksen’s  [1989] 
analysis,  for  example,  showed  a  favorable  feedback  relationship  between  defence 
spending  and  economic  growth  in  Indonesia  for  the  period  1964  to  1985.  Economic 
issues  and  methodology  in  arms  race  analysis  are  surveyed  in  Leidy  and  Staiger 
[1985]. 

The  link  between  strategy  and  procurement  is  historically  lacking  in  the  US 
(and  Singapore).  The  chairman  of  a  Defense  Resource  Management  Study  in  1979, 
for  example,  declared  that  there  is  broad  agreement  that  the  first  P  in  PPBS  is 
silent  [Rice,  1979].  The  Packard  Commission  in  1986  stated: 

. . .  there  is  a  great  need  for  improvement  in  the  way  we  think  through  and  tie 
together  our  security  objectives,  what  we  spend  to  achieve  them,  and  what  we 
buy.  The  entire  undertaking  for  our  nation’s  defence  requires  more  and  better 
long-range  planning.  This  will  involve  concerted  action  by  our  professional 
military,  the  civilian  leadership  of  the  Department  of  Defense,  the  President, 
and  the  Congress. 
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Today,  there  is  no  rational  system  whereby  the  executive  branch  and  the 
Congress  reach  coherent  and  enduring  agreement  on  national  military  strategy, 
the  forces  to  carry  it  out,  and  the  funding  that  should  be  provided  -  in  light  of 
the  overall  economy  and  competing  claims  on  national  resources.  The  absence 
of  such  a  system  contributes  substantially  to  the  instability  and  uncertainty 
that  plagues  our  defence  program.  These  cause  imbalances  in  our  military 
forces  and  capabilities  and  increases  the  cost  of  procuring  military  equipment. 

This  thesis  seeks  to  improvise  a  link  between  strategy  and  force/resource  plan¬ 
ning  by  incorporating  strategic  outlook  into  system  renewal  plans.  The  cost  and 
effectiveness  of  alternatives  of  interest  are  projected  into  the  future  to  make  for 
more  forward-looking  plans.  Figure  2.1  illustrates  the  links  between  strategic  issues 
and  system  renewal  planning.  It  is  hoped  that  at  some  future  time  system  renewal 
plans  can  be  aggregated  and  linked  to  arms  race  models,  strategic  analysis  models, 
econometric  models,  etc.  to  provide  a  higher  level  of  optimisation. 

Technology  pushes  system  renewal  because  it  holds  the  key  to  one’s  edge  over 
the  enemy.  However,  excessive  emphasis  on  performance  exacts  an  inordinate  cost. 
Figures  2.2  and  2.3  [Gansler,  1989]  show  the  cost  histories  and  performance  growth 
respectively  of  a  series  of  U.S.  fighter  and  attack  aircraft.  The  performance  assess¬ 
ment  methodology  used  in  Figure  2.3  is  explained  in  Chapter  4.  As  can  be  seen, 
the  difference  in  unit  cost  between  high-end  aircraft  (e.g.  F-15A,  F-18)  and  low-end 
ones  (e.g.,  F-16A,  A-7A)  had  grown  significantly  through  the  years.  Performance 
growth  was  however  much  more  modest,  indicating  that  cost-effectiveness  of  high- 
end  aircraft  had  deteriorated  relative  to  that  of  the  low-end  ones.  Figure  2.3  suggests 
that  the  performance  of  the  first  operational  units  of  a  low-end  aircraft  lags  that  of 
high-end  ones  bv  only  five  years.  It  also  implies  that  technological  superiority  over 
Soviet  counterparts  could  have  been  maintained  using  lower-risk,  lower-cost  aircraft. 
The  lower  cost  translates  into  more  aircraft,  which  may  matter  more  in  terms  of 
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military  results.  Expensive  high-performance  or  multi-role  assets  also  require  in¬ 
tricate  defences  as  they  are  too  valuable  to  lose.  Quality  must  also  be  traded  off 
against  quantity  when  selecting  alternatives  in  a  system  renewal  plan.  This  tradeoff 
can  be  accomplished  within  the  model  to  be  introduced  in  the  following  chapters  by 
denoting  quantity  as  a  performance  parameter.  The  impact  of  future  innovation  on 
the  cost-effectiveness  of  alternatives  of  interest  is  projected  so  that  optimal  numbers 
of  the  right  alternatives  will  be  procured  at  the  right  time. 

The  question  of  how  far  to  push  technology  in  a  procurement  program  has 
major  cost  and  effectiveness  implications.  Figure  2.4  illustrates  how  the  technolog¬ 
ical  frontier  for  a  certain  type  of  system  may  advance  with  R&D  effort.  Point  A  is 
typical  of  a  target  set  to  ensure  the  ability  to  counter  probable  threats  over  a  long 
horizon.  Due  to  the  risks  involved,  it  is  likely  in  this  case  to  end  up  somewhere  inside 
the  area  bounded  by  the  points  A,  B  and  C.  The  less  ambitious  and  better  plan 
would  be  to  seek  to  quickly  and  cheaply  arrive  at  say  point  D  first,  then  improve 
incrementally  to  point  A.  Having  such  pre-planned  product  improvements  ( P3I ) 
requires  good  forecasts  of  likely  future  upgrades.  The  model  to  be  introduced  facil¬ 
itates  the  selection  and  scheduling  of  such  upgrades.  The  reader  may  be  interested 
in  the  1984  report  by  the  U.S.  Defense  Science  Board  entitled  “Improved  Defense 
through  Equipment  Upgrades:  The  U.S.  and  Its  Security  Partners”. 

Requirement  specifications  largely  determine  the  technological  content  of  a 
program.  These  requirements  may  not  be  the  result  of  careful  study.  In  fact,  they 
understandably  tend  to  reflect  the  psychological  comfort  levels  of  those  who  have  to 
deal  with  the  exigencies  of  war.  After  all,  it  never  hurts  to  be  able  to  go  the  extra 
mile.  Accordingly,  the  best  of  various  candidate  systems  often  becomes  the  “re¬ 
quirement”.  The  importance  of  specifying  objective,  mission-oriented  requirements 
cannot  be  overs  tressed. 
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Figure  2.1  :  System  Renewal  and  Related  Issues 
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100th  UNIT  COST  (MILLIONS  87$) 


YEAR  OF  INITIAL  OPERATIONAL  CAPABILITY 


Figure  2.2  :  Cost  Growth  in  U.S.  Fighter/Attack  Aircraft 
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PERFORMANCE  POTENTIAL 


YEAR  OF  INITIAL  OPERATIONAL  CAPABILITY 


Figure  2.3  :  Performance  Comparison  of  U.S. 

&  Soviet  Fighter/Attack  Aircraft 
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C.  ECONOMIC  ISSUES 

Cost-effectiveness  of  alternatives  are  dependent  on  future  interest  rate,  foreign 
exchange  rates,  inflation  rate,  resale/salvage/scrap  market  condition  and  budget 
allocations.  The  U.S.  Department  of  Defense  publishes  forecasts  of  inflation  indices 
for  various  categories  of  military  hardware  which  can  be  used  in  system  renewal 
planning.  DoD  would  of  course  like  to  estimate  conservatively  (high  inflation)  so 
as  to  get  more  money  to  play  safe  in  case  it  actually  is  high.  However,  in  times  of 
high  inflation,  there  is  in  fact  considerable  political  pressure  for  lower  estimates  in 
order  to  suppress  defence  budget  growth  projections.  The  cumulative  inflation  from 
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1978  to  1981  was  45  percent  ,  while  the  inflation  estimated  in  the  DoD  budget  was 
only  28  percent  [Gansler,  1989J.  Therefore  one  must  exercise  discretion  in  using  the 
inflation  indices. 

Assuming  a  common  interest  rate  for  all  programs  does  not  make  for  equitable 
comparison  as  alternatives  with  different  life-spans  and  cost  profiles  are  affected 
differently.  One  should  try  to  forecast  the  rate  as  accurately  as  possible  to  avoid 
biased  results.  As  economic  factors  such  as  interest  rates,  foreign  exchange  rates 
and  market  conditions  are  very  difficult  to  estimate,  sensitivity  of  solutions  to  their 
likely  fluctuations  should  be  tested  in  actual  system  renewal  planning. 

Adoption  of  a  longer  term  perspective  in  system  renewal  planning  should  be 
matched  by  the  adoption  of  a  multi-year  budget.  Procurement  of  needed  quantities 
of  a  system  is  usually  most  economical  at  a  certain  rate.  An  optimal  system  renewal 
plan  may  involve  humps  in  budgetary  outlays.  The  minimum  efficient  rate  of  pro¬ 
duction  for  F-15  aircraft,  for  example,  is  120  per  year.  However,  DoD  was  allowed 
only  41  per  year  [Gansler,  1989].  This  U.S.  practice  of  stretching  out  programs  to 
smooth  out  the  humps  tends  to  lead  to  gross  inefficiencies  in  procurement.  The 
floating  of  treasury  bonds,  for  example,  may  be  a  better  way  to  achieve  a  more 
efficient  procurement  program.  For  multi-year  budgets  to  work,  accurate  cost  esti¬ 
mates  are  necessary.  Buying  off-the-shelf  items,  adopting  low-risk  approaches  (e.g., 
by  emphasizing  P3I)  and  designing  to  cost  are  collaborative  means  towards  this 
end. 

D.  MODELLING  ISSUES 

Military  system  renewal  planners  have  to  grapple  with  problems  of  multi¬ 
facetedness,  vagueness,  uncertainty  and  subjectivity.  Even  if  the  effectiveness  of  a 
given  alternative  can  be  summarised  by  a  single  MOE,  one  has  still  to  consider  its 
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cost  and  probable  service  life.  The  system  renewal  problem  has  therefore  at  least 
three  dimensions  -  effectiveness,  cost  and  useful  life.  Often,  service  life  is  simply 
assumed  and  net  cost  is  amortised  into  a  rate  over  this  period  for  evaluations.  How¬ 
ever,  advocates  of  an  alternative  will  tend  to  assume  a  longer  life  than  opponents. 
Instead  of  assuming  when  an  alternative  will  become  obsolete,  the  model  to  be  pre¬ 
sented  derives  the  service  life  as  an  output  based  on  economic,  technological  and 
military  (requirement)  forecasts.  This  should  make  for  more  realistic  and  equitable 
evaluations. 

Risks  due  to  uncertainty  can  be  assessed  by  testing  the  sensitivity  of  solutions 
to  possibilities  learned  from  experience.  Dominance  or  robustness  of  a  solution 
would  provide  some  assurance  of  its  validity  under  a  range  of  expected  scenarios. 
The  quality  of  experience  is  of  critical  importance  for  such  risk  assessment.  As  such, 
qualified  experts  should  be  consulted  in  the  formulation  of  scenarios. 

Besides  multi-facetedness  and  uncertainty,  a  planner  has  to  contend  with  prob¬ 
lems  of  vagueness.  Strategic  analysis  may  highlight  the  need  to  heighten  readiness  in 
the  coming  years  but  such  qualitative  goals  do  not  readily  lend  themselves  to  math¬ 
ematical  modelling.  Even  if  readiness  can  be  precisely  quantified  by  a  parameter, 
♦here  may  not  be  consensus  concerning  what  constitutes  the  proper  level  of  readi¬ 
ness.  Fuzzy  set  and  logic  theory  offer  means  to  bridge  the  gap  between  quantitative 
and  qualitative  reasoning,  thus  linking  strategy  formally  to  procurement,  planning. 
The  efforts  of  J.  Dockery  [1985]  is  an  example  of  work  done  in  this  area.  Karprzyk 
and  Oriovski  [1987]  and  Zimmerman,  ct  ah,  [1984]  also  discuss  new  paradigms  more 
suited  to  situations  where  vagueness  and  uncertainty  are  clearly  present.  The  prob¬ 
lem  of  vagueness  is  not  considered  in  this  thesis  but  may  be  addressed  in  future 
embellishments. 
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Although  the  exercise  of  human  judgment  is  essential  in  system  renewal  plan¬ 
ning,  the  subjectivity  involved  must  be  managed  well  to  reduce  the  possibility  of  bias 
(conscious  or  otherwise).  The  problems  encountered  in  assessing  utility  functions, 
for  example,  were  surveyed  by  Hershey  et  al.,  (1982)  and  Farquhar  [  198-4] .  It  should, 
however,  be  noted  that  intuitive  (subjective)  judgment  of  experts  may  in  certain 
cases  prove  superior  to  analytical  studies  (Hammond,  1987],  Human  behavioral  as¬ 
pects  are  important  in  a  real  application,  but  they  are  not  considered  in  this  thesis 
because  the  cases  studied  are  hypothetical  ones  meant  for  illustrating  the  model 
that  was  developed.  Utility  functions,  for  example,  were  assumed  to  be  explicitly 
defined  even  though  such  functions  are  difficult  to  assess  accurately  in  practice. 

The  horizon  of  optimisation  is  an  important  issue  in  modelling  system  renewal. 
Very  often,  this  horizon  is  either  fixed  at  some  arbitrary  value  or  taken  to  be  in¬ 
finity.  Both  practices  suffer  from  serious  flaws.  Choosing  horizons  that  are  a  small 
non-integral  number  of  life-spans  would  produce  incorrect  results.  This  occurs,  for 
instance,  if  the  horizon  for  a  system  with  a  typical  lifespan  of  ten  years  is  set  to, 
say,  twenty-five  years.  Assuming  an  infinite  horizon  is  clearly  unrealistic.  There  are 
procedures  that  allow  one  to  find  the  horizon  beyond  which  the  decision  for  the  first 
change  remains  optimal.  However,  these  so-called  “planning  horizon  procedures’’ 
are  essentially  also  based  on  the  infinite  horizon  assumption.  This  thesis  takes  a 
different  approach,  namely  the  finite-stage  approach  (from  dynamic  programming), 
in  which  the  optimisation  horizon  is  a  variable  output  given  by  the  years  spanned 
by  a  sequence  of  systems. 

Modelling  the  effectiveness  of  a  military  system  is  sometimes  confused  with 
predicting  the  outcome  of  its  use  in  battle.  The  two  are  related  but  not  synonmous. 
As  an  example,  let  us  suppose  that  a  certain  armour  unit  has  the  MOE  structure 
shown  in  Figure  2.5.  The  outcome  in  a  given  battle  scenario  may  be  described 
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by  the  number  of  own  and  enemy  losses.  Neither  of  these  numbers  expresses  fully 
the  effectiveness  of  the  force.  An  abstract  figure  of  merit  would  have  to  be  com¬ 
puted  from  these  numbers  based  on  judgment  concerning  what  constitutes  success. 
Therefore,  one  might  say  that  battle  outcomes  are  predicted  by  combat  modelling 
whereas  effectiveness  is  evaluated  by  “combat  value”  modelling.  The  abstract  figure 
of  merit  may  be  said  to  be  the  MOE,  and  performance  parameters  that  determine 
this  value  may  be  said  to  be  the  measures  of  performance  (MOPs).  This  distinction 
is,  however,  an  artificial  one  as  the  MOPs  may  collectively  be  considered  to  be  a 
vector-valued  MOE.  Accordingly,  MOE  and  MOP  in  this  thesis  are  both  taken  to 
mean  the  same  thing  -  key  performance  parameters  that  determine  the  effective¬ 
ness  of  a  system.  In  the  example  shown  in  Figure  2.5,  the  MOEs  would  be  range, 
accuracy  and  lethality  if  these  parameters  are  the  ones  that  have  major  effects  on 
effectiveness. 

E.  ORGANISATIONAL  ISSUES 

System  renewal  plans  may  be  developed  for  many  different  levels,  ranging  from 
a  single  equipment  to  an  entire  task  force.  Unless  they  are  coordinated  such  that  the 
objectives  of  upper  echelons  are  considered  by  lower  echelon  planners,  sub-optimal 
or  even  counterproductive  plans  are  likely  to  be  made.  Most  resource  allocation 
models  do  not  consider  the  organisation  as  a  factor  in  the  decision  process;  that 
is,  they  assume  a  monolithic  decision-making  entity.  The  same  decision-makers  are 
assumed  whether  the  subject  of  planning  is  a  single  artillery  gun  or  an  infantry  divi¬ 
sion.  Ruefli  [1974]  surveyed  several  economic  and  management  science  models  that 
provide  remedy  for  this  unrealistic  assumption  by  explicitly  reflecting  organisational 
structure  in  the  model.  The  Dantzig- Wolfe  model  (Dantzig  and  Wolfe,  19G0]  is  the 
best  known  of  the  classical  models.  In  more  recent  years,  models  such  as  the  one  by 
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Figure  2.5  :  Hierarchy  of  MOEs 

Source  :  Modified  version  of  the  diagram  from 
"Notes  on  MOEs”  by  E.B.Rockower,  Aug.  1985. 


Cassidy  et  al.,  [1971]  that  are  based  on  behavioral  theories  have  emerged.  Organ¬ 
isational  issues  are  not  considered  in  this  thesis  as  the  model  is  still  a  conceptual 
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one  that  is  not  ready  for  organisation-wide  implementation. 

F.  LITERATURE  SURVEY 

No  literature  suitable  for  military  system  renewal  planners  could  be  found  as 
all  those  surveyed  were  commercially  oriented  ones  that  subjugated  performance 
considerations  under  cost.  Most  of  them  consider  only  replacement  alternatives. 
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The  earliest  of  these  were  based  on  simple  calculations  and  differential  and  integral 
calculus.  Examples  are  Preinrich’s  Constant  Chain  Model  [1940]  and  Terborgh’s 
MAPI  system  [1949].  In  more  recent  years,  methods  such  as  dynamic  program¬ 
ming,  Markov  chains  and  Pontryagin’s  maximum  principle  have  become  more  pop¬ 
ular.  The  first  such  model  was  introduced  by  R.  Bellman  [1955]  using  dynamic 
programming.  Howard  [1960]  used  a  policy-iteration  method  to  solve  multi-stage 
decision  problems  in  connection  with  Markov  processes.  Numerous  articles  conern- 
ing  stochastic  models  were  published  in  the  1960’s.  A  survey  of  the  treatment  of 
stochastically  failing  equipment  was  compiled  by  Jorgenson  et  al.,  [1967].  Naslund 
[1966]  was  the  first  to  apply  the  maximum  principle  to  the  replacement  problem. 
Thomson  [1968],  Rapp  [1969],  Sethi  [1970]  and  many  others  used  these  so-called 
“optimal  control”  theory  models  that  uses  the  maximum  principle.  There  are  nu¬ 
merous  variants  in  replacement  models,  the  assumptions  of  which  were  grouped  by 
Luxhoj  and  Jones  [1986]  in  a  general  framework.  Some  of  the  works  not  referenced 
in  this  survey  may  be  found  in  the  bibliography. 

Some  of  the  more  noteworthy  recent  publications  that  explicitly  considered 
technological  advances  are  as  follows: 

1.  R.C.  Stapleton  et  ah,  [1972]  examined  the  effects  of  different  forms  of  tech¬ 
nological  change  (change  in  purchase  price,  resale  value,  operating  costs  and 
rate  of  deterioration)  on  costs  using  dynamic  programming.  This  work  pro- 
vided  the  author  of  this  thesis  with  .basic  ideas  on  how  to  model  the  impact 
of  technological  progress  on  cost. 

2.  S.  Chand  and  S.  Sethi  [1979  and  1982]  applied  the  planning  horizon  procedure 
developed  by  H.M.  Wagner  and  T.M.  Whitis  [1958]  for  the  dynamic  lot  size 
model  to  the  replacement  problem.  It  allows  one  to  obtain  the  first  replacement 
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time  for  an  infinite  planning  horizon  based  on  the  optimal  policy  for  a  finite 
horizon.  Technological  advances  is  forecasted  for  this  equivalent  finite  horizon. 
The  1982  paper,  which  extends  the  first  to  the  case  of  multiple  alternatives, 
treated  upgrades  within  the  same  framework  as  replacement  alternatives.  This 
idea  led  to  the  adoption  of  a  common  framework  for  all  alternatives  in  this 
thesis. 

3.  T.  Goldstein  et  al.,  [1988]  examined  the  influence  of  a  forthcoming  new  and 
better  machine  (based  on  an  expected  technological  breakthrough)  on  the 
system  renewal  decision.  Each  year,  if  the  new  machine  does  not  appear,  a 
decision  is  made  to  either  continue  with  the  incumbent  and  wait  for  the  bet¬ 
ter  machine,  or  buy  a  new  version  of  the  incumbent.  The  model  determines 
the  optimal  age  for  replacing  the  existing  machine.  This  model  ignores  evolu¬ 
tionary  changes  that  characterises  technological  innovation  and  thus  was  not 
adopted  in  this  thesis. 

4.  Y.  Kusaka  [1985]  showed  a  method  for  deciding  whether  to  keep  or  replace 
the  existing  machine  at  the  present  time  without  determining  the  sequence  of 
subsequent  replacement  times.  A  joint  paper  with  II.  Suzuki  [1988]  expanded 
on  the  earlier  evaluation  system  by  introducing  the  concept  of  “control  limit” 
policy.  Kusaka’s  method  was  not  adopted  in  this  thesis  as  the  provision  of 
estimates  concerning  future  renewals  is  deemed  to  be  a  useful  feature. 
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III.  THE  MODEL 


A.  INTRODUCTION 

As  evident  in  the  previous  chapter,  many  complex  issues  impinge  upon  de¬ 
cisions  concerning  systems  renewal  planning.  Literature  searches  have  failed  to 
uncover  any  model  suitable  for  system  renewal  planning  in  the  military.  The  dy¬ 
namic  system  renewal  planning  (DSRPM)  model  to  be  presented  in  this  chapter 
is  a  first-cut  attempt  to  weave  a  unifying  framework  spanning  some  of  the  more 
quantitative  issues  involved  in  systems  renewal  planning.  The  overall  approach  is 
outlined  first,  followed  by  theoretical  underpinnings  of  the  model. 

B.  APPROACH 

The  adopted  approach  is  to  use  a  versatile  methodology,  namely  dynamic  pro¬ 
gramming  augmented  by  the  Lagrangian  relaxation  technique,  to  explicitly  model 
the  effects  of  aging,  technological  advancement,  and  changes  in  capability  require¬ 
ments.  The  methodology  operates  on  an  acyclic  directed  network  of  the  kind  shown 
in  Figure  3.1.  This  network  representation  of  the  system  renewal  process  is  unique 
to  this  thesis.  As  this  network  lies  at  the  heart  of  the  model,  a  brief  description  is 
in  order  before  we  move  on  to  the  other  parts  of  the  model. 


20 


1996 


Note  1  :  Nodes  denote  possible  transitions  while  arcs 
represent  possible  alternatives  and  system  life  spans 
The  dotted  arc  from  node  3  to  7,  for  instance,  stands 
forthe  possibility  of  acquiring  alternative  2  in  1990  and 
using  it  as  System  8  until  1992 
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Note  2  .  This  diagram  is  just  an  illustration  In  real 
applications,  the  number  of  arcs,  nodes  and  years 
spanned  would  be  higher.  The  years  spanned  by 
the  network  also  need  not  be  contiguous. 
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Note  3  :  System  A  and  Alternative  0  are 
synonymous  with  the  present  system. 
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Figure  3.1  :  Network  Representation  of  System  Renewal  Process 
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The  network  represents  the  system  renewal  process  involving  the  selection  of  a 
succession  of  systems  to  fulfill  some  capability  requirements.  The  root  node  1  repre¬ 
sents  the  present  system  (System  A)  in  the  current  year.  Each  of  the  arcs  emanating 
from  it  represent  a  possible  residual  life-span  of  this  system.  These  arcs  terminate 
on  a  column  of  nodes  representing  possible  points  in  time  for  transition  to  the  next 
system  (System  B).  Transitions  may  be  replacements  or  upgrades.  Similarly,  arcs 
extending  from  these  nodes  to  an  adjacent  column  of  nodes  represent  the  possible 
years  for  transition  from  System  B  to  System  C.  The  number  of  arcs  connecting  each 
pair  of  nodes  corresponds  to  the  number  of  alternatives  that  may  be  selected  for 
System  B;  likewise  for  System  C.  Thus  if  cumulative  cost  and  effectiveness  of  using 
a  given  alternative  over  a  given  life-span  are  assigned  as  attributes  to  the  associated 
arc,  optimisation  for  the  three  systems  reduces  to  one  of  finding  the  “best”  path 
through  the  network. 

Finding  the  “best”  path  constitutes  solving  a  three  stage  dynamic  program¬ 
ming  problem.  If  cost  is  set  to  a  very  large  number  whenever  an  alternative  ages 
such  that  one  or  more  performance  parameters  at  the  end  of  a  given  life  span  fails 
to  meet  stipulated  requirements,  then  dynamic  programming  yields  the  solution 
that  meets  the  requirements  at  minimum  cost.  However,  if  one  seeks  the  optimal 
trade-off  between  cost  and  effectiveness,  then  the  problem  becomes  a  multicriteria 
one.  This  problem  can  be  solved  using  the  Lagrangian  relaxation  technique  by  for¬ 
mulating  this  problem  as  one  of  maximising  effectiveness  given  a  certain  budgetary 
limit  on  cost.  Details  of  this  technique  are  given  in  section  D  of  this  chapter.  The 
“best”  path  can  also  be  solved  using  linear  programming  formulation  of  the  network 
flow  problem.  However,  the  dynamic  programming  formulation  has  more  intuitive 
appeal  and  lends  itself  easily  to  incorporation  of  lead-time,  discounting  and  other 
considerations. 
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The  aging  process  is  modeled  by  discounting  performance  and  cost  character¬ 
istics  over  time,  based  on  empirical  parametric  relationships.  What  these  charac¬ 
teristics  will  be  for  a  given  alternative  acquired  new  in  a  certain  year  is  projected 
based  on  technological  trends.  Forecasts  of  inflation  and  foreign  exchange  rates  are 
also  factored  into  the  projected  acquisition,  maintenance  and  other  costs.  All  out¬ 
lays  are  converted  to  current  year  dollars  based  on  forecasted  interest  rates.  The 
projected  rates  may  vary  from  year  to  year  but  a  steady  trend  is  assumed  beyond  an 
arbitrary  forecast  horizon.  Performance  requirement  levels  may  also  be  stipulated 
for  each  year  based  on  analysis  of  politico-military  trends.  A  unique  set  of  life-cycle 
cost  and  effectiveness  measures  is  thus  associated  with  the  network  for  each  scenario 
expressing  a  given  set  of  expectations  concerning  relevant  future  trends.  The  mini¬ 
mum  cost  solution  and  efficiency  boundary  for  each  scenario  can  then  be  obtained 
by  applying  dynamic  programming  and  Lagrangian  relaxation  to  the  corresponding 
network.  The  efficiency  boundary  expresses  the  upper  bound  on  effectiveness  over 
the  entire  spectrum  of  life-cycle  cost  spanned  in  a  given  scenario  (Figure  3.2).  The 
point  on  this  boundary  that  yields  the  highest  effectiveness- to-cost  ratio  is  the  so¬ 
lution  that  gives  the  most  bang  for  the  buck.  As  the  boundary  is  piecewise  linear, 
all  points  along  a  line  segment  may  qualify  if  extrapolation  brings  the  line  through 
the  origin.  Otherwise,  a  Lagrangian  solution  point  at  the  intersection  of  two  line 
segments  will  qualify. 

The  reader  may  wonder  why  a  three  stage  process  is  adopted  in  this  model. 
1  he  model  itself  does  not  restrict  the  user  to  any  number.  Rather,  the  matter  is 
one  of  judgment  concerning  the  credibility  of  long-term  forecasts.  Although  it  is 
desirable  to  optimise  over  as  long  a  stream  of  systems  as  possible,  it  is  also  futile 
to  try  to  visualise  the  distant  future.  A  two-stage  model  comprising  the  present 
system  and  its  immediate  successor  is  satisfactory  as  far  as  a  reasonable  optimisation 
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horizon  goes.  A  three-stage  model  would  be  preferable  as  it  allows  realistic  bounds 
to  be  developed  for  the  immediate  successor’s  economic-technological  life.  The  cost 
of  competing  alternatives  can  therefore  be  more  equitably  compared  since  their 
amortisation  periods  may  be  more  accurately  estimated.  Infinite  horizon  models 
meant  for  static  environment  problems  are  clearly  not  suitable  for  our  purposes. 
Accordingly,  assumptions  of  monotonicity,  contraction,  boundedness,  etc.,  which 
are  required  for  asymptotic  stability  of  such  models  [Henig  1978,  Beckmann  1968, 
Boudarel  1971]  need  not  be  made. 

It  may  be  noted  that  a  long  lead  time  often  precedes  transitions,  be  they 
upgradings  or  replacements.  Lead  time  of  an  alternative  is  defined  to  mean  the 
minimum  number  of  years  by  which  the  initiation  of  acquisition  must  precede  the 
achievement  of  initial  operating  capability  (IOC).  By  setting  to  infinity  the  cost  of 
alternatives  that  cannot  be  introduced  because  of  insufficient  preceding  lead  time, 
we  can  avoid  generating  solutions  contradicting  lead  time  requirements.  This  check 
is  made  not  only  for  transitions  from  the  present  system  to  its  immediate  successor, 
but  also  for  subsequent  ones.  The  latter  requires  the  assumption  that  initiation 
of  a  transition  never  begins  before  the  preceding  one  even  takes  place.  A  path  in 
the  network  is  thus  deemed  infeasible  if  there  is  any  arc  representing  a  system  life 
shorter  than  the  lead  time  required  by  its  successor. 

While  no  one  can  reliably  predict  the  future,  one  ignores  future  possibilities 
and  trends  only  to  one’s  disadvantage.  Taking  calculated  risks  after  examining 
a  broad  range  of  scenarios  is  the  only  realistic  approach  to  handling  uncertainties 
about  the  future.  For  this  purpose,  the  flexibility  of  the  model  presented  is  a  valuable 
asset.  Long-range  projections  can  be  varied  more  than  near-term  ones  to  account 
for  greater  uncertainty  in  these  projections. 
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It  should  be  emphasized  that  the  focus  of  analysis  is  on  the  timing  and  choice 
of  alternatives  for  the  first  transition,  the  optimality  of  which  is  affected  less  by 
inaccuracies  in  long-range  forecasts  than  short-term  ones.  By  updating  the  scenarios 
annually  for  changes  in  the  perception  of  qualified  experts,  rolling  plans  can  be 
developed  for  system  renewal  decisions  to  be  made  optimally  and  in  a  timely  fashion. 
For  this  reason,  the  model  is  named  the  Dynamic  System  Renewal  Planning  Model 
(DSRPM). 
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Figure  3.2  :  Domination  of  Solution  Space  by  Efficiency  Boundary 


C.  ASSUMPTIONS 

For  the  purposes  of  this  thesis,  the  following  assumptions  have  been  used: 

•  A  steady  state  exists  beyond  an  arbitrary  forecast  horizon  where  all  projections 
are  either  constant,  linear  or  harmonic.  This  assumption  expresses  the  lack  of 
detailed  knowledge  about  the  remote  future. 

•  Cost  and  performance  parameters  can  be  measured  deterministically  and  util¬ 
ity  functions  can  be  obtained  explicitly.  (Note:  This  assumption  would  be 
relaxed  in  future  embellishments  incorporating  stochastic  models  and  implicit 
approaches  to  utility  modeling.) 

•  An  upgraded  system  is  a  new  system  with  acquisition  cost  equal  to  the  direct 
upgrading  costs  plus  the  salvage  or  resale  value  of  the  original  system  in  the 
previous  year.  A  new  acquisition  with  no  precedent  is  simply  one  in  which 
the  incumbent  system  is  a  null  or  dummy  system.  With  this  assumption, 
upgrading,  replacement  and  new  acquisition  alternatives  are  integrated  within 
a  common  framework. 

•  Acquisitions  (upgrading  or  replacement)  and  their  associated  outlays  occur 
at  the  beginning  of  a  transition  year.  The  old  system  is  disposed  of  at  the 
end  of  the  year.  This  definitive  assumption  provides  for  proper  costing  of 
the  metamorphic  process  by  which  a  new  system  achieves  initial  operating 
capability  (IOC)  while  the  old  is  prepared  for  disposal.  Disposal  docs  not 
really  occur  in  an  upgrading,  but  the  new  capability  is  still  assumed  to  be 
achieved  at  the  end  of  the  transition  year. 

•  No  transition  is  initiated  before  the  preceding  one  takes  place,  as  was  explained 
in  the  previous  section. 
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One  should  note  the  small  number  and  unrestrictive  nature  of  the  assumptions 
made  in  the  model.  This  is  a  direct  consequence  of  the  versatility  of  the  methodology 
employed.  There  are  other  assumptions  needed  to  model  particular  alternatives  and 
scenarios  such  as  those  pertaining  to  utility  functions,  nature  and  impact  of  trends, 
etc.,  but  these  are  not  generic  to  the  model  discussed. 

D.  METHODOLOGY 

Dynamic  programming  and  Lagrangian  relaxation  are  the  cornerstones  of  the 
model's  integrity  and  the  source  of  much  of  its  power.  This  section  explains  these 
key  techniques  as  applied  to  finite  discrete  decision  processes.  As  methodologies 
concerning  forecasting,  cost  estimation  and  other  components  of  the  model  may 
vary  according  to  the  peculiarities  of  the  application,  these  will  not  be  discussed. 
The  reader  is  referred  to  the  bibliography  for  materials  on  these  aspects  of  the  model. 

1.  Dynamic  Programming 

Borrowing  the  notation  used  by  the  pioneering  developer  of  dynamic  pro¬ 
gramming,  Richard  E.  Bellman,  [Bellman  and  Kalaba,  1965],  we  begin  by  examining 
the  important  notion  of  policy  in  the  context  of  system  renewal.  If  a  decision  c/,  in 
the  ith  stage  results  in  a  transition  from  a  state  y,  to  a  state  y,+ j,  then  an  n-stage 
process  may  be  represented  as 


yi  =  T(y0,d0) 

V2  =  7’(yi,d,) 

yn  =  7  (yn-,,dn_i)  (3.1) 

In  the  3-stage  system  renewal  process  depicted  in  Figure  3.1,  for  instance, 
the  process  of  operating  the  present  system  for  a  year,  then  acquiring  alternative 
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1  to  use  for  two  years  before  acquiring  anew  another  alternative  1  for  a  further 
two  years  use  would  be  represented  by  a  path  Px< 3i 7iI2  of  unbroken  line  arcs  passing 
through  nodes  1,  3,  7  and  12. 

Suppose  we  are  concerned  with  maximizing  a  prescribed  scalar  function 
of  the  states  and  the  decisions: 

■ft(yO,2/l,J/2,-’-,yn,  dQ,  d\ 1  '  ‘  ‘  »  ^n-l  )  (3-2) 

This  so-called  criterion  function,  in  the  case  of  a  system  renewal  process,  may  be 
the  cumulative  expected  return  accrued  from  acquiring  and  operating  a  particular 
succession  of  alternatives,  i.e.,  traversing  a  given  path  in  the  network.  Without  an 
optimal  rule  as  to  how  decisions  should  be  made,  the  only  way  to  obtain  the  max¬ 
imum  of  R  and  its  associated  path  is  by  exhaustive  enumeration  along  all  possible 
paths.  Effort  is  much  reduced  with  the  use  of  rules  in  the  form 

dk  =  J/i ,  •  •  • ,  yk t  do,  d j ,  •  •  • ,  dk-\)  (3.3) 

These  functions  are  called  policy  functions,  or  simply  policies.  A  policy  that  maxi¬ 
mizes  the  criterion  function  R  is  called  an  optimal  policy. 

Let  us  consider,  in  particular,  policies  which  depend  only  on  the  current 
state,  k ,  that  is  policies  of  the  form 


dk=dk(yk)  (3.4) 

This  additional  simplification  is  possible  only  when  the  structure  of  the  criterion 
function  R  has  the  important  property  of  divorcing  the  past  from  the  future.  Ex¬ 
amples  of  criterion  functions  with  this  “separable  criterion"  property  include 
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(3.5) 


k 


£C(i n,d,) 

i= 0 


k 

and  llr(y*^0  (3.6) 

»=n 

The  first  example  is  of  the  most  common  type  of  criterion.  The  usual  criterion  of 
cumulative  cost  or  cumulative  returns  used  in  system  renewal  is  of  this  variety.  The 
reliability  of  series  systems  would  perhaps  come  to  mind  in  the  case  of  the  second 
one.  An  example  of  a  criterion  function  that  does  not  have  the  separability  property 
is  that  of  an  overall  ratio,  that  is, 


¥ -  (3.7) 


i=0 

The  use  of  overall  cost-efFectiveness  ratio  as  a  criterion  in  the  renewal  process  would 
thus  require  decisions  to  be  based  not  only  on  the  current  state  but  also  on  the  path 
leading  to  that  state. 

When  the  criterion  function  used  has  this  vital  separability  property,  the 
optimal  policy  is  characterized  very  simply  by  the  following  principle  of  optimality: 


(yo,yi,--- ,ym,ym+u-- ■  ,yn)  is  optimal  =J>  {y0,  yu  •••,t/m)  is  optimal  (3.8) 

This  intuitive  result  can  be  proved  by  contradiction.  Suppose  that  for  the  (m  -4-  1) 
stage  process,  (t/0,J/i,- •  • , ym)  is  not  optimal.  Specifically  say  (r0, Xj,  •  •  • , zm)  is 
optimal.  Then  (x0,  Zi,  •  •  ■ ,  zm,  ym+\,  •  ■  ■  *  Vn)  is  better  than  (t/o,  t/i ,  •  •  • ,  ym,  ym+u 
-••,t/n).  But  this  contradicts  the  original  premise  that  (yo,yi,  •  •  ■ ,  ym ,  ym+ii  •  ■  • ,  y„) 
is  optimal.  Therefore,  expression  3.8  must  be  true.  Using  the  principle,  one  obtains 
for  the  examples  given  by  expressions  3.5  and  3.6  the  following  optimal  policies: 


fk+\(T(yk,(lk))  =  max[c(yjt,f/fc)  +  fk(yk)}  (3.9) 
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(3.10) 


fk+\{T(yk,dk))  =  max  [c(yfc,  4)/*(j/fc)] 

a* 

where  fi(yi,d,)  is  the  criterion  function  in  state  i  obtained  using  an  optimal  policy. 

Dynamic  programming  hinges  around  the  formulation  of  optimal  policies 
to  enable  solutions  to  be  found  by  systematic  recursion  [Beckmann,  1968].  A  solution 
may  be  given  in  terms  of  the  sequence  of  criterion  functions  {/,(y,)}  or  the  sequence 
of  policy  functions  {(4(y*)}  .  Each  determines  the  other,  although  one  should  note 
that  there  may  be  several  optimal  policies  yielding  the  same  sequence  of  criterion 
functions.  A  policy  may  also  be  satisfied  by  one  or  more  solutions.  Only  when  the 
solution  is  unique  is  the  policy  a  necessary  and  sufficient  condition  for  optimality. 
Otherwise,  it  is  necessary  to  compare  the  solutions  using  additional  considerations 
such  as  risk,  expected  life-span,  etc.,  to  define  the  optimal  solution. 

Being  a  flexible  and  potent  tool,  dynamic  programming  has  many  vari¬ 
ations.  In  stochastic  dynamic  programming,  the  next  state  is  chosen  based  on 
conditional  transition  probabilities  and  the  criterion  function  is  given  as  a  random 
variable,  usually  the  mean  of  possible  rewards  [Ross,  1983].  Variations  with  vector¬ 
valued  criterion  functions  have  been  developed  [Henig,  1978],  but  these  multicriteria 
dynamic  programming  models  are  still  not  as  widely  applied  as  the  stochastic  vari¬ 
ants.  Although  these  multicriteria  models  could  be  used  in  system  renewal  problems 
to  obtain  dominating  (Pareto-efficient)  solution  sets,  it  is  simpler,  more  useful,  nat¬ 
ural  and  elegant  to  use  the  Lagrangian  relaxation  technique  in  such  problems,  as 
shall  be  shown  in  the  following  section. 

2.  Lagrangian  Relaxation  Technique 

Two  conceptual  approaches  used  in  systems  analysis  are  either  to  fix 
effectiveness  and  minimize  cost  or  fix  the  budget  and  maximize  effectiveness  [Fisher, 
1970].  Suppose  we  choose  (he  latter  and  formulate  the  system  renewal  problem  as 
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n 

max 

do,  -,dn  )  =  0 


(3.11) 


subject  to 

jrC{yi,di)  <  Budget  (3.12) 

1=0 

where  £(j/;,d, )  and  C(y;,<L)  are  the  life-cycle  effectiveness  and  cost  respectively 
of  making  a  decision  d{  in  state  j/j.  Units  for  the  former  are  utiles,  while  that  of 
costs  are  dollars.  In  the  system  renewal  process  shown  in  Figure  3.1,  for  example, 
if  the  current  state  is  represented  by  node  3,  and  the  life-cycle  effectiveness  and 
cost  of  acquiring  alternative  1  and  operating  it  for  two  years  are  2,100  utiles  and 
$7m  respectively,  then  £(3,2)  =  2,100  utiles  and  C(3,2)  =  $7m.  The  constraint 
in  expression  3.12  may  be  relaxed  by  moving  it  to  the  objective  function  with  the 
inclusion  of  a  Lagrangian  multiplier  A,  thus 


na5  Ie 

"  ,dn  v  t=0 


III  X  r£(y„^)-A 


^2C(y{,dt)  -  Budget 


Li=0 


(3.13) 


Noting  that  expressions  3.5  and  3.13  are  of  the  same  form,  we  know  that  the  optimal 
policy  for  maximizing  this  criterion  function  is  given  by  the  principle  of  optimality 


as: 

A+i (T(yk,dk))  =  max  {[E(yk,dk)  -  \C(yk,dk)}  +  fk(yk)}  (3.14) 

a* 

Therefore,  for  any  fixed  A,  conventional  dynamic  programming  may  be  used  to 
obtain  the  optimal  solution  to  expression  3.13.  Lagrangian  relaxation  is  thus  a 
technique  for  converting  constrained  optimisation  problems  into  unconstrained  max¬ 
imization  (or  minimization)  problems  (Everett,  1963].  It  may  also  be  interpreted  as 
a  technique  to  solve  the  dual  of  the  original  problem  given  by  Equations  3.11  and 
3.12,  as  will  be  apparent  in  the  following  discussion. 

To  consider  the  effect  of  varying  A,  let  us  define  the  Lagrangian  function 


L( A)  to  be  equal  to  the  objective  function  given  in  expression  3.13.  T  hat  is, 


1(A)  =  max  d,)  -  A 

do,",dn  (,=0 


5Z  C{yt,d,)  -  Budget 


1=0 


(3.15) 


This  piecewise  linear  function  is  illustrated  in  Figure  3.3. 


Figure  3.3:  Lagrangian  Function  L(A) 


When  A  =  0,  the  budgetary  constraint  is  ignored  and  we  obtain  the 
maximum  effectiveness  solution  to  the  (n  +  1)  system  renewal  problem  given  in 
expression  3.13.  Denoting  the  cumulative  life-cycle  effectiveness  and  cost  of  this 
solution  as  LCEq  and  LCC’o  respectively,  we  note  that  L( 0)  =  LCEq  and  that  the 
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slope  at  this  vertical  axis  intercept  is  —(LCCo  —  Budget).  For  marginal  increases  in 
A  beyond  zero, 


L{ A)  =  LCEq  -  \{LCC0  -  Budget)  (3.16) 

which  is  the  equation  of  the  leftmost  line  segment.  When  A  reaches  a  threshold 
value,  the  solution  changes  to  one  with  smaller  cumulative  life-cycle  effectiveness 
and  cost.  The  next  line  segment  will  thus  have  a  lower  intercept  value  and  an 
increased  slope.  This  therefore  implies  that  the  Lagrangian  function  will  always  be 
convex  for  system  renewal  problems  formulated  as  in  expression  3.13. 

The  Lagrangian  multiplier  A  has  a  physical  interpretation,  namely  that 
of  the  relative  value  of  cost  and  effectiveness,  i.e.,  shadow  or  dual  price.  Its  unit  is 
utiles  per  dollar.  As  A  is  increased  from  zero  to  infinity,  the  Lagrangian  function 
steps  through  the  entire  solution  space,  from  the  maximum  effectiveness  solution  to 
the  minimum  cost  solution.  The  latter  corresponds  to  the  rightmost  line  segment 
which  extends  to  infinity.  What  then  is  the  optimal  value  of  A?  To  see  this,  we  first 
note  that  for  A  >  0,  L( A)  is  an  upper  bound  on  the  optimal  value  F*  of  the  original 
objective  function  given  in  expression  (3.11).  This  is  because  nunL(A)  is  the  dual 
of  the  original  problem.  The  minima  of  L{ A)  therefore  yields  F*.  That  is, 

mmI(A)  =  I(r)  =  r  (3.17) 

This  is  akin  to  seeking  the  constrained  stationary  points  in  analog  applications  of 
Lagrange’s  method  to  differentiable  functions. 

The  minima  of  L( A)  do  not  necessarily  exist.  In  a  commonly  observed 
case,  LCCo  is  higher  than  the  budget  and  the  slope  at  the  vertical  intercept  is 
negative.  As  L{ A)  is  convex,  A*  is  always  finite  and  non-zero  in  such  cases.  If, 
however,  LCCo  is  lower  than  the  budget,  that  is  to  say  the  maximum  effectiveness 
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solution  is  affordable,  then  the  slope  at  the  intercept  will  be  positive  and  we  obtain 
A*  =  0.  Should  LCC0  be  higher  than  the  budget  and  there  are  no  affordable 
solutions,  then  there  will  be  no  line  segments  of  non-negative  slopes,  no  minima  and 
hence  no  feasible  solution. 

The  value  of  A*  can  be  found  most  expeditiously  and  accurately  by  the 
supporting  plane  method  illustrated  in  Figure  3.4.  Given  the  upper  and  lower 
bounds  on  A*,  the  method  involves  the  extrapolation  of  the  line  segments  corre¬ 
sponding  to  these  bounds  to  obtain  the  intersection  value  A,.  The  intercept  and 
slope  of  these  line  segments  are  given  by  the  cumulative  life-cycle  effectiveness  and 
cost  respectively  of  the  solution  obtained  by  solving  expression  (3.13)  for  the  bound¬ 
ing  values  of  A.  The  lower  or  upper  bound  is  replaced  by  A,  if  the  slope  of  the  line 
segment  corresponding  to  A,  if  the  slope  of  the  line  segment  corresponding  to  A^  is 
negative  or  non-negative  respectively.  This  process  is  iterated  until  A,  equals  one  of 
the  bounds.  This  equality  occurs  because  L{ A)  is  piecewise  linear  and  convex,  which 
enables  A*  to  be  obtained  precisely.  The  starting  lower  bound  is  set  to  be  zero.  Any 
suitably  large  value  of  A  that  yields  a  line  segment  with  non-negative  slope  may  be 

.  .  .  H 

used  as  the  initial  upper  bound  A^^j.,  where 

=  No.  of  nodes  x  LC Emar  4-  LCCmin  (3.18) 

Experience  of  the  author  with  similar  problems  shows  this  supporting  plane  method 
to  be  faster  and  more  accurate  than  the  commonly  used  bisection  method,  and  A* 
may  be  found  in  only  a  few  iterations  even  for  very  large  networks  with  thousands 
of  nodes  and  arcs. 

For  each  scenario,  using  the  methodology  described,  we  are  able  to  obtain 
the  solution  points  marking  the  efficiency  boundary  shown  in  Figure  3.2.  This 
is  begun  by  artificially  setting  a  sufficiently  high  budget  to  obtain  the  maximum 
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0  ^2  A3A4  Aj  Amu* 


Figure  3.4:  Iterations  for  Derivation  of  A* 

effectiveness  solution,  which  is  marked  1  in  Figure  3.5.  The  corresponding  value 
of  A*  is  obviously  zero.  The  budget  is  then  set  to  just  under  the  value  of  the 
cumulative  life-cycle  cost  C\  of  the  maximum  effectiveness  solution.  This  produces 
the  next  solution  point,  marked  2  in  the  diagram,  and  the  process  is  repeated  until 
we  reach  the  minimum  cost  solution,  which  is  the  sixth  solution.  No  solutions  will  be 
obtained  for  budgets  less  than  Ce-  The  optimal  Lagrangian  multiplier  A*  obtained  in 
each  case  determines  the  slope  of  the  line  segment  to  the  right  of  the  corresponding 
solution  point.  It  may  be  noted  that  the  efficiency  boundary  consists  of  bounding 
hyperplanes  which  envelope  the  solution  space  but  do  not  pass  through  all  the 
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Pareto-efficient  points.  There  are,  for  instance,  three  such  points  denoted  by  circles 
for  budgets  between  C3  and  C4  which  will  not  be  discovered  by  the  Lagrangian 
relaxation  method.  The  basic  cause  of  these  inaccessible  regions  (called  gaps)  is 
non-concavity  in  the  relationship  between  the  optimal  effectiveness  and  cost.  These 
gaps  will  not  exist  if  the  law  of  diminishing  marginal  utility  holds  strictly.  H. 
Nakayama,  et  al . ,  [Nakayama,  1975],  has  a  good  survey  of  methods  to  deal  with  this 
characteristic  of  the  Lagrangian  relaxation  method.  The  existence  of  these  gaps  does 
not  create  a  problem  as  far  as  the  application  of  the  method  in  systems  analysis  is 
concerned.  This  is  because  the  solution  offering  a  maximum  overall  effectivenss-to- 
cost  ratio  is  always  a  Lagrangian  solution  point.  Also,  a  straightforward  application 
of  dynamic  programming  will  produce  the  minimum  cost  solution  even  if  it  does  not 
lie  on  the  efficiency  boundary.  Note  that  as  system  lives  are  modeled  explicitly  for 
economic,  technological  and  politico-military  changes,  the  net  cost  of  alternatives 
can  be  amortized  over  the  proper  period  to  reveal  the  true  minimum  cost  solution 
as  defined  by  amortized  annual  cost,  as  will  be  amplified  in  the  following  paragraph. 
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Figure  3.5  :  Derivation  of  Efficiency  Boundary  hy 
Lagrangian  Relaxation  Technique 
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We  have  seen  how  characterization  of  the  systems  renewal  problem  as  a 
constrained  optimization  problem  (in  which  one  seeks  to  maximize  effectiveness  for  a 
given  budget)  enables  us  to  obtain  the  efficiency  boundary  corresponding  to  a  given 
scenario.  A  band  of  these  boundaries  is  found  for  each  set  of  scenarios  investigated. 
If  probabilities  are  attached  to  scenarios,  our  “expected”  efficiency  boundary  and 
its  associated  optimal  solutions  can  be  obtained.  Furthermore,  bearing  in  mind  that 
the  cost  of  an  arc  in  the  network  is  set  to  infinity  if  there  is  insufficient  lead  time 
preceding  it  or  if  any  performance  parameter  deteriorates  below  its  requirement  level, 
we  know  that  any  solution  with  finite  cost  is  feasible  from  the  effectiveness  point  of 
view.  The  minimum  cost,  in  current  year  dollars,  for  the  various  feasible  solution 
paths,  as  well  as  the  associated  period  of  use,  can  then  be  translated  into  amortized 
costs  flows  for  given  interest  rate  assumptions.  On  this  basis,  we  can  find  the  solution 
that  allows  us  to  minimize  equivalent  annual  cost  while  meeting  a  given  set  of  lead 
time  and  performance  requirements.  Therefore,  both  the  fixed  budget  and  the  fixed 
effectiveness  approaches  are  embodied  within  the  model,  delivering  the  maximum 
overall  efFectiveness-to-cost  ratio  solution  and  the  minimum  cost  satisficing  solution 
respectively.  Furthermore,  it  should  be  noted  that  A*  has  the  interpretation  of 
shadow  or  dual  prices  in  economic  analysis.  Therefore,  efficiency  boundaries  of 
systems  with  similar  MOEs  can  be  meaningfully  compared.  Also,  choice  of  system 
mix  for  various  capabilities  may  be  guided  by  the  fact  that  shadow  prices  should 
be  roughly  the  same  for  all  the  systems  in  an  optimal  mix.  The  combined  use 
of  dynamic  programming  and  the  Lagrangian  relaxation  technique  thus  provide  us 
with  a  conceptually  simple,  vet  potent  and  useful  methodology  that  is  completely 
congruent  with  and  derivable  from  the  tenets  of  systems/economic  analysis. 
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IV.  IMPLEMENTATION 

A.  INTRODUCTION 

The  dynamic  systems  renewal  planning  model  (DSRPM)  presented  in  the  pre¬ 
vious  chapter  provides  a  formal  basis  by  which  diverse  economic,  technological  and 
politico-military  concerns  may  be  integrated  within  a  common  framework  (Figure 
4.1)  to  make  for  more  rational  management  of  defense  assets.  Credible  and  valid 
results,  however,  can  be  obtained  only  with  the  active  participation  of  relevant  logis¬ 
tics,  engineering,  finance,  intelligence  and  operational  agencies  in  the  development 
of  proper  effectiveness  and  cost  sub-models  as  well  as  realistic  scenarios.  High-level 
management  support  is  thus  essential  for  the  successful  implementation  of  the  model. 
To  foster  a  better  understanding  of  the  model  and  to  demonstrate  its  plausibility 
and  usefulness,  a  small  prototype  was  developed  based  on  a  simple  hypothetical 
weapon  system.  This  prototype  is  documented  in  the  next  chapter  along  with  case- 
studies.  Although  the  addition  of  selected  successful  applications  would  give  greater 
credence  to  the  model,  the  limited  time  available  for  work  on  this  thesis  precluded 
full-scale  implementation  for  any  real  systems.  Nevertheless,  efTort  shall  be  made 
to  outline  conceptually  how  implementation  is  envisaged.  Only  aspects  of  cost  es¬ 
timation,  effectiveness  modelling  and  requirement  formulation  peculiar  to  DSRPM 
is  discussed.  For  a  full  discussion,  the  reader  is  referred  to  relevant  materials  in  the 
bibliography. 
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Figure  4.1:  System  Renewal  Planning  Framework 
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B.  CONCEPTUAL  FRAMEWORK 

Implementation  of  the  model  is  best  explained  in  terms  of  its  constituent  el¬ 
ements  and  their  relationships  to  the  processes  involved  in  defense  planning.  As 
depicted  in  Figure  4.1,  DSRPM  has  five  elements,  all  of  which  are  part  of  the 
force/resource  planning  process.  A  model  with  these  five  elements  is  developed 
for  each  entity  of  force  planning,  which  may  be  weapon  systems,  task  forces,  bat¬ 
tle  groups,  etc.  In  renewal  planning  parlance,  all  these  entities  are  simply  called 
systems,  and  candidate  systems  are  called  alternatives. 

The  cost  estimation  element  projects  the  life-cycle  costs  for  both  the  present 
system  as  well  as  future  alternatives.  The  effectiveness  modelling  element  concerns 
the  appraisal  of  the  value  of  a  system  as  a  function  of  its  measures  of  effectiveness. 
The  deterioration  modelling  element  accounts  for  the  change  in  effectiveness  and 
cost  over  the  life  of  a  system  caused  by  wear-out,  operational  attrition  and  enemy 
response.  The  other  two  elements  are  self-explanatory.  All  of  these  elements  and 
their  implementation  are  elaborated  in  sections  C,  D,  and  E.  Only  deterioration 
modelling  and  the  selection  and  scheduling  of  programs  are  organic  to  the  DSRPM. 
The  other  elements  are  on-going  force/resource  planning  processes  adapted  with 
future  projections  for  use  in  the  DSRPM. 

Formulation  of  future  alternatives  is  a  part  of  the  force/resource  planning  pro¬ 
cess  that  lies  wholly  outside  the  province  of  the  DSRPM.  It  involves  definition  of 
the  physical  and  functional  forms  of  alternatives  identified  in  the  strategy  formu¬ 
lation  process  and  is  shaped  by  technological  outlook.  Alternatives  encompass  the 
following: 

1.  Conversion  of  incumbent  systems  to  serve  different  roles 

2.  Upgrading  of  incumbent  systems  to  serve  expanded  or  higher  roles 
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3.  Modernisation  of  incumbent  systems  to  better  serve  their  assigned  roles 

4.  Replacement  of  incumbent  systems  with  newer  and  probably  better  alterna¬ 
tives 

5.  Activation  of  reserve  systems  for  active  deployment  (incumbent  systems  are 
the  ones  in  dormant  or  preserved  status) 

6.  Reduction  of  number  of  incumbent  systems  through  disposal/sale 

7.  Retirement  of  incumbent  systems  into  the  reserves 

8.  Acquisition  of  new  (no  precedent)  systems  to  meet  qualitative  and/or  quanti¬ 
tative  deficiencies  or  new  (threat)  requirements. 

No  present  system  exists,  i.e.,  system  A  in  the  network  shown  in  Figure  3.1 
does  not  exist,  in  the  case  of  new  acquisitions. 

C.  EFFECTIVENESS  MODELING 

The  effectiveness  modeling  adopted  in  the  DSRPM  is  completely  congruent 
with  the  way  military  systems  are  selected  and  managed.  Each  system  has  one  or 
more  missions  that  represent  the  raison  d'etre  of  that  system.  These  missions  are 
usually  phrased  in  terms  that  are  too  broad  for  concrete  guidance.  What  it  takes  to 
carry  out  the  mission(s)  are  stipulated  clearly  as  required  operational  (capabilities) 
requirements.  These  are  usually  lower  bounds  set  on  key  performance  parameters 
such  as  speed,  payload,  survivability,  operational  range,  etc.  The  effectiveness  of 
candidate  systems  is  evaluated  based  on  those  key  parameters  and  some  criterion. 
The  decision  to  choose  a  particular  system  need  not  be  based  cxclnsi vely  on  objec¬ 
tive  parameters,  and  the  criterion  used  is  not  necessarily  a  formal  one.  I  he  views 
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of  qualified  experts  and  decision  makers  are  also  important.  In  fact,  the  politi¬ 
cal  process  leading  to  the  decision  may  even  be  colored  by  elements  that  detract 
from  the  military  contribution  to  national  security,  such  as  amicable  relations  with 
another  state.  However,  assuming  rational  arguments  hold  sway,  an  effectiveness 
model  representing  the  embodiment  of  the  values,  preferences  and  possibly  biases 
of  the  decision  makers,  tempered  with  experience  gained  from  performance  audits, 
can  be  constructed.  Justification  papers,  test  and  evaluation  plans,  trial  reports, 
etc.  are  good  sources  of  information  needed  to  construct  an  effectiveness  model. 
Combat  models  and  wargames  also  provide  valuable  input.  The  aim  of  effectiveness 
modelling  for  this  model  is  to  describe  value  judgements  of  the  decision-makers  and 
not  to  predict  battle  outcomes  (although  predictions  do  tend  to  afTect  judgement). 
Therefore,  the  same  effectiveness  model  used  for  evaluation  of  alternatives  in  the 
acquisition  phase,  and  possibly  altered  by  performance  audits  during  deployment, 
is  incorporated  in  the  DSRPM.  The  only  difference  in  the  DSRPM  case  is  that 
prospective  alternatives  are  evaluated  alongside  the  incumbent,  and  performance 
parameters  are  projected  both  for  future  generations  of  each  alternative  as  well  as 
through  the  years  of  use  after  an  alternative  is  acquired.  Figure  4.2  shows  how  the 
jth  performance  parameter  of  an  alternative  as  acquired  in  year  y  may  change  over 
the  years  of  usage  t. 

As  an  example,  let  us  suppose  that  a  weapon  system  was  evaluated  based  on 
an  effectiveness  model  with  the  following  additive  structure: 


TE„  =  ^  Hjf/j  {A'Jia}  provided  XJ[a  >  for  Vj 

j= i 

=  0  otherwise 


(4.1) 
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where 


•  TEa  is  the  total  effectiveness  of  alternative  a. 

•  Wj  is  the  weight  of  the  jth  performance  parameter  expressing  its  importance 
relative  to  the  other  parameters. 

•  XjA  is  the  value  of  the  jth  parameter  for  alternative  a. 

•  Uj  {Xj^a}  is  the  utility  function  that  transforms  XjA  into  a  utility  value. 

•  2Lj  is  the  minimum  level  required  of  the  jth  parameter  (operational  require¬ 
ment). 

In  this  case,  the  effectiveness  sub  model  used  in  the  application  of  DSRPM  to  the 
system  will  be: 

-  f;  WjUj  (4.2) 

j= i 

where 

•  TEa'<y(t )  is  the  total  effectiveness  of  alternative  a'  (not  necessarily  the  same  as 
a)  if  acquired  new  in  year  y  during  the  tth  year  of  deployment. 

•  Xjta\y  is  the  value  of  the  jlh  performance  parameter  for  alternative  a'  if  ac¬ 
quired  in  year  y. 

•  daj(t)  is  the  deterioration  function  for  the  jth  parameter  of  alternative  a'  in 
the  tlh  year  of  deployment. 

The  alternatives  a  and  a'  need  not  be  the  same.  The  latter  typically  comprises 
not  only  the  replacement  alternative,  but  also  upgrading  options.  Note  that  the 
structure  of  the  effectiveness  model,  weights  and  utility  functions  are  retained  in 


the  DSRPM.  TEa',y(t)  is  the  effectiveness  associated  with  all  arcs  representing  al¬ 
ternative  a',  leading  out  of  nodes  of  year  y,  with  life-spans  of  t  years  in  the  network 
shown  in  Figure  3.1.  If  Xa',j(t)  <  X.j,v+t  lor  any  ji  i-e.,  if  any  parameter  fails  to 
meet  the  requirement  level  stipulated  for  the  year  of  operation  (y-M),  then  the  cost 
of  the  corresponding  arc  is  set  to  infinity. 

The  effectiveness  model  used  here  is  similar  to  that  of  the  TASCOM  method¬ 
ology  developed  by  The  Analytic  Sciences  Corporation  for  OSD  over  several  years 
in  the  late  1970’s  and  the  early  1980’s.  The  TASCOM  model  provides  a  first-order 
indication  of  an  aircraft  performance  potential  by  computing  a  score  based  on  at¬ 
tributes  such  as  range,  speed,  maneuverability,  payload,  basing  mode,  target  acqui¬ 
sition  and  fire  control  capabilities.  These  factors  are  normalised  against  a  consistent 
standard  or  assigned  a  graduated  evaluation  score.  The  relative  weights  are  based 
on  the  consensus  judgment  of  experienced  pilots  as  well  as  military  analysts.  The 
final  performance  are  scores  based  on  weighted  calculations.  The  reader  is  referred 
to  Hildebrandt,  et  al.,  1986,  and  Congressional  Budget  Office,  “ Tactical  Combat 
Forces  of  the  USAF,  Issues  and  Alternatives,"  April  1985,  for  further  details  on  the 
TASCOM  methodology. 

The  values  of  the  parameters  for  future  versions  of  each  alternative  are  pro¬ 
jected  as  Xj,a',y.  These  are  discounted  (for  deterioration)  through  the  years  of  usage 
by  the  function  d3{t).  Technical  expertise  is  required  for  such  forecasting.  Require¬ 
ment  levels  X-j,y+t  need  not  be  static  but  may  be  varied  through  the  future  years 
in  accordance  with  projected  threats.  Figure  4.2  illustrates  how  this  requirement 
level  may  change  with  time.  Projections  are  required  only  for  the  span  of  the  un¬ 
derlying  network  (Figure  3.1)  and  a  steady  state  should  be  assumed  beyond  the 
forecast  horizon.  These  projections  form  the  effectiveness  aspect  of  DSRPM  scenar¬ 
ios.  Examples  of  U}  {  }  and  r/a<;(  )  arc  shown  in  Figures  4.3  and  4.4  respectively. 
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The  utility  function  Uj  {  }  is  the  instrument  that  allows  us  to  express  non-linear 
relationships  such  as  that  of  diminishing  marginal  utility.  Adjusting  it  regularly  for 
the  effect  of  changes  in  enemy  force  structure,  size,  tactics,  strategy,  etc.,  also  serves 
to  account  for  the  interactive  nature  of  defense  “posturing”.  The  example  shown 
here  is  of  explicit  utility  functions  and  single  mission  systems.  Implicit  approaches 
may  be  used  if  utility  functions  can  be  specified  explicitly  only  with  considerable 
difficulty.  Multiple  mission  systems  may  be  modelled  by  combining  utility  values 
using  weights  or  other  multi-criteria  modelling  schemes. 
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Note  :  Utility  functions  need  not  be  continuous 


?igure  4.3  :  Examples  of  Utility  Functions 
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Figure  4.4:  Examples  of  Deterioration  Functions 

D.  COST  ESTIMATION 

Resource  requirements  are  evaluated  for  both  incumbent  systems  and  prospec- 
tive  alternatives.  Cost  models  for  the  former  deal  with  the  projection  of  ownership 
costs  based  on  actual  cost  experience  as  well  as  forecasts  of  disposal  cost  or  salvage 
value.  Those  models  developed  for  future  alternatives  are  either  statistical  ones 
based  on  the  physical  and  functional  characteristics  of  these  alternatives  (Figure 
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4.1)  or  analogy  types  based  on  comparison  of  similar  systems  and  judgement  con¬ 
cerning  complexity  factors.  Procurement  costs  are  included  in  cost  estimation  of 
alternatives  but  are  considered  sunk  costs  in  the  case  of  incumbent  systems. 

Cost  models  used  in  DSRPM  differ  from  the  ones  typically  used  in  evaluations 
in  that  various  types  of  costs  such  as  procurement  cost,  operating  cost,  etc.,  are 
estimated  not  just  for  one  possible  life-cycle  but  for  multiple  possible  ones  in  the 
future.  The  projection  of  costs  needed  to  adapt  cost  models  for  use  in  DSRPM  is 
usually  accomplished  using  time-series  or  regression  analysis.  Adjustments  may  be 
made  if  significant  and  relevant  technological  breakthroughs  are  anticipated.  For  ex¬ 
ample,  if  major  advances  in  certain  automation  technologies  are  expected  five  years 
from  now,  then  manpower  costs  may  be  adjusted  accordingly.  Such  adjustments 
should  be  made  conservatively  and  in  a  not  too  abrupt  manner.  Figure  4.5  shows 
an  example  of  the  cost  parameter  forecasts  that  are  used  in  the  model. 

To  examine  how  costs  may  be  projected  into  the  future,  let  us  consider  the 
following  simplified  li  cycle  cost  model: 

LCCa,y{t)  =  ACa,y  +  +  OCay^{t~l)  +  DCa,y(i)  (4.3) 

DCa,y(t)  =  Ua,y  -  IVW  (4.4) 

where: 

•  LCCa<y{t)  is  the  life-cycle  cost  of  alternative  a  if  procured  in  year  y  and  used 
for  i  years 

•  ACa>y  is  the*  acquisition  cost  of  alternative  a  if  procured  in  year  y 

•  MCa,y  is  the  maintenance  cost  of  alternative  a  in  the  first  year  if  procured  in 


Figure  4.5:  Examples  of  Cost  Parameter  Forecast 


•  0Ca,v  is  the  operating  cost  of  alternative  a  in  the  first  year  if  procured  in  year 

y 


•  DCa:y(t)  is  the  disposal  cost  of  alternative  a  if  procured  in  year  y  and  used  for 
t  years  (Note:  This  may  be  negative  if  asset  has  resale  value.) 

•  Ua,y  is  the  upper  bound  of  DCa,y(t) 

•  Va<y  is  a  proportionality  constant  dictated  by  the  resale  market 

•  7a,yt  and  7a,y  arc  exponential  constants 
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If  inflation  and  technological  change  have  been  (and  presumably  will  be)  gradual, 
the  individual  costs  may  be  forecasted  using  non-linear  regression  models  such  as 
the  following: 

ACa,y  =  PACa{l+rACy  +  ev  (4.5) 

where: 

•  PACa  is  the  present  acquisition  cost  of  alternative  a 

•  r^c  is  the  compound  inflation-technological  improvement  rate  for  acquisition 
cost 

•  e„  is  the  error  term 

Solutions  are  obtained  using  either  the  maximum  likelihood  approach  or  Zellner’s 
method  [Judge,  et  al.,  1988]  based  on  historical  data.  If  major  breakthroughs  are 
expected  in  the  production  of  alternative  a,  the  forecasted  figures  for  ACa<v  may 
be  adjusted  to  some  extent.  Anticipated  foreign  exchange  rates  may  also  be  fac¬ 
tored  into  cost  projections  if  they  are  relevant  to  a  system  renewal  problem.  The 
cost  aspect  of  DSRPM  is  manifested  in  these  projections.  There  are  many  differ¬ 
ent  methods  of  forecasting  both  cost  and  effectiveness  [Ayres,  19C9,  Stewart  and 
Wyskida,  1987,  and  Ascher  and  Overhold,  1988],  the  choice  of  which  depends  very 
much  on  the  degree  of  understanding  one  has  of  the  dynamics  of  the  change  process. 

E.  OTHER  ELEMENTS 

The  deterioration  function  dj(l)  introduced  in  the  section  on  effectiveness  mod¬ 
elling  is  an  example  of  how  the  change  in  performance  parameters  over  the  life  of  a 
system  may  be  modelled.  Audit  programs  such  as  OT&E,  inspections  and  exercises 
as  well  as  technical  trials  and  laboratory  tests  provide  the  data  needed  to  derive 
these  functions.  To  sec  how  deterioration  may  be  incorporated  in  a  cost  model,  let 
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us  consider  the  following  maintenance  cost  model: 


=  cMot0 
dt 


(4.6) 


where: 

•  M  is  the  cumulative  maintenance  cost  incurred 

•  t  is  the  age  of  the  system 

•  c  is  a  proportionality  constant 

•  a  and  0  are  exponentiality  constants 

This  model  links  the  rate  of  change  in  maintenance  cost  to  the  age  t  as  well  as  past 
maintenance  effort  measured  in  terms  of  cumulative  cost.  Given  an  initial  condition 
of  M  —  0  at  t  ~  0,  we  have: 


M  =  Ar 


where: 


1+0 

(l-o) 


A  = 


(4.7) 

(4.8) 

(4.9) 


a  ^  1 


Hence  the  rate  of  increase  in  maintenance  cost  is  given  by: 


(4.10) 


(4.11) 


The  parameters  A  and  7  are  obtained  by  using  regression  analysis  on  historical 
maintenance  cost  date.  A  is  the  maintenance  cost  incurred  in  the  first  year  of  use. 
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As  7  is  always  greater  than  1,  the  rate  of  increase  in  maintenance  cost  ~  will 
always  be  increasing  over  time.  The  impact  of  inflation  and  technological  advances 
in  maintenance  cost  is  modelled  by  varying  both  A  and  7  for  future  alternatives.  A 
reduction  in  operating  cost  for  future  generations  may  also  be  modelled  in  a  similar 
fashion. 

Operational  requirements  for  present  systems  are  projected  by  defining  the 
minimum  requirement  level  in  each  year  for  various  performance  parameters  of  a 
system.  These  are  denoted  as  2Lj<y+t  in  section  B.  For  instance,  shifts  in  enemy  beach 
defence  tactics  may  result  in  significant  increases  in  the  speed  required  of  amphibious 
assault  craft.  Requirement  projections  form  an  integral  part  of  a  DSRPM  scenario 
and  are  driven  mainly  by  the  net  assessment  process  (Figure  4.1). 

Force  levels  decline  even  in  peacetime  with  attrition  due  to  accidents,  wear 
and,  in  the  case  of  human  systems,  retirement  and  resignation.  If  the  number  of  an 
alternative  a'  acquired  in  year  y  is  taken  to  be  the  performance  parameter 
the  decline  due  to  attrition  can  be  modelled  by  the  deterioration  function 
The  effectiveness  of  various  force  levels  can  then  be  modelled  using  the  associated 
utility  function  Uj  {},  and  the  force  level  requirement  in  the  year  ( y  +  t)  is  simply 
2Lj,y+t •  Treating  the  number  acquired  as  a  performance  parameter  in  this  way 
allows  us  to  check  the  adequacy  of  the  size  of  an  acquisition.  The  attractiveness  of 
a  phased  build-up  versus  the  economies  of  a  bulk  buy  may  be  studied  by  defining 
alternatives  representing  different  systems  and  force  levels.  DSRPM  produces  an 
acquisition  program  that  meets  projected  force  level  requirements  by  selecting  and 
scheduling  these  alternatives.  Figure  4.6  illustrates  the  kind  of  solution  that  may  be 
generated  by  DSRPM  in  such  studies.  Force  mix  can  be  optimised  by  ensuring  that 
the  lagrangian  multiplier  associated  with  the  levels  of  the  various  force  components 
are  similar.  In  this  way,  force  level  and  mix  issues  can  be  dealt  with  by  the  model. 


Figure  4.6:  Phased  Acquisition  Solution  Example 

The  last  element  we  look  at  concerns  the  selection  and  scheduling  of  alterna¬ 
tives  for  implementation.  The  inputs  for  these  are  the  cost  and  effectiveness  models 
discussed  in  sections  B  and  C  of  this  chapter  as  well  as  budgetary,  temporal,  legal 
and  other  constraints.  The  combined  dynamic  programming  -  Lagrangian  relax¬ 
ation  discussed  in  Chapter  3  is  applied  to  these  inputs  to  produce  solutions  for  the 
fixed  budget  and  fixed  effectiveness  approaches.  Constraints  in  the  form  of  quanti¬ 
tative  rules  can  be  incorporated  within  DSRPM.  An  example  is  that  of  lead  time 
requirements  for  various  alternatives  which  is  implemented  in  the  prototype  that 
was  developed.  The  model  is  exercised  by  iterating  ‘  for  various  scenarios  to  make 
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explicit  uncertainties  faced  by  the  decision  maker.  Strategic  considerations  with 
the  formal  quantitative  results  from  DSRPM  are  then  the  basis  of  system  renewal 
decisions  which  are  the  output  of  the  force/ resource  planning  process. 

DSRPM  is  easily  implemented  in  phases.  The  model  could  first  be  applied  to 
a  couple  of  incumbent  systems  for  which  systems  renewal  decisions  are  expected  in 
the  next  one  or  two  years.  This  would  allow  enough  time  for  a  full  implementation 
of  the  model  for  these  systems.  This  phase  serves  to  validate  the  model  with  a 
real  application  and  the  experience  gained  would  guide  implementation  in  the  next 
phase.  The  model  should  however  be  applied  as  early  as  possible  in  a  system’s 
life-cycle  to  provide  enough  time  for  the  cost  and  effectiveness  models  to  be  refined 
with  inputs  from  operational  trials,  exercises,  wargames,  etc.  Therefore,  the  second 
ph  ase  should  sec  the  model  applied  to  some  new  acquisitions.  When  management  is 
comfortable  with  the  model  and  is  ready  to  introduce  it  on  a  broader  scale,  a  master 
implementation  plan  can  then  be  formulated  to  formalise  the  use  of  the  model. 


V.  PROTOTYPE  DESIGN 


A.  GENERAL 

To  foster  a  better  understanding  of  the  DSRPM,  a  prototype  was  developed 
based  on  a  hypothetical  weapon  system.  This  chapter  documents  the  prototype 
which  was  developed  on  an  IBM  AT-compatible  personal  computer  using  PC  SIM- 
SCRIPT  II. 5.  The  language  was  chosen  for  its  superlative  power  and  ease  of  use, 
as  well  as  its  provisions  for  incorporation  of  simulation  in  later  development.  The 
program  is  designed  to  produce  for  each  run  a  single  solution  path  corresponding  to 
a  point  on  Figure  3.2.  The  program  reads  data  from  an  input  file  “IN. DAT”  and 
channels  results  to  an  output  file  “OUT. DAT”.  The  type  of  solution  depends  on  the 
approach  selected  for  that  run,  namely  Lagrangian  relaxation,  maximum  effective¬ 
ness,  minimum  satisficing  cost  or  maximum  system  cost  effectiveness  ratio.  Solution 
approach  is  specified  by  the  user  in  the  input  file.  The  budget  is  relevant  only  when 
the  Lagrangian  method  is  used.  The  overall  structure  of  the  prototype  is  described 
in  Figure  5.1  and  the  program  is  listed  in  Appendix  A. 
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PREAMBLE 


Declaration  of 
entities  and  global 
variables 


->  denotes  links 
conditional  on 
input  information 


CALCl.NODE 


Structuring  the 
system  renewal 
network 


BEST. PATH 

\ 

\ 

MAIN 

Solution  of 
optimal  path 
using  forward 
dynamic 
programing  for 
given  approach 

** - 

\ 

Initialize  system 
renewal  network 
&  solve  for  optimal 
path  given 
approach,  scenario 

\  / 

INITIALIZE 


Read  approach, 
scenario  &  other 
inputs.  Create 
system  renewal 
network  & 
compute  life-cycle 
effectiveness  & 
cost  of  each  arc 


CALC2.NODE 


Computing 
year  &  system 
associated  with 
each  node  & 
arcs  entering  it 


(’AI.C.AI.T.COI)E 


Assigning 
alternative 
code  to  each  arc 


SUM.LCC.LCE 


Compute 
cumulative  life 
cycle  eff.  &  cost  of 
3  systems  in 
optimal  path 


LAGRANGIAN 


Solves  for  optimal 
path  &  associated 
A*  using 

supporting  plane 
method 


CALC.  ARC 


Computation  of 
lifecycle 
effectiveness  & 
cost  associated 
with  each  arc 


ECHO.  IN  PUT 


Check-pi  int  of 
input  and  key 
computational 
values  when 
required 


Figure  5.1  :  Prototype  DSRPM  Block  Diagram 


B.  ASSUMPTIONS 


The  following  assumptions  are  peculiar  to  the  hypothetical  weapon  system 
upon  which  the  prototype  was  developed: 

•  The  effectiveness  sub-model  is  given  by  Equation  4.3  which  was  discussed  as  an 
example  in  the  previous  chapter.  Deterioration  functions  da  j(t)  are  assumed 
to  be  exponential  for  all  a  and  j.  Furthermore,  utility  functions  Uj  {A}  are 
assumed  to  be  linear  and  bounded  (maximum  value  of  1)  as  follows: 


(/;{*) 


X 

max  {A} 


(5.1) 


Therefore,  the  total  effectiveness  of  an  alternative  a  acquired  in  year  y  in  the 
Ith  year  of  use  is  given  by: 


TE.JI)  = 

;=i 

Tl  IVV 

_  J"  CP,. ■».»(<) 

jix  m ax  (X,.,.,) 


(5.2) 


The  cost  sub-model  used  is  also  similar  to  the  example  discussed  in  the  previ¬ 
ous  chapter.  Annual  maintenance  cost  AMCaj(t)  is  assumed  to  be: 


AMCa,y(t)  =  (5.3) 

At  the  end  of  a  system's  life,  all  parts  that  can  be  salvaged  are  sold  and  those 
parts  that  could  not  be  sold  are  disposed  of  at  a  certain  cost.  Therefore,  the 
model  has  both  disposal  cost  and  salvage  value.  Operating  cost,  disposal  cost 
and  salvage  value  are  all  assumed  to  vary  exponentially  through  the  years. 


Therefore,  the  total  life-cycle  cost  incurred  by  an  alternative  a  acquired  in 

year  y  and  used  for  t  years  is: 

LCCa,y{t)  =  ACa,y  +  DC^e"'’1  + 

+ E  +  oc.y-’‘}  (5.4) 

(=i 

C.  MODULES 

This  section  briefly  explains  the  function  and  workings  of  each  component 
module  of  the  prototype. 

1.  PREAMBLE 

This  module  is  peculiar  to  the  S1MSCRIPT  11.5  language.  It  provides 
for  declaration  of  entities  and  global  variables.  Entities  are  structured  data  items 
that  represent  elements  of  a  model.  Definition  of  an  element  as  a  permanent  entity 
entails  the  automatic  creation  of  separate  arrays  for  each  declared  attribute  when  the 
entity  is  created,  the  length  of  these  arrays  being  the  number  of  that  element.  For 
example,  alternatives  were  defined  to  be  the  permanent  entity  ALT  with  attributes 
ALT. NAME  and  LEAD-TIME.  When  the  entity  ALT  is  created,  two  arrays  with 
the  names  of  the  attributes  will  be  created  simultaneously  and  their  lengt  h  would  be 
equal  to  the  number  of  alternatives  N.ALT  (a  system  variable).  Temporary  entities 
would  be  relevant  only  in  the  modeling  of  dynamic  systems  where  elements  enter 
and  leave. 

2.  MAIN 

As  the  name  suggests,  this  module  controls  the  overall  computational 
flow,  which  is  shown  graphically  in  Figure  5.2.  Note  that  the  modules  BEST. PATH 

f.l 


and  SUM.LCC.LCE  are  executed  several  times  if  the  Lagrangian  relaxation  ap¬ 
proach  is  selected. 

3.  INITIALIZE 

This  module  reads  the  input  data,  controls  the  modules  generating  the 
system  renewal  network  and  creates  the  permanent  entities  and  associated  arrays.  It 
also  reserves  the  memory  space  needed  by  the  arrays.  The  flowchart  for  this  module 
is  shown  in  Figure  5.3.  The  input  data  is  read  in  free  format  in  the  order  showm  in 
the  listing. 

4.  CALCl.NODE 

The  structure  of  the  network  is  determined  in  this  module  based  on  the 
number  of  alternatives  and  the  minimum  and  maximum  life-span  of  each  system. 
The  number  of  permanent  entities  NODE  and  ARC  to  be  created  is  also  computed. 

5.  CALC2.NODE 

This  module  computes  for  each  node  information  such  as  the  year  and 
system  (transition)  it  belongs  to,  as  well  as  the  arcs  entering  it.  The  hierarchical  list 
data  structure  illustrated  in  Figure  5.4  was  used  to  represent  the  network.  Therefore, 
the  indices  of  the  arcs  entering  node  n  runs  from  EP(n)  to  [EP(n  -F  1)  —  1]. 

6.  CALC.  ALT.  CODE 

The  code  of  the  alternative  represented  by  each  arc  is  assigned  in  this 

module. 

7.  CALC. ARC 

This  module  determines  the  tail  node,  life-cycle  cost  and  effectiveness 
associated  with  each  arc.  The  one-time  cost  items  in  Equation  5.1  are  discounted  and 
summed  separately  from  the  cumulative  items  such  as  maintenance  and  operating 
costs.  Acquisition  cost  of  the  present  system  is  ignored  as  sunk  cost.  The  benefit 
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accrued  in  each  year  is  added  to  yield  the  cumulative  or  life-cycle  effectiveness  of 
the  arc  representing  a  given  alternative  and  life-span.  If  the  value  of  any  MOE 
falls  below  that  required  in  any  year  spanned  by  an  arc,  the  life-cycle  cost  of  that 
arc  is  set  to  an  unfavorably  high  value  to  mark  it  as  infeasible  and  to  prevent  its 
inclusion  in  the  selection  path.  The  highest  index  of  the  MOE(s)  that  failed  to  meet 
its  requirement(s)  is  recorded  as  the  limiting  MOE  in  such  cases.  The  flowchart  for 
this  module  is  shown  in  Figure  5.5. 


Figure  5-2:  Flowchart  or  MAIN  Program 
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TUANS.  YEAR: 

SYSCODE 

EI’: 

TAM 

ALT. CODE: 


Transition  year 
Preceding  system  code 
Entry  point 
Tail  node 
Alternative 
represented 


Note:  The  convention  shown  here  is  that 
of  reverse  star. 


NODE 


Trans 

Year 

Sys. 

Code 

1 

1989 

0 

2 

1989 

A 

3 

1990 

A 

4 

1991 

A 

5 

1990 

B 

Hierarchical 
List  Data 
Structure 


EP 


0 


1 


2 


3 


4 


Network  Represented 


6 


I 


^  Return 


Includes  checks  of 
MOE  value  (which 
deteriorates)  in  each 
year  to  the  level 
required  in  that  year 
_  Failure  to  meet 
requirement  for  any 
MOF  in  any  year 
spanned  by  an  arc  results 
in  the  life  cycle  cost 
being  set  to  a  very  high 
value 


Figure  5.5:  Flowchart  for  CALC. ARC  Module 


67 


8.  ECHO. INPUT 


Key  inputs  may  be  echoed  for  debugging  purposes  by  setting  the  flag 
variable  ECHO. ON  in  the  input  file  to  “on”.  This  module  will  then  be  activated 
and  intermediate  results  such  as  the  life-cycle  cost  and  effectiveness  associated  with 
each  arc  are  printed  along  with  the  key  input  data. 

9.  LAGRANGIAN 

This  module  produces  the  optimal  path  for  a  given  budget  using  the 
Lagrangian  relaxation  method.  The  value  of  A*  corresponding  to  the  optimal  path 
is  found  by  iterations  of  the  supporting  plane  method.  Each  of  these  iterations 
entails  calling  the  modules  BEST. PATH  and  SUM.LCC.LCE  to  solve  for  L(\i)  and 
the  slope  of  L( A)  at  A,.  The  limitations  of  binary  representation  will  be  exceeded, 
even  with  double  precision  variables,  if  A,  and/or  the  cost  values  are  too  high.  The 
initial  upper  bound  Amar  given  by  Equation  3.18  works  well  for  cost  values  scaled 
to  no  more  than  106.  Figure  5.6  shows  the  flow  of  the  computations  in  this  module. 
Note  that  even  though  the  slope  at  the  minimum  of  L( A)  itself  may  be  positive  or 
negative  depending  on  whether  A3  =  A2  or  A3  =  Aj  respectively,  the  optimal  solution 
path  is  given  by  the  positively  sloped  line  segment  that  is  adjacent  to  the  minimum 
point.  When  A-  =  Aj,  this  line  segment  corresponds  to  A2  (otherwise,  the  intercept 
of  the  line  segments  at  Aj  and  A2  would  not  have  yielded  A3  =  Aj).  Therefore,  the 
slope  is  readjusted  when  A3  =  Aj. 

10.  BEST. PATH 

1  his  module  produces  the  optimal  solution  path  for  the  specified  approach 
using  forward  dynamic  programming.  The  criterion  function  L(J)  associated  with 
each  arc  J  when  the  approach  is  to  maximize  the  effectiveness  to  cost  ratio  of  the 
individual  systems  is  given  by: 


(5.5) 


where: 


L(J)  = 


LCE{J) 

(LCC(J)  -  MIN.LCC  +  1) 


•  LCC(J)  and  LCE(J)  are  the  life-cycle  cost  and  effectiveness  respectively  of 
arc  J 

•  MIN.LCC  =  min  LCC(J) 

This  is  to  allow  for  the  possibility  of  LCC(J)  being  zero  or  negative. 
Provisions  are  also  made  to  ensure  that  the  solution  path  has  sufficient  lead  time 
for  each  alternative  chosen.  The  computational  flow  is  illustrated  in  Figure  5.7. 


(if) 


Figure  5.7;  F  lowchart  for  BEST. PATH  Module 


11.  SUM.LCC.LCE 


The  solution  path  is  traced  backward  from  the  terminal  node  will)  the 
highest  value  of  the  criterior  function  CRIT.FUNC(I)  and  the  life-cycle  cost  and 
effectiveness  of  the  systems  are  summed  along  the  way  in  this  module.  Note  that 
all  the  terminal  nodes  will  have  very  low  values  (negative)  of  CRIT.FUNC(I)  if  all 
pu.ths  violate  either  lead  time  or  MOE  requirements  or  both,  i.e.,  no  feasible  solution 
exists.  Conversely,  if  CRIT.FUNC(I)  of  a  terminal  node  is  not  a  very  low  value, 
below  10~6,  then  a  feasible  solution  path  exists  that  leads  to  that  node. 

12.  RESULTS. PRINT 

Computational  results  are  printed  by  th.c  module. 

D.  COMPLEXITY  AND  EXECUTION  TIME 

In  terms  of  time  complexity,  the  dominating  piece  of  code  is  the  CALC. ARC 
module.  Let  the  difference  between  the  maximum  and  minimum  life-span  of  the 
present  system  be  nt  and  that  of  successor  systems  be  n2,n3,  •  •  • ,  n4,  where  s  is  the 
number  of  stages  or  systems  in  the  systems  renewal  network.  If  we  assume  n  —  ??, 
for  V?  >  1,  then  the  structure  of  the  network  is  given  by  Table  5.1. 

I  he  total  number  of  arcs  n  is  thus  given  by: 

j  —  2 

r?  =  1  +  7?!  -f  o(s  -  1)(1  +  n,)(l  +  n)  +  an(l  -f  n)  ^  i 

i=i 

=  1  -f  r?!  -f  a(s  —  1)(1  -f  r? j ) ( 1  -f  n)  -f  an(l  4-  n)(s  —  2)(s  —  l)/2  (5.6) 

If  we  let  the  average  life  of  successor  systems  be  denoted  as  /,  and  the  number  of 
MOEs  evaluated  be  c,  then  the  growth  rate  of  execution  time,  ET,  is  of  the  order 
irt(l  +  r).  In  other  words. 
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TABLE  5.1:  System  Renewal  Network  Structure  (n j  =  n  for  V,  >  1) 


Transition  # 

System  #i 

No.  of  Nodes 

No.  of  Arcs 

0 

1 

1 

(1  +  ni) 

1 

(1  +  ni) 

2 

a(l  +  ni)(l  +  n) 

2 

(1  +  ni  -P  n) 

3 

a(l  -f  n\  +  n)(  1  +  n) 

3 

(1  +  til  +  2n ) 

4 

a(l  +  U\  -p  2n)(l  -p  n) 

4 

(1  -P  i  T  3n) 

And 

so  on  until  i  = 

s 

Note:  a  denotes  the  number  of  alternatives  evaluated 


ET  —  0[m(/  +  e)] 


=  0 


1  +  Tij  +  a(s  -  1)(1  +  n,)(l  +  n)  -f 


an2(s  —  2)(s  —  1) 


(*  +  e) 


=  0  m&x(an,  ns,  an2s2)ma\(l,  e)j 


(5.7) 


Therefore,  execution  time  for  the  prototype  DSRPM  is  of  polynomial  order,  the  key 
determinants  of  which  arc  n  and  ,s.  We  can  also  see  that  time  complexity  is  linear 
in  o,n,f  and  c.  For  a  run  with  a  =  2,ii\  =  5, ri  =  5,5  =  3,1=  7.5  and  e  =  3,  the 
execution  time  of  the  prototype  DSRPM  on  the  Zenith  AT-compatible  is  typically 
about  43  seconds. 


VI.  CASE  STUDIES 


Case  studies  were  conducted  using  the  prototype  to  demonstrate  the  plausi¬ 
bility  and  usefulness  of  the  DSRPM.  To  better  serve  the  purposes  of  illustration  and 
plausibility  testing,  the  case  studies  were  designed  to  comprise  a  series  of  scenarios 
of  incremental  complexity,  all  of  which  are  based  on  the  same  hypothetical  system 
discussed  in  Chapter  V.  Although  effort  was  made  to  incorporate  as  much  realism  as 
possible,  the  case  studies  and  the  weapon  system  were  kept  simple  (o  better  portray 
the  workings  of  the  model.  A  good  understanding  of  and  confidence  in  the  model  is 
necessary  before  one  confronts  the  complexities  involved  in  a  real  application. 

In  the  case  studies,  the  key  effectiveness  parameters  or  MOEs  arc  assumed  to 
be  range  (nautical  miles),  accuracy  (percentage  hit)  and  availability  (percentage) 
with  relative  weight  W)  of  0.4,  0.3  and  0.3  respectively.  The  cost  and  effectiveness 
characteristics  of  the  incumbent  system  assumed  arc  shown  in  Table  0.1.  Disposal 
cost  is  assumed  to  be  zero  for  the  incumbent  system  and  alternatives,  i.e.,  DC„,y  =  0 
for  Va . 

TABLE  6.1:  Characteristics  of  the  Incumbent  System 


Current  Value 

Po 

Cost  Parameters  ( 

j  Range  (NM) 

104.6 

-0.03 

Maintenance 

Operating 

Salvage 

ujeesesmesh 

87.2 

-0.02 

h 

70 

OC0 

EH 

S\b 

^0 

93.5 

-0.03 

30 

_ 

1.3 

7 

0.05 

-20 

-0.2 
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As  a  starting  point,  we  shall  assume  that  two  alternatives,  1  and  2,  are  available 
and  that  their  minimum  and  maximum  life  are  five  years  and  ten  years  respectively. 
They  are  assumed  to  require  no  development  or  acquisition  lead  time  (presumably 
available  “off-the-shelf”).  The  incumbent  system  (alternative  0)  is  assumed  to  have  a 
minimum  and  maximum  life  of  zero  years  and  five  years  respectively,  i.e.,  alternatives 
may  be  introduced  immediately,  or  at  the  latest,  five  years  hence. 

Four  solution  approaches  are  used  in  each  of  the  case  studies,  namely  the 
minimum  average  annual  cost  (MAAC)  approach,  the  minimum  annualised  cost 
(MAC)  approach,  the  system  effectiveness  to  cost  ratio  (SECR)  approach  and  the 
Lagrangian  relaxation  (LR)  approach.  The  average  annual  cost  is  obtained  by  di¬ 
viding  the  total  cost  of  a  succession  of  three  systems  (in  current  dollars)  by  the  years 
spanned  by  the  systems.  The  MAAC  approach  thus  ignores  monetary  interest.  The 
MAC  solution  is  obtained  by  amortizing  the  total  cost  of  the  system  over  the  years 
spanned.  It  is  the  solution  obtained  when  one  wishes  to  minimise  cost  needed  to 
meet  effectiveness  requirements.  The  SECR  approach  maximizes  the  cumulative 
total  of  the  elTcctivcncss-to-cost  ratio  of  each  of  the  three  systems  in  a  path.  The 
LR  approach  was  explained  in  Chapter  3.  Solution  points  obtained  using  the  LIt 
approach  are  linked  as  in  Figure  3.5  to  form  the  efficiency  boundary.  The  point  on 
the  boundary  that  yields  the  maximum  overall  effectiveness  to  cost  ratio  would  be 
the  preferred  solution  if  one  adopts  the  approach  of  maximising  returns  on  cost. 

A  notational  point  to  note  is  that  the  transition  years  shown  in  the  results  of 
the  case  studies  are  all  abbreviated  to  the  last  two  digits.  That  is,  97  and  05  mean 
1997  and  2005  respectively.  The  current  year  is  1989. 
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A.  CASE  STUDY  1:  STATIONARY  ENVIRONMENT 


The  first  scenario  we  examine  is  one  in  which  the  only  change  occurring  is 
that  due  to  deterioration.  Technology  and  requirement  levels  are  at  a  standstill  and 
no  inflation  occurs.  The  value  of  money  is  constant  over  time,  i.e.,  interest  rate  is 
zero.  This  highly  improbable  scenario  is  implied  by  the  commonly  used  “equal  life” 
assumption  whereby  an  alternative,  once  selected,  will  be  selected  over  and  over 
again  and  will  always  have  the  same  service  life. 

Supposing  alternatives  1  and  2  have  the  following  cost  and  effectiveness  char¬ 
acteristics  shown  in  Table  6.2  and  6.3.  Alternative  2  is  thus  three  times  as  expensive 
to  acquire  but  has  a  longer  range  and  otters  lower  maintenance  and  operating  costs 
than  Alternative  1.  Requirement  levels  are  assumed  constant  at  95  nni  for  range, 
80%  for  accuracy  and  70%  for  availability.  Results  of  analysis  performed  using  the 
prototype  DSRPM  for  this  case  is  presented  in  Table  6.4  below. 

TABLE  6.2:  Characteristics  of  Alternative  1 


MOE 

|  Value  When  New 

P\ 

Cost  Parameters 

!  Rangel  nm) 

i  120 

-0.03 

i  Acquisition 

Maintenance 

Operating 

Salvage 

Accuracy  (7  i 

1  95 

-0.02 

AC, 

li 

OC,  I  7, 

sr,  j  0, 

Availability  (7  ) 

i  95 

-0.03  1;  100 

30 

1.3 

7  1  0.05 

-50  .  -0.2 

I  he  value  of  the  total  life-cycle  cost  and  effectiveness  are  plotted  in  Figure  61. 
As  interest  rate  is  zero,  MAAC  and  MAC.'  approaches  yielded  the  same  solution.  The 
results  show  that  an  alternative  that  is  selected  will  be  selected  repeatedly  in  the 
future,  the  service  life  each  time  being  the  same.  This  is  only  to  be  expected  as, 
under  st  at  ic  con  d  it  ions,  t  here  ai  e  no  reasons  for  change.  Note  that  the  solution  that 
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TABLE  6.3:  Characteristics  of  Alternative  2 


j  MOE 

Value  When  New 

Pi 

)  Cost  Parameters 

!  Range  (nm) 

140 

-0.03 

1  Acquisition 

Maintenance 

Operating 

Salvage  | 

Accuracy  (7c) 

95 

1  -4C2 

m 

72 

OCi 

K9 

S\\ 

Availability  (V7) 

95 

-0.03 

15 

13 

4 

EES 

TABLE  6.4:  Results  for  Stationary  Environment 


1  Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt. 

Chosen 

Transit 

on  Year 

1 

I  Approach 

(3  Systems) 

(3  Systems) 

Ratio 

a 

B 

3 

D 

B 

3 

|  MAAC 

120G.05 

1G.550 

0.01372 

□ 

B 

2 

m 

07 

MAC 

120G.05 

□ 

B 

E3 

E 

07 

1  SECR 

1097.35 

0.01373 

a 

B 

2 

89 

Ea 

05 

LR 

1 270.64 

17.2G5 

0.01359 

□ 

B 

2 

EEl 

Eil 

OS 

0 

LR 

120G.05 

1G.550 

0.01372 

n 

R 

MB 

m 

0.01104 

LR 

1147.33 

lo.SIG 

•0.01379 

□ 

B 

2 

98 

0.01251 

LR 

801.84 

10. 934” 

0.013G4 

0 

i 

1 

90 

95 

00 

0.01413 

LR 

751. SG 

10.179 

0.01354 

a 

D 

89 

94 

t.J9 

0.01509 

gives  the  maximum  return  per  unit  cost  for  the  three  systems  overall  (marked  with 
a.'  aster  >k  in  Table  0.4  and  a  square  in  Figure  6.1)  is  different  from  that  obtained 
using  t.ie  SECR  approach.  The  overall  LCE  to  OCC  ration  in  the  case  of  the  SECR 
solution  is  not  as  high  even  though  the  solution  also  lies  on  the  efficiency  boundary. 
The  implication  of  this  is  obvious  -  maximizing  return  per  unit  cost  for  individual 
systems  considered  separately  leads  to  a  good  but  not  best  solution  in  terms  of 
overall  return  per  unit  cost. 
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Figure  6.1:  Stationary  Environment  Case  Study  Results 

B.  CASE  STUDY  2:  INFLATION  AND  INTEREST  RATES 

The  next  scenario  we  look  at  is  probably  the  one  most  commonly  encountered 
in  a  typical  system  renewal  stud}'.  Future  inflation  and  interest  rates  are  assumed 
to  be  non-zero  and  constant,  and  future  technological  advancements  and  changes 
in  requirement  levels  are  ignored.  For  illustrative  purposes,  let  us  say  the  inflation 
and  interest  rates  were  projected  to  be  constant  at  four  and  five  percent  per  annum 
respectively.  Cost  projections  based  on  these  rates  are  shown  in  Figure  6.2  in  the 
sample  input  file  for  this  case  study  (Appendix  B).  Analysis  based  on  this  scenario 
produced  the  results  shown  in  Table  6.5. 


Coat  Cost 

200  400  600  800  50  100  150  200  250 


Acquisition  Year 

Figure  6  2b.  Cost  Evolution  toi  New  Alternative  2 

Figure  6.2  :  Evolution  of  Cost  Parameters  for  Alternatives  1  and  2 
With  Inflation  {  No  Technological  Progress  ) 


TABLE  6.5:  Results  for  Inflation  and  Interest  Rates 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt.  Chosen 

Transition  Year 

A’ 

Approach 

3  Systems 

Ratio 

D 

□ 

3 

1 

K9 

3 

MAAC 

1031. 2G 

17.266 

0.01674 

0 

2 

92 

liTl 

mm 

MAC 

811.98 

13.375 

0.01647 

0 

2 

90 

98 

SECR 

921.03 

14.396 

0.01563 

2 

91 

97 

04 

LR 

1031.2G 

17.265 

0.01674 

2 

92 

00 

0 

LR 

9SG.44 

16.550 

★0.01678 

2 

91 

99 

0.01596 

LR 

944.29 

15.816 

0.01675 

2 

90 

98 

0.01742 

LR 

811.9S 

13.375 

0.01647 

0 

2 

1 

90 

98 

03 

0.01845 

LR 

683.34 

10.934 

0.01600 

0 

1 

1 

90 

95 

00 

0.01898 

LR 

G43.73 

10.179 

0.01581 

0 

1 

1 

89 

94 

99 

0.01904 

Total  LCC  and  LCE  from  this  table  are  plotted  in  Figure  6.3.  The  efficiency 
boundary  is  shifted  towards  lower  cost  since  the  interest  rate  is  higher  than  the 
inflation  rate.  An  opposite  movement  would  occur  if  the  converse  was  true.  Note 
that  the  shift  is  greater  for  the  higher  cost  end  because  inflation  and  interest  rates 
have  a  bigger  impact  there.  As  the  interest  rate  is  significantly  non-zero,  we  see 
that  the  solutions  obtained  using  the  MAAC  and  MAC  approaches  are  different. 
Note  also  that  the  SECR  solution  is  clearly  inferior  as  it  lies  somewhat  below  the 
efficiency  boundary,  from  the  above  table,  we  can  sec  that  Alternative  1  is  the 
low-end  alternative,  f  he  last  line  in  the  table  shows  the  solution  that  meets  the 
specified  minimum  requirements  at.  the  lowest  possible  total  cost,  which  is  the  one 
where  the  service  lives  of  all  three  systems  arc  at  their  minimums. 
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m 


-  Efficiency  Boundary 

(formed  by  Lagrongian  Solution  Points) 


Figure  6.3:  Inflation  and  Interest  Rate  Case  Study  Results 

C.  CASE  STUDY  3:  TECHNOLOGICAL  ADVANCEMENT 

For  this  scenario,  let  us  suppose  that  technological  progress  has  the  efTect  of 
suppressing  cost  growth  due  to  inflation  from  four  to  two  percent  per  annum.  The 
revised  cost  projections  are  shown  in  Figure  6.4.  The  sample  input  file  for  this  case 
study  is  shown  in  Appendix  D.  In  addition,  technological  advances  are  also  reflected 
in  improvement  in  performance  characteristics.  For  the  purpose  of  this  case  study, 
let  us  suppose  that  the  various  MOEs  evolve  as  shown  in  Figure  6.5  and  the  sample 
input  hie  in  Appendix  B.  Improvements  are  primarily  in  range,  with  Alternative  1 
becoming  marginally  better  than  Alternative  2  about  ten  years  from  now.  Other 
conditions  arc  as  assumed  in  Case  Study  2.  The  results  for  this  case  study'  are  shown 
in  Figure  6.6  and  Table  6.6. 
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Contrasting  Table  6.5  and  Table  6.6,  we  note  that  Alternative  1  is  now  the 
dominantly  favored  alternative  because  it  became  competitive  performance-wise  as 
well  as  cost-wise.  Technological  advancement  resulted  not  only  in  a  further  shift 
in  the  efficiency  boundary  towards  lower  cost  but  also  an  extension  into  higher 
levels  of  cumulative  effectiveness.  The  lower  end  of  the  boundary  is  not  affected 
as  much  as  solutions  in  the  higher  end  that  extends  further  into  the  future.  Note 
that  technological  advancement  also  results  in  lengthening  of  the  maximum  service 
life  (three  systems),  as  evident  in  the  maximum  effectiveness  solution  identified  by 
A*  =  0.  In  Case  Studies  I  and  2,  no  solutions  extend  beyond  the  year  2 008  whereas 
in  this  case  study,  some  solutions  extend  to  the  year  2009. 

TABLE  6.6:  Results  for  Technological  Advancement 


j  Solution 

Total  LCC 

Total  LCE 

LCE /LCC 

Alt.  Chosen 

Transition  Teat 

— -•  r, 

A!  1 

i  Approach 

3  Systems 

3  Systems 

Ratio 

1 

2 

3 

1 

2 

3 

■ 

|  MAAC 

SS6.1I 

18.491 

0.020S2 

0 

nr 

2 

92 

00 

09 

• 

|  MAC 

736.01 

14.509 

0.01971 

0 

i 

1 

91 

99 

04 

j  SECR 

804.94 

16.220 

0.02015 

0 

i 

1 

92 

99 

06 

i  LR 

907.83 

18.77o 

0.02068 

0 

2 

2 

92 

00 

09 

0 

i  LR 

888.11 

18.491 

*0.02082 

0 

1 

1 

92 

00 

09 

.  0.01444 

LR 

634.53 

12.219 

0.01926 

0 

2 

1 

92 

96 

01 

0.024"  : 

j  LR 

601.75 

11.366 

0.018S9 

0 

1 

1 

90 

95 

00 

0.0260! 

!  LR 

574.70 

10.506 

0.01828 

0 

1 

1 

89 

94 

99 

'  (1.03180 
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— • —  Volue  of  1st  MOE  for  Alternative  1 

. * Volue  of  1st  MOE  for  Alternotlve  2 

--H —  Volue  of  2nd  end  3rd  MOES  for  both  Alt. 


a 

0 


Note  :  Technological  progress  Is  primarily  In 
the  let  MOE  for  both  alternatives.  Improvement 
in  the  other  MOEe  ore  morginol. 


o  1  J  I _ I _ I - 1 - 1 - 1 - 1 - 1 - 1 - 

1000  1085  2000  2005  2010 

Acquisition  Year 

Figure  6.5:  Evolution  of  MOEs  With  Technological  Progress 
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Figure  6.G:  Technological  Advancement  Case  Study  Results 
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D.  CASE  STUDY  4:  REQUIREMENT  ESCALATION 


Let  us  suppose  that  everything  is  as  before,  except  that  requirement  levels  set 
for  range  and  availability  increase  over  time  as  shown  in  Figure  6.7.  The  results  for 
this  scenario  are  shown  in  Table  6.7  and  illustrated  in  Figure  6.8. 


TABLE  6.7:  Results  for  Requirement  Escalation 


j  Solution 

Total  LCC 

Total  LCEl 

LCE/LCC 

Alt.  Chosen 

Transition  Year 

A* 

i  Approach 

3  Systems 

3  Systems 

D 

B 

3 

1 

H 

3 

1  MAAC 

741.G2 

14.498 

0.01955 

0 

i 

1 

91 

9G 

04 

1  MAC 

741.G2 

14.498 

0.01955 

0 

i 

1 

91 

9G 

04 

!  SECR 

771.83 

14.8G5 

0.0192G 

0 

2 

1 

91 

98 

04 

!  LR 

872.43 

1G.431 

0.01883 

0 

2 

2 

91 

98 

0G 

0 

j  LR 

835.93 

1G.395 

*0.01961 

0 

2 

1 

91 

98 

0G 

0.00097 

LR 

741.G2 

14.498 

0.01955 

0 

1 

1 

91 

9G 

04 

0.02011 

LR 

G34.53 

12.219 

0.0192G 

0 

1 

1 

91 

9G 

01 

0.02129 

LR 

574.70 

10.506 

0.01828 

0 

1 

1 

89 

94 

99 

0.028G3 

Comparing  Figures  6.6  and  6.8,  we  can  see  that  the  escalation  in  requirement 
levels  resulted  in  a  downward  depression  of  the  higher  cost  end  of  the  efficiency 
boundary.  This  happens  because  the  escalation  in  requirement  levels  results  in  the 
shortening  of  service  lives  only  if  the  solution  path  extends  far  enough  to  be  affected. 
Note  that  the  maximum  effectiveness  solution,  for  instance,  now  extends  only  until 
the  year  2006  whereas  in  the  previous  case  study  it  reaches  the  year  2009.  An 
interesting  result  in  the  equality  of  the  M  AAC  and  MAC  solutions  even  though 
the  interest  rate  is  five  percent  per  annum.  Note  also  that  solutions  calling  for  the 
first  transition  to  be  made  in  1991  arc  clearly  dominating. 


Total  Eff.  (  3  Sys.  Sequence  )  MOE  Req.  Lerel  (  n.m.  or  percent  ) 
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Figure  6.7:  Escalation  in  Requirement  Levels 


Figure  6.8:  Requirement  Escalation  Case  Study  Results 


SO 


E.  CASE  STUDY  5:  SCHEDULING  OF  UPGRADES 


Supposing  Alternative  1  is  an  upgrading  option,  the  exercise  of  which  is  sig¬ 
nificantly  cheaper  if  carried  out  along  with  major  scheduled  overhauls.  If  such 
maintenance  is  carried  out  every  three  years,  (e.g.,  refit  of  a  ship),  and  the  savings 
is  say  twenty  percent,  then  the  acquisition  cost  of  a  new  Alternative  1  varies  as 
illustrated  in  Figure  6.9.  In  this  scenario,  everything  else  is  as  assumed  for  Case 
Study  4. 

The  results  for  this  scenario  are  given  in  Table  6.8  and  Figure  6.10.  As  evident 
from  the  table,  there  are  no  saving  realised  in  the  higher  cost  end  where  the  discounts 
are  not  taken  up.  The  lower  end  of  the  efficiency  boundary  shifts  to  the  left  as  a 
result  of  the  savings.  We  see  that  the  solution  offering  maximum  effectiveness  per 
unit  cost  suggests  that  the  incumbent  should  be  replaced  with  Alternative  2  in 
1990,  then  upgraded  with  Alternative  1  in  1997  (a  non-maintenance  year  with  no 
discount)  for  use  until  the  year  2005.  The  SECR  solution  also  indicates  an  upgrade 
carried  out  in  a  non-maintenace  year  (1998).  The  overall  desirability  of  matching  a 
major  upgrading  program  to  some  maintenance  schedule  to  minimise  downtime  or 
to  secure  some  savings  can  thus  be  evaluated  using  the  model. 
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Total  Eff.  (  3  Sys.  Sequence  )  £  Cost 


Not*  :  Ths  acquisition  cost  of  ths  - Acquisition  Cost 

upgrading  Is  lawsr  only  In  thoss  .  1st  Year  Matntsnoncs  Cost 


_  l . _j . J  _ .  _.J _ 1  _  J _ I _ l  — -l  -  1  < 

1690  1695  2000  2D05  2010  2015 


Year  of  Upgrading 

:  6.9:  Evolution  of  Cost  Parameters  With  Scheduled  Overhauls 


Figure  6.10:  Scheduling  of  Upgrading  Case  Study  Results 
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TABLE  6.8:  Results  for  Upgrade  Scheduling 


|  Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt.  Chosen 

Transition  Year 

a: 

;  Approach 

3  Systems 

3  Systems 

Ratio 

D 

a 

3 

1 

2 

3 

|  MAAC 

741.62 

14.498 

0.01955 

0 

1 

1 

91 

96 

04 

MAC 

653.39 

12.734 

0.01949 

0 

1 

1 

89 

94 

02 

SECR 

771.83 

14.865 

0.01926 

0 

2 

1 

91 

9S 

04 

LR 

872.43 

16.431 

0.01883 

0 

2 

2 

91 

98 

06 

0 

LR 

835.93 

16.395 

0.01961 

0 

2 

91 

98 

06 

0.00097 

LR 

791.73 

15.532 

*0.01962 

0 

2 

90 

97 

05 

0.01953 
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Note:  In  this  case  stud}-,  a  1  for  alternative  chosen  represents  an  upgrade  of  the 
preceding  system.  Thus  the  sequence  0,2,1  means  the  incumbent  is  replaced  with 
alternative  2  and  then  upgraded  with  alternative  1. 

F.  CASE  STUDY  6:  DEVELOPMENT  LEAD  TIME 

Let  us  suppose  that  everything  is  as  in  Case  Study  5  except  that  Alternative 
1  is  an  upgrade  that  requires  one  year  of  development  and  cannot  be  introduced 
immediately.  The  impact  of  such  lead-time  requirements  can  be  seen  in  Table  6.9 
and  Figure  6.11.  From  the  table,  we  can  sec  that  the  lowest  cost  solution  (bottom  of 
table)  still  involves  an  immediate  transition  as  in  the  other  case  studies.  However, 
the  transition  is  to  Alternative  2  in  this  case  as  the  model  now  bars  the  immediate 
execution  of  Alternative  1.  T  he  diagram  clearly  shows  that  the  lead-time  require¬ 
ment  pushes  up  cost  at  the  lower  cost  end  of  the  efficiency  boundary.  The  MAAC 
solution,  the  SECR  solution  and  the  maximum  effectiveness  per  unit  cost  solution 
are  all  unafTccled.  The  MAC  solution,  however,  became  equal  to  the  MAAC  solu¬ 
tion  (curiously).  Note  that  the  MAC  solution  in  Case  Study  5  involved  immediate 
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introduction  of  Alternative  1  (upgrade  of  alternative  0  preceding  it);  a  solution  that 
is  not  admissible  in  this  case  study. 


TABLE  6.9:  Results  for  Development  Lead-Time 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Al 

it.  Chosen 

Transition  Year 

A* 

Approach 

3  Systems 

3  Systems 

Ratio 

D 

B 

3 

1 

a 

3 

MAAC 

741.G2 

14.498 
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Figure  6.11:  Development  Lead  Time  Case  Study  Results 


G.  GENERAL  OBSERVATIONS 

From  the  series  of  six  case  studies  analysed  using  the  prototype  DSRPM,  we 
saw  the  impact  of  various  issues  as  they  were  integrated  into  the  model.  DSRPM  of¬ 
fers  a  common  framework  for  decision  makers  to  explore  the  consequences  of  various 
possibilities.  It  does  not  advocate  any  single  solution  for  a  given  scenario,  but  rather 
maps  out  the  efficiency  boundary  spanning  the  entire  solution  space  and  indicates 
the  solutions  obtained  using  various  approaches.  A  plot  of  the  results  for  various 
scenarios,  as  is  done  in  Figure  G.12  for  the  case  studies,  gives  the  decision  maker  a 
good  overview  of  the  problem.  The  picture  can  be  enriched  by  varying  the  shade  of 
the  result  for  each  scenario  according  to  its  likelihood  of  occurrence.  Furthermore, 
a  third  axis  indicating  overall  service  time  for  the  three  systems  can  be  added  as 
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in  Figure  6.13  to  give  fuller  expression  to  the  problem.  As  explained  in  Chapter  2, 
the  system  renewal  problem  is  basically  a  three-dimensional  one  -  cost,  effectiveness 
and  time. 
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Figure  6.12:  Consolidated  Results  for  Case  Studies 


Although  three  systems  are  analysed  in  the  model,  the  focus  is  really  on  when 
the  first  transition  should  take  place  and  what  alternative  should  be  selected  to 
succeed  the  incumbent  system.  Analysis  may  reveal  a  dominating  (good  under 
all  scenarios)  or  robust  (favored  under  most  scenarios)  solution  to  the  liming  and 
choice  involved  in  the  first  transition.  For  example,  the  escalation  of  requirement  in 
Case  Study  4  resulted  in  hastening  of  the  first  transition  by  a  year  for  the  maximum 
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Figure  6.13:  3-D  Plot  of  Efficiency  Boundaries 

overall  effectiveness  to  cost  ratio  solution.  In  Case  Study  6,  however,  we  saw  that  this 
solution  was  not  affected  by  the  development  lead  time  of  one  year.  By  varying  the 
scenarios  for  various  requirement  escalation  trends,  development  lead  time  estimates, 
etc.,  we  can  test  the  sensitivity  of  various  solutions  to  variations  in  estimates  and 
trends.  Even  if  analysis  proved  inconclusive,  a  rough  idea  of  how  long  the  incumbent 
is  likely  to  be  kepi  will  emerge.  No  commitment  needs  to  be  made  if  the  time  to 
the  likely  date  of  the  first  transition  is  more  than  the  lead  time  required  to  effect 
the  likely  transition.  For  example,  if  analysis  shows  that  an  E\V  suite  upgrade  of  a 
frigate  should  take  place  three  years  down  the  road  and  such  an  upgrade  has  a  lead 
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time  of  two  years,  we  needn’t  make  a  commitment  yet,  but  we  know  the  decision 
is  due  a  year  from  now.  At  that  time,  the  model  can  be  run  again  with  updated 
information.  In  this  way,  we  can  ensure  decisions  are  made  in  a  timely,  fashion  using 
the  latest  information. 

Besides  facilitating  timely,  consistent  and  rational  decision-making,  the  model 
also  looks  forward  to  the  future  to  provide  estimates  of  cost,  effectiveness  and  service 
life  that  facilitates  good  planning.  If  aggregated  to  sufficiently  high  levels,  the  model 
allows  future  year  procurement  budgets  to  be  forecasted.  Designs  and  support  plans 
for  succeeding  systems  can  be  optimised  to  service  lives  that  are  more  realistic 
in  terms  of  broad  economic,  technological  and  operational  outlook.  Amortization 
periods  would  no  longer  need  to  be  arbitrarily  assumed  and  budgets  can  be  better 
planned  knowing  cost  estimates  of  succeeding  systems.  The  anticipated  system 
selections  also  serves  to  better  guide  planning  for  defense  infrastructure,  R&D  and 
other  strategic  concerns.  The  explicitness  of  the  model  fosters  better  consensus 
building  and  coordination  among  the  various  agencies  responsible  for  system  renewal 
planning.  As  discussed  in  Chapter  4,  Section  E,  force  level  and  mix  issues  as  well 
as  phased  build-up  versus  bulk  procurement  decisions  can  also  be  dealt  with  by 
the  DSRPM.  The  model  thus  provides  a  comprehensive  and  systematic  framework 
for  system  renewal  decisions  to  be  made  optimally  and  in  a  consistent  and  timely 
manner  based  on  current  forecasts. 

Lastly,  it  should  be  noted  that  inflation  and  interest  rates  were  assumed  con¬ 
stant  and  linear  utility  functions  were  used  in  the  case  studies  for  reason  of  conve¬ 
nience  and  not  necessity.  The  model  allows  these  rates  and  effectiveness  requirement 
levels  to  be  varied  for  future  years  in  a  scenario.  Long-range  projections  can  be 
varied  more  than  near-term  ones  to  account  for  greater  uncertainty  in  these  projec¬ 
tions.  The  impact  of  technological  advancement  is  also  freely  interpreted  in  terms 
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of  changes  in  cost  and  effectiveness  parameter  projections  and  no  restrictions  are 
placed  on  the  definition  of  utility  functions.  The  DSRPM  thus  affords  very  flexible 
means  to  model  the  factors  influencing  systems  renewal  problems. 


VII.  CONCLUSIONS  AND 
RECOMMENDATIONS 


A.  CONCLUSIONS 

The  dynamic  system  renewal  planning  model  (DSRPM)  provides  a  system¬ 
atic  framework  within  which  optimal  system  renewal  decisions  may  be  made  in  a 
consistent  and  timely  manner,  under  realistic  conditions  of  changing  environments. 
Conceptually,  it  embodies  a  concept  of  obsolescence  that  is  more  comprehensive, 
better  formalised  and  more  relevant  to  military  planners  than  is  available  in  current 
literature.  Prediction  of  what  future  renewals  are  likely  to  be  made  and  when,  and 
estimation  of  the  associated  cost  and  effectiveness  would  allow  plans  to  be  made 
anticipatively  with  a  long-term  perspective.  Through  simple  reformulation  and  use 
of  shadow  prices,  bulk  versus  phased  procurement  decisions  as  well  as  force  level 
and  mix  issues  can  also  be  examined  using  the  model. 

The  plausibility  and  usefulness  of  the  model  was  demonstrated  by  a  series  of 
case  studies  using  the  prototype  that  was  developed.  The  equal  life  assumption 
often  used  in  system  renewal  models  was  illustrated  by  the  first  case  study.  Other 
case  studies  showed  the  impact  of  interest  and  inflation  rates,  technological  progress, 
requirement  escalation  and  development  lead  time.  Application  of  the  model  to  the 
scheduling  of  upgrades  that  are  cheaper  to  carry  out  during  overhaul  periods  was 
also  shown. 

The  prototype  DSRPM  illustrates  how  the  model  may  be  used  as  a  descriptive 
decision  aid  to  better  understand  the  complex  effects  of  a  diverse  range  of  issues  on 
the  system  renewal  decision.  It  was  shown  how  the  model  organises  financial,  oper¬ 
ational,  intelligence  and  other  staff  projections  within  an  integrated  framework  to 
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provide  decision  makers  with  a  broad,  consistent,  long-term  perspective  of  relevant 
issues.  The  decision  maker  is  presented  with  a  graphic  picture  of  the  entire  solu¬ 
tion  space  as  structured  by  the  system  renewal  scenarios  considered.  Scenarios  are 
flexibly  modelled  and  are  used  to  point  to  solutions  that  are  robust  in  the  face  of  un¬ 
certainty.  By  continually  refining  scenarios  based  on  lessons  learned  from  exercises, 
audits,  and  other  new  information  on  trends,  past  experience  is  brought  to  bear  on 
the  present  so  that  the  basis  of  decisions  is  kept  current  and  realistic.  The  commu¬ 
nication  among  all  parties  involved  in  this  process  is  facilitated  by  the  explicitness 
of  the  model.  This  in  turn  fosters  consensus  building  necessary  for  cohesive  effort 
in  defense  planning. 

The  model  that  was  developed  is  a  conceptual  prototype  that  is  as  yet  untested 
by  a  real  application.  Much  remains  to  be  done  before  it  is  ready  for  organisation¬ 
wide  implementation.  However,  it  holds  much  promise  and  should  be  examined 
further  for  future  implementation. 

13.  RECOMMENDATIONS 

The  prototype  that  was  developed  is  a  rudimentary  one  that  could  be  extended 
in  many  ways  as  future  research  topics.  It  is  recommended  that  research  proceed  in 
the  following  order  of  tractability: 

1.  Factoring  foreign  exchange  rales  into  the  model  is  a  desirable  goal.  The  treat¬ 
ment  of  interest  and  inflation  rates  would,  however,  require  some  adjustments. 

2.  Optimal  force  mix  can  be  determined  if  shadow  prices  of  different  components 
of  a  task  force  are  made  compatible.  Means  of  integrating  the  utility  functions 
of  these  components  under  a  common  mission  structure  needs  to  be  examined. 
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3.  Simulation  and  forecasting  capabilities  can  be  incorporated  to  better  explore 
the  range  of  system  renewal  scenarios  that  may  be  encountered  in  the  future. 

4.  Stochastic  models  can  be  constructed  to  better  portray  uncertainty  in  cost, 
effectiveness  and  requirement  levels. 

5.  Utility  modelling  may  be  facilitated  by  the  use  of  implicit  approaches  that  do 
not  require  utility  functions  to  be  explicitly  defined  over  the  entire  domain  of 
the  MOE  value. 

6.  System  renewal  models  developed  for  various  systems  could  be  linked  and  ag¬ 
gregated  to  higher  entities  using  schemes  that  incorporate  organisation  struc¬ 
ture  into  the  solution  algorithm  to  avoid  sub-optimisation.  Plans  developed 
for  ships  could,  for  example,  be  tied  to  plans  for  their  operational  task  forces 
and  administrative  units.  The  ideal  would  be  to  achieve  a  level  of  aggregation 
high  enough  to  interface  with  arms  race  and  econometric  models  so  that  very 
high  levels  of  optimisation  may  be  achieved. 

7.  Fuzzy  set  and  logic  theory  may  be  introduced  to  better  model  elements  that 
are  inherently  vague.  Requirement  levels,  for  example,  may  be  made  elastic  by 
using  what  is  known  as  fuzzy  numbers.  This  may  allow  some  form  of  abstract 
qualitative  reasoning  to  be  represented  in  the  model. 


APPENDIX  A:  PROGRAM  LISTINGS 


X 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


preamble 

normally,  mode  is  real,  dimension  is  0 
permanent  entities 

every  YEAR  has  a  DISCOUNT. FACTOR  and  a  CUM.DF 
normally,  mode  is  integer 

every  NODE  has  a  TRANS. YEAR,  a  SYS. CODE, 

a  BEST. PRED. NODE  and  a  BEST. PRED. ARC 

every  ARC  has  a  TAIL  and  an  ALT. CODE 

every  SYS  has  a  MIN. LIFE,  a  MAX. LIFE,  a  NUM. NODES,  a  NUM.ARCS  and 
a  BASE. YEAR 

every  ALT  has  an  ALT. NAME  and  a  LEAD-TIME 

define  MOE,  MOE.DF  as  3-dimensional  real  arrays 
define  AC,  MC.P1,  MC.P2,  OC,  OC.EF,  DC,  DC.EF,  SV,  SV.EF  and 
MIN . MOE . REQ  as  2-dimensional  real  arrays 
define  MOE. NAME  as  a  1-dimensional  text  array 
define  L  and  CRIT.FUNC  as  1-dimensional  double  arrays 
define  LCE ,  LCC,  MOE.WT,  O.MOE  and  O. MOE.DF  as  1-dimensional 
real  arrays 

define  MAX . PRED. NODES ,  LIM.MOE  and  EP  as  1-dimensional  integer  arrays 
define  ARGMAX.CF,  ARG.MAC,  ARG.AAC,  CURR. YEAR,  NUM. MOE  and  NUM. ALT  as 
integer  variables 

define  LAMDA ,  LAMDA1 ,  LAMDA2 ,  LAMDA3 ,  SUM. LCC,  SUM. LCC. 1,  SUM. LCC. 2, 
SUM. LCE,  SUM. LCE. 1,  SUM. LCE. 2,  MAX. LCE,  MIN. LCC,  MOE. VALUE, 
SLOPE1 ,  SLOPE2,  SLOPE3 ,  LI,  L2 ,  L3 ,  BUDGET,  AMC,  AOC,  DCY , 

SVY,  O.MC.P1,  0.MC.P2,  O.OC,  0. OC.EF,  O.DC,  O. DC.EF,  O.SV, 

O. SV.EF,  MAX. MAC,  MAX.AAC,  AAC.CF,  INT . RATE  and  INT.DF  as 
real  variables 

define  MAX.CF  as  a  double  variable  "  Double  means  double  precision, 
define  APPROACH,  ECHO. OH,  CASE. NAME  and  ALT. NAME  as  text  variables 


end 
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1  main 

2 

3  call  INITIALIZE 

4 

5  if  APPROACH  =  "Lagrangian" 

6  call  LAGRANGIAN 

7  else 

8  call  BEST. PATH 

9  call  SUM.LCC.LCE 

10  always 

11 

12  call  RESULTS. PRINT 

13 

14  stop 

15 

1 6  end 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


routine  INITIALIZE 

open  1  for  input,  name  is  "IN. DAT" 
use  1  for  input 

open  2  for  output,  name  is  "OUT. DAT" 
use  2  for  output 

read  CASE. NAME 

read  N.SYS,  N.ALT,  NUM.M0E,  CURR. YEAR,  INT . RATE 
read  ECHO. ON 


create  every  SYS  and  ALT 


N . YEAR  =  1 

for  I  =  1  to  N.SYS 

do 

read  MIN.LIFE(I)  and  MAX.LIFE(I)  • •  Computing  number  of  years  to 
N . YEAR  -  N . YEAR  +  MAX.LIFE(I)  • •  be  covered  in  the  analysis, 

loop 

create  every  YEAR 

for  I  =  1  to  N.ALT 

read  ALT. NAME (I)  and  LEAD . TIME { I ) 
read  DISCOUNT. FACTOR ( 1) 

CUM.DF(l)  -  DISCOUNT. FACTOR(l) 
for  I  ■=  2  to  N .  YEAR 
do 

read  DISCOUNT . FACTOR ( I ) 

CUM.DF(I)  -  CUM.DF(I-l)  *  DISCOUNT . FACTOR ( I ) 
loop 

read  APPROACH  ' '  Determines  the  type  of  solution  to  be  generated 

read  BUDGET  ' '  Used  only  in  the  Lagranglan  approach 

1 1  Allocating  memory  to  arrays 
reserve  MAX . PRED. NODES ( * )  as  N.SYS 

reserve  AC(*,*),  MC.P1(*,*),  MC,P2(*,*),  OC(*,*),  OC.EF(*,*), 

DC (*,*),  DC . EF ( * , * ) ,  SV(*,*),  SV . EF ( * , * )  as  N.ALT  by  N . YEAR 
reserve  MOE(*,*,*)  and  MOE . DF (*,*,* )  as  N.ALT  by  NUM.MOE  by  N . YEAR 
reserve  HOE. NAME (*),  MOE.WT(«),  O.MOE(*),  O.MOE.DF(*)  as  NUM.MOE 
reserve  MIN . MOE . REQ ( * , * )  as  NUM.MOE  by  N . YEAR 

call  CALC1 . NODE  ' 1  Structuring  the  network. 


create  every  NODE  and  ARC 
reserve  CRIT.FUNC(*)  as  N . NODE 
reserve  EP<*)  as  (N . NODE  +  1) 

reserve  LCE(*),  LCC(*),  LIM.M0E(«)  and  L(*)  as  N . ARC 


call  CALC2 . NODE  '• 

i  i 

call  CALC. ALT. CODE  ' • 
read  0.MC.P1,  0.MC.P2, 


Computing  the  year  &  system  of  each  node 
&  the  start  index  of  arcs  entering  it. 

Assigning  code  of  alternatives. 

0.0C,  0.0C.EF,  0 . DC ,  0.DC.EF,  O.SV,  O.SV.EF 
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56  for  I  •=  1  to  N . ALT 

57  for  J  =  1  to  N . YEAR 

58  read  AC(I, J) ,  MC.P1(I,J),  MC.P2(I,J),  OC(I,J),  OC.EF(I,J), 

59  DC ( I ,  J ) ,  DC .  E  F  ( I ,  J ) ,  SV(I,J)  and  SV.EF(I,J) 

60  for  J  -  1  to  NUM.MOE 

61  read  MOE . NAME ( J ) ,  MOE.WT(J),  O.MOE(J)  and  O.MOE.DF(J) 

62  for  J  =  1  to  NUM.MOE 

63  for  K  =  1  to  N . YEAR 

64  read  MIN .MOE . REQ( J , K) 

65  for  I  =  1  to  N . ALT 

66  for  J  =  1  to  NUM.MOE 

67  for  K  =  1  to  N. YEAR 

68  read  MOE(I,J,K)  and  MOE. DF(I , J , K) 

69 

70  call  CALC. ARC  "  Computes  tall  node,  LCC  and  LCE  for  each  arc. 

71 

72  print  3  line  with  CASE. NAME  thus 

73  Dynamic  End-Of-Life-Cycle  Planning  Model 

74 

75  Case  of  ************************************** 

76 

77  skip  1  line 

78  print  4  line  with  APPROACH  and  BUDGET  thus 

79  Approach  -  ******************  Budget  =  $*********.** 

80 

81  Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 

82  ignored  in  the  other  approaches. 


83 

84 

if  (  ECHO. 

ON  »  "On" 

85 

call  ECHO. INPUT 

86 

87 

skip  1 

line 

88 

print  9 

line  thus 

89  Candidate  Arcs  For  Solution  Path 

90 

91  Note  -  LCC  is  made  to  be  unfavorably  high  if  the  minimum  reguirement 

92  level  for  any  component  MOE  of  LCE  has  not  been  met.  The  limiting  MOE 

93  is  given  as  the  highest  index  of  the  MOE(s)  that  failed  to  meet  the 

94  minimum  requirement  level (s)  during  the  life-span  covered  by  the  arc. 

95 

96  Arc  Alt.  Tail  Head  Life-Cycle  Limiting 

97  No.  Code  Node  Node  Effectiveness  Cost  MOE 

98 

99  skip  1  line 

100  for  I  =  2  to  N . NODE 

101  for  K  =  EP ( 1 )  to  (  EP ( 1+ 1 )  -  1  ) 

102  print  1  line  with  K,  ALT . CODE (K) ,  TAIL(K) ,  I,  LCE(K) ,  LCC(K) 

103  and  LIM.MOE(K)  thus 

} 0  4  ***  ***  ***  ***  ******.****  *****,**  *** 

105  always 

106 

107  return 

108 

109  end 
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routine  CALC1.N0DE 

• '  Structuring  the  network  by  determining  the  number  of 
1 ’  nodes  and  arcs  in  each  system 

BASE . YEAR ( 1 )  =  CURR. YEAR  +  MIN.LIFE(l) 

MAX. PRED. NODES (1)  =  1 

NUM . NODES ( 1 )  =  MAX . LIFE ( 1 )  -  MIN.LIFE(l)  +  1 
NUM . ARCS ( 1 )  =  NUM. NODES (1) 

N . NODE  =  NUM. NODES (1)  +  1 
N . ARC  =  NUM . ARCS ( 1 ) 
for  I  =  2  to  N.SYS 
do 

J  =  I  -  1 

BASE . YEAR ( I )  =  BASE . YEAR ( J )  +  MIN.LIFE(I) 

MAX. PRED. NODES (I)  =  NUM . NODES ( I- 1 ) 

NUM . NODES ( I )  =  NUM . NODES ( J)  +  MAX.LIFE(I)  -  MIN.LIFE(I) 
NUM. ARCS (1)  =  N . ALT  *  NUM. NODES (J)  * 

(  MAX . LIFE ( I )  -  MIN . LIFE ( I )  +  1  ) 

N . NODE  =  N . NODE  +  NUM.NODES(I) 

N . ARC  =  N . ARC  +  NUM . ARCS ( I ) 
loop 

return 


end 
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1  routine  CALC2 . NODE 

2 

3  1 •  Computes  for  each  node  information  such  as  the  year  and  system 

4  ' 1  it  belongs  to,  as  well  as  the  lead  index  (  contained  in  entry 

5  1 1  point  array  EP  )  of  the  arcs  entering  the  node 

6 

7  TRANS. YEAR(l)  =  CURR . YEAR 

8  S YS . CODE ( 1 )  =  0 

9  EP ( 1 )  =  0 

10  K  =  2 

11  JJ  =  1 

12  for  J  =  1  to  NUM. NODES (1) 

13  do 

14  TRANS. YEAR (K)  =  BASE. YEAR (1)  +  J  -  1 

15  SYS . CODE (K)  =  1 

16  EP ( K)  =  EP(K-l)  +  1 

17  K  =  K  +  1 

18  loop 

19  for  I  =  2  to  N . SYS 

20  do 

21  for  J  =  1  to  NUM. NODES (I) 

22  do 

23  TRANS. YEAR (K)  =  BASE. YEAR (I)  +  J  -  1 

24  S YS . CODE ( K)  =  I 

25  EP ( K)  =  EP(K-l)  +  JJ 

26  M  =  MAX . LIFE ( I )  -  MIN.LIFE(I)  +  1 

27  JJ  =  N.ALT*min.f ( J , MAX . PRED . NODES ( I ) ,M, (NUM. NODES ( I ) -J+l ) ) 

28  K  =  K  +  1 

29  loop 

30  loop 

31  EP ( N . NODE+ 1 )  »  EP (N . NODE)  +  N .ALT 

32 

33  return 

34 

3  5  end 
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1  routine  CALC. ALT. CODE 

2 

3  1 '  Assigns  alternative  codes  to  each  arc. 

4 

5  for  K  =  1  to  NUM. NODES (1) 

6  ALT. CODE (K)  =  0 

7  for  I  =  (  NUM . NODES ( 1 )  +  2  )  to  N . NODE 

8  for  K  =  EP ( I )  to  (  EP ( 1+ 1 )  -  1  ) 

9 

10  ALT . CODE (K)  =  (  mcd.f  (  K  -  EP(I)  , 

11 

12  return 

13 

14  end 


N .ALT  )  +  1  ) 
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routine  CALC. ARC 

1 1  Determines  the  tall  node,  life-cycle  cost  &  effectiveness 
• '  associated  with  each  arc. 

1 '  Initialising  LCE  &  Limiting  HOE  arrays 
for  K  =  1  to  N . ARC 
do 

LCE ( K)  =  0 
LIM.MOE(K)  «=  0 
loop 

' '  Computations  for  the  present  system  (  system  1  ) . 

' '  All  arcs  here  have  the  alternative  code  0. 
for  K  -  1  to  NUM. NODES (1) 
do 

TAIL(K)  *»  1 

M  -  TRANS. YEAR(K+1)  -  CURR . YEAR 

LCC(K)  -  O . DC  *  exp . f (  O.DC.EF  *  M  )  -  O.SV  *  exp.f(  O.SV.EF  *  H  ) 
for  I  -  1  to  NUM.MOE 
do 

MOE. VALUE  =  O.MOE(I)  "  LCC  is  made  very  high 

if  (  MOE. VALUE  <  MIN . MOE . REQ ( I , M+ 1 )  )  ••  if  any  MOE  require- 

LCC(K)  =  9999999  ' '  ment  level  is  not  met 

LIM.MOE(K)  =  I 
always 

' '  Note  :  Utility  fn.  is  Value  /  Max. Value 

if  (  I  -  1  ) 

LCE (K)  =  LCE (K)  +  MOE.WT(I)  *  MOE. VALUE  /  180 
else 

LCE ( K)  -  LCE(K)  +  MOE.WT(I)  *  MOE. VALUE  /  100 
always 
loop 

for  J  «  l  to  M 
do 

AMC  =  O.MC.P1  *  O.MC.P2  *  (  J  **  (  O.MC.P2  -  1  )  ) 

AOC  =  O.OC  *  exp. f (  O.OC.EF  *  J  ) 

LCC { K)  =  LCC(K)  +  (  AMC  +  AOC  )  *  CUM.DF(J) 
for  I  -  1  to  NUM.MOE 
do 

MOE. VALUE  «  O.MOE(I)  *  exp.f(  O.MOE.DF(I)  *  J  ) 
if  (  MOE. VALUE  <  MIN . MOE . REQ ( I , J+ 1 )  )  •'  Checking  if  MOE 

LCC ( K)  «  9999999  "  req.  were  met. 

LIM.MOE(K)  «=  I 
always 

' '  Note  :  Utility  fn.  is  Value  /  Max. Value 

if  (  I  -  1  ) 

LCE (K)  «=  LCE ( K)  +  MOE.WT(I)  *  MOE. VALUE  /  180 
else 

LCE(K)  -  LCE (K)  +  MOE.WT(I)  *  MOE. VALUE  /  100 
always 
loop 
loop 
loop 
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Computations  for  later  systems,  represented  by 
• '  arcs  with  alternative  code  higher  than  0. 
for  S  =  2  to  N.SYS 

for  J  =  3  to  N . NODE 
with  SYS . CODE ( J)  =  S 
for  I  =  2  to  N . NODE 
with  S YS . CODE ( I )  =  S  -  1 
do 

M  =  TRANS. YEAR(J)  -  TRANS . YEAR ( I )  ' *  Number  of  years  used, 

if  (  (  M  >=  MIN . LIFE ( S )  )  and  (  M  <«  MAX.LIFE(S)  )  ) 

for  N  *=  1  to  N .  ALT 
do 

TAIL(K)  =  I  11  Start  node  of  the  arc  K. 

DCY  =  DC ( N ,  I )  * 

exp . f (  DC . EF (N , I )  *  M  )  ' '  Disposal  cost  S 

SVY  =  SV(N,I)  *  1 '  salvage  value  in  the 

exp. f (  SV . EF (N , I )  *  M  )  "  last  year. 

YR  =  TRANS. YEAR(I)  -  CURR . YEAR  "  No.  of  years 

"  from  the  present. 

LCC(K)  =  AC (M , I )  *  CUM . DF ( YR+ 1 )  + 

(  DCY-SVY  )  *  CUM. DF ( YR+ 1 ) 
for  E  =  1  to  NUM.MOE 
do 

MOE. VALUE  =  MOE (N, E, YR+ 1) 

"  Note  :  Utility  fn.  is  Value  /  Max. Value 
if  (  E  =  1  ) 

LCE(K)  -  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  180 
else 

LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  100 
always 

if  (  MOE. VALUE  <  MIN . MOE . REQ ( E , YR+ 1 )  ) 

I,CC(K)  =  9999999  '  1  Checking  if  MOE 

LIM.MOE(K)  =  E  "  req.  were  met. 

always 
loop 

for  Y  =  1  to  M 
do 

AMC  =  MC. PI (M, I)  *  MC . P2 (N, I )  *  ' '  Annual  maint. 

(  Y  *  *  (  MC . P2 ( N , I }  -  1  )  )  ' ’  £  operating 
AOC  >=  OC  (N ,  I  >  *  exp .  f  (  OC.EF(N.I)  *  Y  )  "  costs. 

LCC(K)  =  LCC(K)  +  (  AMC  +  AOC  )  * 

CUM. DF(YRIM) 
for  E  »  1  to  NUM.MOE 
do 

MOE. VALUE  =  MOE ( N , E , YR+ 1 )  * 

exp.f  (  MOE . DF (N , E , YR+ 1 )  *  Y  ) 

1  1  Note  :  Utility  fn.  is  Value  /  Max. Value 
if  (  E  “  1  ) 

LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  180 
else 

LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOF. VALUE  /  100 
always 

if  (  MOE. VALUE  <  MIN .MOE. REQ (E , YR+Y+1 )  ) 
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LCC(K)  =  9999999 
LIM.MOE(K)  =  E 
always 
loop 
loop 

compute  MAX.LCE  as  the 
maximum  of  LCE(K) 
if  (  LCC(K)  <  9999999  ) 
compute  MIN.LCC  as  the 
minimum  of  LCC(K) 

always 
K  =  K  +  1 
loop 
always 
loop 


return 


end 


1 •  Checking  if  MOE 
1 '  req.  were  met. 


"  MAX.LCE  &  MIN.LCC  are 
•  1  needed  for  computing 
' '  the  upper  bound  on  the 
1 '  Lagrangian  multiplier. 
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1  routine  ECHO. INPUT 

2 

3  ' 1  Echo  of  important  inputs 

4 

5  skip  1  line 

6  print  1  line  thus 

7  Sys.Code  Min. Life  Max. Life  Base.Yr.  No. Nodes  No. Arcs 


8 

For  I  =  1  to  N . SYS 

9 

print  1  line  with 

I,  MIN . LIFE ( I ) , 

MAX . LIFE ( I ) ,  BASE. YEAR(I) , 

10 

NUM . NODES ( I )  ,  ! 

NUM . ARCS ( I )  thus 

11 

***  *** 

*  *  *  *  *  *  ★ 

****  **** 

12 

1  3 

skip  1  line 

14 

print  1  line  thus 

15 

Alternative  Name 

Alt. Code 

Lead  Time 

16 

For  1  =  1  to  N . ALT 

17 

print  1  line  with 

ALT. NAME (I),  I, 

LEAD-TIME (I)  thus 

18 

******************** 

*** 

*  *  * 

19 

20 

skip  1  line 

21 

print  2  line  thus 

22 

Measure  of  Effectiveness 

Present  System 

23 

Code  Name  Weightage 

Level  Exp. RateOf Change 

24 

for  I  =  1  to  NUM.MOE 

25 

print  1  line  with 

I,  MOE.NAME(I), 

MOE . WT (I) ,  o . MOE ( I ) ,  O . MOE . DF ( I ) 

26 

*  *  *  A************** 

*  .  *  *  * 

******.***  ****.*** 

27 

28 

return 

29 

30 

end 

109 


1  routine  LAGRAHGIAN 

2 

3  ' '  Solution  by  Lagrangian  relaxation  method 

4 

5  ' 1  The  starting  lower  bound  for  the  Lagrangian  multiplier  lamda  is  0. 

6  1 1  The  upper  bound  is  determined  such  that  it  is  barely  sufficient  to 

7  •'  ensure  that  the  corresponding  slope  will  be  positive.  Variables 

8  ' '  may  be  overwhelmed,  even  when  double  precision  is  used,  if  lamda 

9  1 1  or  the  cost  values  are  too  high. 

10 

11  LAMDA 1  =  0 

12  LAMDA2  =  (  N . NODE  *  MAX.LCE  /  MIN.LCC  )  +  1 

13 

14  skip  1  line 

15  print  3  line  thus 

16  Iterations  of  the  Lagrangian  Relaxation  Method 

17 

18  Lamdal  Lamda]  Lamda2  Slopel  Slope!  Slope2  LI  L3  L2 

19  Skip  1  line 

20 

21  1 1  Iterations  by  intercept  of  supporting  planes  method. 

22 

23  'Iterate' 

24 

25  LAMDA  =  LAMDA1 

26  call  BEST. PATH 

27  call  SUM.LCC.LCE 

28  SUM.LCC.l  =  SUM.LCC 

29  SUM.LCE.l  =  SUM.LCE 

30  SLOPE1  =  -1  *  (  SUM.LCC.l  -  BUDGET  ) 

31  LI  -  SUM.LCE.l  -  LAMDA 1  *  (  SUM.LCC.l  -  BUDGET  ) 

32 

33  if  (  SLOPE1  >  0  ) 

34  skip  1  line 

35  print  2  line  thus 

36  The  budget  is  sufficient  for  cost  to  be  no  longer  a  constraint, 

37  ie.  the  solution  is  that  of  maximum  effectiveness  over  time. 

38  return 

39  always 

40 

41  LAMDA  =  LAMDA2 

42  call  BEST. PATH 

43  call  SUM.LCC.LCE 

44  SUM.LCC. 2  =  SUM.LCC 

45  SUM.LCE. 2  =  SUM.LCE 

46  SLOPE2  =  -1  *  (  SUM.LCC. 2  -  BUDGET  ) 

47  L2  -  SUM.LCE. 2  -  LAMDA2  *  (  SUM.LCC. 2  -  BUDGET  ) 

48 

49  if  (  SL0PE2  <  0  ) 

50  skip  1  line 

51  print  4  line  with  BUDGET  thus 

52  Ho  solution  is  possible  because  of  one  of  the  following  reasons: 

53  (1)  the  given  budget  of  ********.**  is  insufficient 

54  (2)  one  or  more  MOE  requirement  level  is  too  high 

55  (3)  cost  values  are  too  high  and  needs  to  be  scaled  down. 
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56  stop 

57  always 

58 

59  if  (  (  SLOFE1  =  SLOPE2  )  and  (  LI  =  L2  )  ) 

60  skip  1  line 

61  print  1  line  thus 

62  The  cost  values  are  too  high  for  computation  &  should  be  scaled  down. 

63  stop 

64  always 

65 

66  *  1  Intercept  of  the  line  segments  at  lamdal  &  lamda2  yields  lamda3. 

67  LAMDA3  ■=  (  SUM.LCE.l  -  SUM .  LCE .  2  )  /  (  SUM.LCC.l  -  SUM.LCC.2  ) 

68 

69  LAMDA  =  LAM DA 3 

70  call  BEST. PATH 

71  call  SUM . LCC . LCE 

72  SLOPE3  ■=  -1  *  (  SUM.  LCC  -  BUDGET  ) 

73  L3  =  SUM. LCE  -  LAMDA 3  *  (  SUM. LCC  -  BUDGET  ) 

74 

75  print  1  line  with  LAMDA1 ,  LAM DA 3 ,  LAMDA 2 ,  SLOPE1,  SLOPE3,  SLOPE2 , 

76  LI,  L3 ,  L2  thus 

77  ***.***  ****.***  ****.***  ****,**  ****.**  ******  ****.*  ****,*  ****.* 

78 

79  • '  The  minima  of  the  piece-wise  linear  Lagrangian  function  L 

80  11  is  reached,  ie.  the  solution  is  found,  when  the  intercept 

81  1 1  lamda3  is  equal  to  any  of  the  two  bounds 

82  if  (  IAMDA3  =  LAMDA2  )  ' '  Minima  is  found  6  slope  there  is  positive 

83  return  '  '  ie.  minima  itself  '  he  optimal  solution 

84  always 

85  if  (  LAMDA3  =  LAMDA 1  )  ' '  Minima  is  found  but  slope  is  negative. 

86  LAMDA  =  LAMDA 2  * '  Lamda  lies  on  the  solution  line  segment, 

87  call  BEST. PATH  1 1  which  is  the  one  with  positive  slope  next 

88  call  SUM. LCC. LCE  ' '  to  the  minima. 

89  SLOPE3  =  -1  *  (  SUM. LCC  -  DUDGFT  5 

90  L3  =  SUM.  LCE  -  IAMDA3  *  (  SUM. LCC  -  bl'DGET  , 

91  print  1  line  with  LAMDA 1 ,  LAM DA 3 ,  LAMDA 2 ,  SLOPE1,  SLOPE3 ,  SL0TE2, 

92  LI,  L3,  L2  thus 

93  ***.***  ****,***  **«*.***  ****.**  ****,**  *«*«.**  ***«,*  ****.*  ****.* 

94  return 

95  always 

96 

97  '  '  Lamda3  replaces  the  redundant  bound  for  the  next  iteration. 

98  if  (  SLOPE3  >=  0  ) 

99  let  LAMDA 2  ■=  LAMDA 3 

100  else 

101  let  LAMDA 1  =  LAM  DA  3 

102  always 

103 

104  go  to  'Iterate' 

105 

106  end 
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routine  BEST. PATH 

"  Forward  dynamic  programming  with  criteria  function  dependent  on 
1 •  the  given  approach.  Provisions  are  made  to  ensure  that  the  solu- 
' '  tion  path  has  sufficient  lead  time  for  each  alternative  chosen. 

CRIT.FUNC(l)  =  0 

for  I  =  2  to  N . NODE  "  Initialising  crit.  func.  with  a  very 

CRIT.FUNC(I)  -  -99999999  "  low  value  as  0  is  too  high  a  default 

1 '  value  to  prevent  adverse  selections. 

for  I  »  2  to  N . NODE 
do 

MAX.CF  =  -999999999 

for  J  =  EP ( I )  to  ( EP ( 1  +  1) -1 ) 

do 

• ’  Lagrangian  relaxation  approach  to  obtain  the 
"  max.  eff.  for  the  given  budget, 
if  (  APPROACH  =  "Lagrangian"  ) 
if  (  LCC(J)  <  9999999  ) 

L(J)  -  (  LCE(J)  -  LAMDA  *  LCC(J)  ) 
else 

L(J)  =  -9999999 
always 
always 

'  *  The  approach  of  maximising  the  cost-eff. 

' 1  ratio  while  meeting  given  requirements  for 
1 1  effectiveness.  The  budget  is  ignored, 
if  (  APPROACH  •=  "Cost-Eff .Ratio"  ) 
if  (  LCC(J)  <  9999999  ) 

L(J)  -  LCE(J)  /  (  LCC(J)  -  MIN.LCC  +  1  ) 
else 

L(J)  =  -9999999 
always 
always 

' '  The  approach  of  satisfying  effectiveness 
' '  requirements  at  minimum  cost. 

'  '  Yields  the  lower  bound  on  budget, 
if  (  APPROACH  =  "Satisficing"  ) 

L(J)  =  -1  *  LCC(J) 
always 

1  '  The  approach  of  maximising  effectiveness  without 
'  '  concern  for  budget.  Yields  the  upper  bound  on  budget, 
if  (  APPROACH  ■=  "Max.  Eff."  ) 
if  (  LCC(J)  <  9999999  ) 

L(J)  =  LCE(J) 
else 

L(J)  =  -9999999 
always 
always 

1 '  Preventing  the  selection  of  an  arc  that  cannot 
' '  be  used  because  of  insufficient  lead  time. 
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56 

If  (  ALT. CODE ( J)  >  0  ) 

57 

if  (  BEST . PR ED. NODE (TAIL( J ) )  -  0  ) 

58 

L(J)  =  -9999999 

59 

else 

60 

K  =  TRANS. YEAR(TAIL(J) )  - 

61 

TRANS . YEAR (BEST. PRED.NODE (TAIL( J) ) ) 

62 

if  (  K  <  LEAD. TIME ( ALT . CODE ( J ) ) 

) 

63 

L(J)  =  -9999999 

64 

always 

65 

always 

66 

always 

67 

68 

if  (  CRIT. FUNC (TAIL(J ) )  >  -9999999  ) 

69 

if  (  CRIT . FUNC (TAIL (J ) )  +  L(J)  )  > 

MAX.CF 

70 

MAX.CF  •=  CRIT.  FUNC(TAIL(J)  )  +  L(J) 

71 

CRIT . FUNC ( I )  =  MAX.CF 

72 

BEST. FRED. NODE (I)  -  TAIL(J) 

73 

BEST. PRED. ARC ( I )  -  J 

74 

always 

75 

always 

76 

77 

loop 

78 

loop 

79 

80 

1 '  The  terminal  node  with  the  highest  value 

in  the 

81 

' '  criteria  function  yields  the  solution  path  that 

82 

' '  is  optimal  for  the  given  approach. 

83 

MAX.CF  =  -999999999 

84 

MAX. MAC  =  -999999999 

85 

MAX.AAC  =  -999999999 

86 

for  I  =  1  to  N . NODE 

87 

with  SYS  .  CODE  ( I )  =>  SYS. CODE(N. NODE) 

88 

do 

89 

if  (  APPROACH  =  "Satisficing"  ) 

90 

N  =  TRANS. YEAR (I)  -  CURR . YEAR  +  1 

91 

INT.DF  =  (  (  INT . RATE  /  100  )  +  1  ) 

**  (  N  - 

92 

AAC.CF  =  (  INT.DF  *  ( INT. RATE/100) 

)  /  (  INT 

93 

if  (  (  CRIT . FUNC ( I )  /  N  )  >  MAX. MAC 

) 

94 

MAX. MAC  =  CRIT . FUNC ( I )  /  N 

95 

ARG.MAC  -  I 

96 

always 

97 

if  (  (  CRIT . FUNC ( I )  *  AAC.CF  )  >  MAX.AAC  ) 

98 

MAX.AAC  -  CRIT . FUNC ( I )  *  AAC.CF 

99 

ARG.AAC  -  I 

100 

always 

101 

always 

102 

if  (  CRIT . FUNC ( I )  >=  MAX.CF  ) 

103 

MAX.CF  =  CRIT . FUNC ( I ) 

104 

APGMAX.CF  =  I 

105 

always 

106 

loop 

107 

108 

return 

109 

110 

end 
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1  routine  SUM.LCC.LCE 

2 

3  ' 1  Totals  up  the  life-cycle  cost  &  effectiveness 

4  ’ 1  for  the  solution  path. 

5 

6  I  =  ARGMAX.CF 

7  SUM.LCC  =  0 

8  SUM.LCE  =  0 

9  until  1=1 

10  do 

11  if  (  BEST . PRED . ARC ( I )  =  0  ) 

12  print  2  lines  thus 

13  All  paths  violates  either  lead  time  or  MOE  level  requirements 

14  or  both.  No  feasible  solution  can  be  found. 

15  stop 

16  always 

17  SUM.LCC  =  SUM.LCC  +  LCC ( BEST. PRED. ARC ( I ) ) 

18  SUM.LCE  =  SUM.LCE  +  LCE ( BEST . PRED . ARC ( I ) ) 

19  I  =  BEST. PRED. NODE(I) 

20  loop 

21 

22  return 

23 

24  end 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
4  1 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


routine  RESULTS . PRINT 

' '  Print-Out  of  Results 


skip  1  line 
print  6  line  thus 

Network  of  Possible  Solutions  in  Nodal  Form 

Note  -  The  criterion  function  of  a  node  will  be  very  low  if  all  paths 
leading  to  it  fail  to  satisfy  lead  time  req.  arid/or  minimum  MOE  levels. 

No.  Sys.  Year  EP  Crit.Fn.  Best. Pred. Node  Best . Pred . Arc  Best. Alt. 
skip  1  line 

print  1  line  with  SYS.CODE(l),  TRANS . YEAR ( 1 ) ,  EP(1)  and 
CRIT.FUNC(l)  thus 
1  *  *  *  **  *  aaaa.aaa 

for  I  -  2  to  N . NODE 

print  1  line  with  I,  SYS.CODE(I),  TRANS . YEAR ( I ) ,  EP(I), 
CRIT.FUNC(I) ,  BEST. PRED. NODE ( I )  ,  BEST . PRED . ARC ( I )  and 
ALT. CODE (max. f (1, BEST . PRED. ARC ( I ) ) )  thus 

***  ***  A*  A  AAAA.AAA  AAA  AAA  AAA 


if  (  CRIT . FUNC ( ARGMAX . CF)  <=  -99999  ) 
print  1  line  thus 

Warning  :  Lead-time  requirements  cannot  be  met. 
always 

if  <  SUM.LCC  >=  9999999  ) 
print  1  line  thus 

Warning  s  Minimum  MOE  level  requirements  cannot  be  met. 
always 


if  (  APPROACH  -  "Satisficing"  ) 
skip  1  line 
print  2  line  thus 

Minimum  Average  Annual  Cost  Solution  Path 
(  Time  value  of  money  is  ignored,  ie.  interest 
skip  1  line 
print  2  line  thus 

Transition  Alternative  Life-Cycle 

Node  Year  System  Employed  Cost 

skip  1  line 
I  =  ARG.MAC 
until  1=1 
do 


'  '  There  are  two  possible  solutions 
'  '  for  the  satisficing  approach. 


rate  is  assumed  to  be  0  ) 


Life-Cycle 

Effectiveness 


Criterion 
Func .Value 


print  1  line  with  I,  TRANS . YEAR ( I ) ,  SYS.CODE(I), 

ALT. C0DE(BEST. PRED. ARC(I) ) ,  LCC ( BEST. PRED . ARC ( I ) ) , 

LCE ( BEST . PRED . ARC { I ) ) ,  CRIT.FUNC(I)  thus 

AAA  AA  *  ***  AAAAA.AAA  AAAAA.AAA  AA AAA .AAA 

I  =  BEST . PRED . NODE ( I ) 
loop 

skip  1  line 

print  2  line  with  I NT . RATE  thus 
Minimum  Annualised  Cost  Solution  Path 

(  Monetary  interest  assumed  constant  at  aaa.aaa  \  per  annum  ) 


56  sKip  1  line 

57  print  2  line  thus 

58  Transition  Alternative  Life-Cycle  Life-Cycle  Criterion 

59  Node  Year  System  Employed  Cost  Effectiveness  Func. Value 

60  skip  1  line 

61  I  =  ARG.AAC 

62  until  1=1 

63  do 

64  print  1  line  with  I,  TRANS . YEAR ( I ) ,  SYS.CODE(I), 

65  ALT. CODE (BEST. FRED. ARC(I) ) ,  LCC ( BEST . PRED . ARC ( I ) ) , 

66  LCE ( BEST . PRED . ARC ( I )  )  ,  CRIT.FUNC(I)  thus 

67  *«*  **  *  ***  *****.***  *****.***  *****_*** 

68  I  =  BEST. PRED.NODE(I) 

69  loop 

70  return 

71  always 

72 

73  skip  1  line 

74  print  4  line  thus 

75  Solution  Path 

76 

77  Transition  Alternative  Life-Cycle  Life-Cycle  Criterion 

78  Node  Year  System  Employed  Cost  Effectiveness  Func. Value 

79 

80  skip  1  line 

81  I  =  ARGMAX.CF 

82  until  1=1 

8  3  do 

84  print  1  line  with  1,  TRANS . YEAR ( I ) ,  SYS.CODE(I), 

85  ALT. CODE (BEST. PRED. ARC(I) ) ,  LCC ( BEST . FRED . ARC ( I ) ) , 

86  LCE ( BEST . PRED . ARC ( I ) ) ,  CRIT.FUNC(I)  thus 

87  ***  **  *  ***  *****.***  ********  ******** 

88  I  =  BEST. PRED. NODE(I) 

89  loop 

90 

91  skip  1  line 

92  if  (  (  CRIT . FUNC ( ARGMAX . CF )  >  -99999  )  and  (  SUM. LCC  <  9999999  )  ) 

93  print  2  line  with  SUM. LCC  and  SUM. LCE  thus 

94  The  total  life-cycle  cost  and  effectiveness  for  the  optimal  solution 

95  are  *******,**  and  ******,***  respectively. 

96  always 

97 

98  skip  1  line 

99  if  (  APPROACH  =  "Lagrangian"  ) 

100  print  1  line  with  LAM DA 3  thus 

101  Optimal  value  of  the  Lagrangian  multiplier  is  *****.***** 

102  skip  1  line 

103  always 

104 

105  return 

106 

107  end 
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APPENDIX  B:  SAMPLE  INPUT  FILES 


Stationary . Environment 
323  1989  0 
Off 
0  5 
5  10 
5  10 

Alt. A  0 
Alt.B  0 

1111111111111 

1111111111111 

l.agrangian 

1200 

30  1.3  7  .05  0  0  20  - . 2 
100  30  1.3  7  .05  0  0  50  - . 2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  - . 2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  - . 2 

100  30  1.3  7  .05  0  0  50  - . 2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

100  30  1.3  7  .05  0  0  50  -.2 

300  15  1.3  4  .05  0  0  90  - .  2' 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1 . 3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 
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300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 
300  15  1.3  4 


05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 
05  0  0  90  -.2 


Range(nm)  0.4  104.6  -.03 
Accuracy(%)  0.3  87.2  -.02 
Availability (%)  0.3  93.5  -.03 


95 

95  95  95 

95 

95 

95 

95 

95 

95  95  95 

95 

95 

95  95  95 

95 

95 

95 

95 

95 

95  95  95 

95 

80 

80  8 C 

)  80 

80 

80 

80 

80 

80 

80  80  80 

80 

80 

80  80  80 

80 

80 

80 

80 

80 

80  80  80 

80 

70 

70  70  70 

70 

70 

70 

70 

70 

70  70  70 

70 

70 

70  70  70 

70 

70 

70 

70 

70 

70  70  70 

70 

120 

-.03 

i  120 

-  . 

03 

120 

i  — 

.  03 

120 

-.03 

120 

- . 

03 

120 

-.03 

120 

-.03 

i  120 

03 

120 

i  - . 

.03 

120 

-.03 

120 

03 

120 

-.03 

120 

-.03 

120 

03 

120 

i  - . 

.03 

120 

-.03 

120 

03 

120 

-.03 

120 

-.03 

120 

-. 

03 

120 

i  - 

.03 

120 

-.03 

120 

03 

120 

-.03 

120 

-.03 

120 

- . 

03 

95 

-.02 

95  - 

.02 

95 

_  ' 

02 

95 

-.02 

95  - 

-.02 

95 

— 

.  02 

95 

-.02 

95  - 

.02 

95 

02 

95 

-.02 

95  -.02 

95 

- 

.02 

95 

-.02 

95  - 

.02 

95 

02 

95 

-.02 

95  -.02 

95 

- 

.02 

95 

-.02 

95  - 

.02 

95 

02 

95 

-.02 

95  - 

-.02 

95 

- 

.02 

95 

-.02 

95  - 

.02 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  -.03 

95 

— , 

.03 

95 

-.03 

95  - 

.03 

95 

“  • 

03 

95 

-.03 

95  - 

-.03 

95 

.03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

-.03 

95 

.03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

-.03 

95 

.03 

95 

-.03 

95  - 

.03 

140 

-.03 

140 

- . 

03 

140 

— , 

,03 

140 

-.03 

140 

- . 

03 

140 

-.03 

140 

-.03 

140 

-. 

03 

140 

.  03 

140 

-.03 

140 

03 

140 

-.03 

140 

-.03 

140 

-. 

03 

140 

,  03 

140 

-.03 

140 

-. 

03 

140 

-.03 

140 

-.03 

140 

03 

140 

,  03 

140 

-.03 

140 

-. 

03 

140 

-.03 

140 

-.03 

140 

— . 

03 

95 

-.02 

95  - 

.  02 

95 

- . 

02 

95 

-.02 

95  - 

-.02 

95 

- , 

,02 

95 

-.02 

95  - 

.02 

95 

02 

95 

-.02 

95  - 

-.02 

95 

,02 

95 

-.02 

95  - 

.02 

95 

-. 

02 

95 

-.02 

95  - 

-.02 

95 

02 

95 

-.02 

95  - 

.02 

95 

02 

95 

-.02 

95  - 

■.02 

95 

02 

95 

-.02 

95  - 

.02 

95 

-.03 

95  - 

.03 

95 

_  _ 

03 

95 

-.03 

95  - 

-.03 

95 

_  m 

03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

.  03 

95 

-. 

03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

.03 

95 

03 

95 

-.03 

95  - 

.03 

95 

03 

95  -.03 

95  - 

.03 

1 18 


Inflation. (.Interest. Rates 
323  1989  5 
Off 
0  5 
5  10 
5  10 

Alt. A  0 
Alt . B  0 


.95  . 

95  .95 

i  .95 

.95 

.  95 

.95  . 

95  .95 

i  .95 

.95 

.95 

.95  . 

95  .95 

-  .95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95 

Satisficing 

900 

30.0 

1 . 3 

7.0 

.  05 

0 

0 

20 

-.2 

100.0 

31.2 

1 . 3 

7.3 

.05 

0 

0 

50.0 

-.2 

104.0 

32.4 

1.3 

7.6 

.05 

0 

0 

52 . 0 

-.2 

108.2 

33.7 

1 . 3 

7.9 

.  05 

0 

0 

54.1 

-.2 

112.5 

35.1 

1.3 

8.2 

.  05 

0 

0 

56 . 2 

-.2 

117.0 

36.5 

1 . 3 

8 . 5 

.05 

0 

0 

58 . 5 

-  .  2 

121.7 

38 . 0 

1.3 

8.9 

.05 

0 

0 

60.8 

-.2 

126.5 

39.5 

1.3 

9.2 

.05 

0 

0 

63.2 

-.2 

131.6 

41.1 

1.3 

9.6 

.05 

0 

0 

65.8 

-.2 

136.9 

42.7 

1.3 

10.0 

.  05 

0 

0 

68.4 

-.2 

142.3 

44.4 

1.3 

10.4 

.  05 

0 

0 

71.2 

-.2 

148.0 

46.2 

1.3 

10.8 

.  05 

0 

0 

74.0 

-.2 

153.9 

48.0 

1.3 

11.2 

.  05 

0 

0 

77.0 

-.2 

160.1 

50.0 

1.3 

11 . 7 

.05 

0 

0 

80.1 

-.2 

166.5 

52.0 

1.3 

12.1 

.05 

0 

0 

83.3 

-.2 

173.2 

54.0 

1 . 3 

12.6 

.  05 

0 

0 

86.6 

-.2 

180.1 

56.2 

1.3 

13.1 

.05 

0 

0 

90.0 

-.2 

187.3 

53.4 

1 . 3 

13.6 

.05 

0 

0 

93.6 

-.2 

194.8 

60.8 

1.3 

14.2 

.05 

0 

0 

97.4 

-.2 

202 . 6 

63 . 2 

1.3 

14.7 

.05 

0 

0 

101.3 

-.2 

210.7 

65.7 

1.3 

15.3 

.05 

0 

0 

105.3 

-.2 

219.1 

68.4 

1 . 3 

16.0 

.05 

0 

0 

109.6 

-.2 

227 . 9 

71.1 

1 . 3 

16.6 

.  05 

0 

0 

113.9 

-.2 

237.0 

73.9 

1.3 

17.2 

.05 

0 

0 

118.5 

-.2 

246.5 

76.9 

1.3 

17.9 

.  05 

0 

0 

123.2 

-.2 

256.3 

80.0 

1,3 

18.7 

.05 

0 

0 

128.2 

-.2 

266.5 

83.2 

1.3 

19.4 

.05 

0 

0 

133.3 

-.2 

300.0 

15.6 

1 . 3 

4 . 2 

.05 

0 

0 

90.0 

-.2 

312.0 

16.2 

1 . 3 

4 . 3 

.  05 

0 

0 

93 . 6 

-.2 

324.5 

16.9 

1.3 

4 . 5 

.  05 

0 

0 

97.3 

-.2 

337.5 

17.5 

1 . 3 

4 . 7 

.  05 

0 

0 

101.2 

-.2 

351.0 

18 . 2 

1 . 3 

4.9 

.05 

0 

0 

105.3 

-.2 

365.0 

19.0 

1.3 

5.1 

.05 

0 

0 

109 . 5 

-.2 

379.6 

19.7 

1.3 

5.3 

.  05 

0 

0 

113.9 

-.2 

394.8 

20.5 

1 . 3 

5.5 

.  05 

0 

0 

118.4 

-.2 

410.6 

21.3 

1.3 

5.7 

.05 

0 

0 

123.2 

-.2 

427.0 

22.2 

1.3 

5.9 

.05 

0 

0 

128.1 

-.2 

444  .  1 

23.1 

1 . 3 

6.2 

.  05 

0 

0 

133.2 

-.2 

461.8 

24.0 

1 . 3 

6.4 

.05 

0 

0 

138.6 

-.2 

480.3 

25.0 

1  .  3 

6.7 

.  05 

0 

0 

144.1 

-.2 

499 . 5 

26.0 

1 . 3 

6.9 

.05 

0 

0 

149.9 

-.2 

519.5 

27.0 

1.3 

7.2 

.05 

0 

0 

155.9 

-.2 

Ilf) 


540.3 

28.1 

1.3 

7.5  .05 

0 

0 

162.1 

-.2 

561.9 

29.2 

1.3 

7.8  .05 

0 

0 

168.6 

-.2 

584.4 

30.4 

1.3 

8.1  .05 

0 

0 

175.3 

-.2 

607 . 7 

31.6 

1.3 

8.4  .05 

0 

0 

182.3 

-.2 

632.1 

32.9 

1.3 

8.8  .05 

0 

0 

189.6 

-.2 

657 . 3 

34.2 

1.3 

9.1  .05 

0 

0 

197.2 

-.2 

683.6 

35.5 

1.3 

9.5  .05 

0 

0 

205.1 

-.2 

711.0 

37.0 

1.3 

9.9  .05 

0 

0 

213.3 

-.2 

739.4 

38 . 4 

1.3 

10.3  .05 

0 

0 

221.8 

-.2 

769.0 

40.0 

1.3 

10.7  .05 

0 

0 

230.7 

-.2 

799.8 

41.6 

1.3 

11.1  .05 

0 

0 

239.9 

-.2 

Range (nm) 

0.4  104 

.6 

-.03 

Accuracy (%) 

0.3  87 

.2 

-.02 

Availability (%) 

0.3  93 

.  5 

-.03 

95 

95  95 

95  95 

95  95 

95 

95 

95  95 

95 

95 

95 

95  95 

95  95 

95  95 

95 

95 

95  95 

95 

95 

80 

80  80 

80  80 

80  80 

80 

80 

80  80 

80 

80 

80 

80  80 

80  80 

80  80 

80 

80 

80  80 

80 

80 

70 

70  70 

70  70 

70  70 

70 

70 

70  70 

70 

70 

70 

70  70 

70  70 

70  70 

70 

70 

70  70 

70 

70 

120 

-.03 

120  -. 

03  120 

_ 

.  03 

120  - 

.03 

120 

-.03 

120 

-.03 

120 

-.03 

120  -. 

03  120 

- 

.  03 

120  - 

.  03 

120 

-.03 

120 

-.03 

120 

-.03 

120  -. 

03  120 

- 

.  03 

120  - 

.03 

120 

-.03 

120 

-.03 

120 

-.03 

120  -. 

03  120 

- 

.03 

120  - 

.03 

120 

-.03 

120 

-.03 

120 

-.03 

120  -. 

03 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.  02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.  02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.02 

95 

-.02 

95  -.02 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.  03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.03 

95 

-.03 

95  -.03 

140 

-.03 

140  -. 

03  140 

_ 

.03 

140  - 

.03 

140 

-.03 

140 

-.03 

140 

-.03 

140  -. 

03  140 

- 

.  03 

140  - 

.03 

140 

-.03 

140 

-.03 

140 

-.03 

140  -. 

03  140 

- 

.  03 

140  - 

.  03 

140 

-.03 

140 

-.03 

140 

-.03 

140  -. 

03  140 

- 

.  03 

140  - 

.  03 

140 

-.03 

140 

-.03 

140 

-.03 

140  -. 

03 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.  02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.02 

95 

-.02 

95  -.02 

95  -. 

02 

95 

-.02 

95 

-.02 

95  - 

.  02 

95 

-.02 

95  -.02 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.  03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.  03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-.03 

95  - 

.  03 

95 

-.03 

95  -.03 

95  -. 

03 

95 

-.03 

95 

-  .  03 

95  - 

.03 

95 

-.03 

95  -.03 

120 
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Techno logical. Change 
323  1989  5 
Off 
0  5 
5  10 
5  10 

Alt. A  0 
Alt.B  0 


.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.  .95 

.95 

.95 

.95  . 

95  .95 

i  .95 

.95 

.95 

.95  . 

95  .95 

■  .95 

.95 

.95 

.95  . 

95 

Satisficing 

650 

30.0 

1.3 

7.0 

.  05 

0 

0 

20 

-.2 

100.0 

30.3 

1.3 

7.3 

.05 

0 

0 

50.0 

-.2 

102 . 0 

30.6 

1 .  3 

7 . 6 

.05 

0 

0 

51.0 

-.2 

104 . 4 

30.9 

1.3 

7.9 

.05 

0 

0 

52.2 

-.2 

106.1 

31.2 

1 . 3 

8 . 2 

.05 

0 

0 

53.0 

-.2 

108.2 

31.5 

1.3 

8.5 

.  05 

0 

0 

54.1 

-.2 

110.4 

31.8 

1.3 

8.9 

.  05 

0 

0 

55.2 

-.2 

112.6 

32.2 

1.3 

9.2 

.05 

0 

0 

56 . 3 

-.2 

114.9 

32 . 5 

1.3 

9.6 

.05 

0 

0 

57 . 4 

-.2 

117.2 

32 . 8 

1.3 

10.0 

.05 

0 

0 

58.6 

-.2 

119.5 

33.1 

1.3 

10.4 

.  05 

0 

0 

59.7 

-.2 

121.9 

33.5 

1.3 

10.8 

.05 

0 

0 

60.9 

-.  2 

124 . 3 

33.8 

1 . 3 

11.2 

.05 

0 

0 

62.1 

-.2 

126.8 

34 . 1 

1.3 

11.7 

.  05 

0 

r 

63 . 4 

-.2 

129.4 

34.5 

1.3 

12.1 

.O'" 

,• 

J 

64.7 

-.2 

131.9 

34.8 

1.3 

12.6 

.  05 

0 

65.9 

-.2 

134.6 

35.2 

1.3 

13.1 

.  '  , 

0 

0 

67.3 

-.2 

137.3 

35.5 

1 . 3 

13.6 

.05 

0 

0 

68.6 

-.2 

140.0 

35.9 

1.3 

14  .2 

.05 

0 

0 

70.0 

-.2 

142.8 

36.2 

1.3 

14.7 

.05 

0 

0 

71.4 

-.2 

145.7 

36.6 

1.3 

15.3 

.05 

0 

0 

72.8 

-.2 

148.6 

37.0 

1.3 

16.0 

.05 

0 

0 

74.3 

-.2 

151 . 6 

37.3 

1 . 3 

16.6 

.05 

0 

0 

75.8 

-.2 

154 . 6 

37  7 

1 . 3 

17.2 

.05 

0 

0 

77 . 3 

-.2 

157.7 

38 . 1 

1 . 3 

17.9 

.05 

0 

0 

78 . 8 

-.2 

160.8 

38.5 

1 . 3 

18.7 

.05 

0 

0 

80.4 

-.2 

164 . 1 

38.9 

1.3 

19.4 

.05 

0 

0 

82 . 0 

-.2 

300.0 

15.6 

1 . 3 

4 . 2 

.  05 

0 

0 

90.0 

-  .  2 

303 . 0 

16.2 

1.3 

4 . 3 

.  05 

0 

0 

90.9 

-.2 

306.0 

16.9 

1.3 

4.5 

.05 

0 

0 

91.8 

-.2 

309. 1 

17 . 5 

1 . 3 

4.7 

.05 

0 

0 

92 . 7 

-.2 

312.2 

18 . 2 

1 . 3 

4.9 

.  05 

0 

0 

93 . 7 

-.2 

315.3 

19.0 

1  .  3 

5.1 

.  05 

0 

0 

94 . 6 

-.2 

318.5 

19.7 

1.3 

5.3 

.05 

0 

0 

95.5 

-.2 

321.6 

20.5 

1 . 3 

5.5 

.05 

0 

0 

96.5 

-.2 

324.9 

21.3 

1 . 3 

5.7 

.  05 

0 

0 

97 . 5 

-.2 

328 . 1 

22 . 2 

1  .  3 

5.9 

.05 

0 

0 

98 . 4 

-.2 

331.4 

23.1 

1  .  3 

6 . 2 

.  05 

0 

0 

99 . 4 

-.2 

334.7 

24.0 

1  .  3 

6 . 4 

.  05 

0 

0 

100.4 

-.  2 

338.0 

25.0 

1 .  3 

6.7 

.05 

0 

0 

101.4 

-.2 

341.4 

26.0 

1  .  3 

6.9 

.05 

0 

0 

102.4 

-.2 

344.8 

27 . 0 

1.3 

7.2 

.05 

0 

0 

103 . 5 

-.2 
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348.3 

28.1 

1.3 

7.5  .05 

0 

0 

104.5 

-.2 

351 . 8 

29.2 

1 . 3 

7.8  .05 

0 

0 

105.5 

-.2 

355. 3 

30.4 

1.3 

8.1  .05 

0 

0 

106.6 

-.2 

358.8 

31 . 6 

1.3 

8.4  .05 

0 

0 

107.7 

-.2 

362.4 

32.9 

1.3 

8.8  .05 

0 

0 

108.7 

-.2 

366 . 1 

34.2 

1.3 

9.1  .05 

0 

0 

109.8 

-.2 

369.7 

35.5 

1.3 

9.5  .05 

0 

0 

110.9 

-.2 

373.4 

37.0 

1.3 

9.9  .05 

0 

0 

112.0 

-.2 

377 . 1 

38.4 

1.3 

10.3  .05 

0 

0 

113.1 

-.2 

380.9 

40.0 

1.3 

10.7  .05 

0 

0 

114.3 

-.2 

384.7 

41.6 

1.3 

11.1  .05 

0 

0 

115.4 

-.2 

Range (nm) 

0.4  104. 

6 

-.03 

Accuracy (%) 

0.3  87. 

2 

-.02 

Availability (%) 

0.3  93. 

5 

-.03 

95  95 

95 

95  95 

95  95 

95  95 

95 

95  95 

95 

95  95 

95 

95  95 

95  95 

95  95 

95 

95  95 

95 

80  80 

80 

80  80 

80  80 

i  80  80 

80 

80  80 

80 

80  80 

80 

80  80 

80  80 

80  80 

80 

80  80 

80 

70  70 

70 

70  70 

70  70 

i  70  70 

70 

70  70 

70 

70  70 

70 

70  70 

70  70 

i  70  70 

70 

70  70 

70 

120  - 

.03 

123  - 

03  128  -.03 

134 

-.03 

140 

-.03 

145 

-.03 

150  - 

.03 

154  - 

03  158  -.03 

162 

!  -.03 

166 

-.03 

169 

-.03 

172 

.03 

175  - 

03  178  -.03 

179 

i  -.03 

180 

-.03 

180 

-.03 

180 

.  03 

180  - 

03  180  -.03 

180 

1  -.03 

180 

-.03 

180 

-.03 

180  - 

.03 

180  - 

03 

95.0 

-.02 

95.1 

-.02 

95.2  - 

.02 

95.3  - 

.02 

95.4 

-.02 

95.5 

-.02 

95.6 

-.02 

95.7 

-.02 

95.8  - 

.02 

95.9  - 

.02 

96.0 

-.02 

96.1 

-.02 

96.2 

-.02 

96.3 

-.02 

96.4  - 

.02 

96.4  - 

.02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02 

96.8  - 

.02 

96.8  - 

.02 

96.9 

-.02 

96.9 

-.02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  - 

.03 

95.3  - 

.03 

95.4 

-.03 

95.5 

-.03 

95.6 

-.03 

95.7 

-.03 

95.8  - 

.  03 

95.9  - 

.03 

96.0 

-.03 

96.1 

-.03 

96.2 

-.03 

96.3 

-.03 

96.4  - 

.  03 

96.4  - 

.03 

96.5 

-.03 

96.6 

-.03 

96.7 

-.03 

96.7 

-.03 

96.8  - 

.  03 

96.8  - 

.03 

96.9 

-.03 

96.9 

-.03 

96.9 

-.03 

97.0 

-.03 

140  - 

.03 

144  -. 

03  147  -.03 

150 

i  -.03 

153 

-.03 

155 

-.03 

158  - 

.  03 

160  -. 

03  162  -.03 

164 

-.03 

165 

-.03 

167 

-.03 

168  - 

.  03 

169  -. 

03  170  -.03 

171 

-.03 

172 

-.03 

173 

-.03 

173  - 

.03 

174  -. 

03  175  -.03 

175 

-.03 

176 

-.03 

176 

-.03 

176  - 

.  03 

177  -. 

03 

95.0 

-.02 

95.1 

-.02 

95.2  - 

.  02 

95.3  - 

.02 

95.4 

-.02 

95.5 

-.02 

95.6 

-.02 

95.7 

-.02 

95.8  - 

.02 

95.9  - 

.  02 

96.0 

-.02 

96. 1 

-.02 

96.2 

-.02 

96.3 

-.02 

96.4  - 

.  02 

96.4  - 

.  02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02 

96.8  - 

.02 

96.8  - 

.  02 

96 . 9 

-.02 

96.9 

-.02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  - 

.  03 

95.3  - 

.  03 

95 . 4 

-.03 

95.5 

-  .  03 

95.6 

-.03 

95.7 

-.03 

95.8  - 

.03 

95.9  - 

.03 

96 . 0 

-.03 

96.1 

-.03 

96. 2 

-.03 

96 . 3 

-.03 

96.4  - 

.  03 

96.4  - 

.03 

96.5 

-.03 

96.6 

-  .  03 

96.7 

-.03 

96.7 

-.03 

96.8  - 

.  03 

96.8  - 

.03 

96.9 

-.03 

96.9 

-.03 

96.9 

-.03 

97.0 

-.03 
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Requirement . Escalation 
323  1989  5 
Off 
0  5 
5  10 
5  10 

Alt. A  0 
Alt. B  0 


.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.  95 

.  95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95 

Lagrangian 

700 

30.0 

1.3 

7 . 0 

.05 

0 

0 

20 

-.2 

100.0 

30.3 

1.3 

7.3 

.  05 

0 

0 

50.0 

-.2 

102.0 

30 . 6 

1 . 3 

7 . 6 

.05 

0 

0 

51.0 

-.2 

104.4 

30.9 

1 . 3 

7 . 9 

.  05 

0 

0 

52.2 

-.2 

106.1 

31.2 

1 . 3 

8 . 2 

.05 

0 

0 

53 . 0 

-.2 

108.2 

31.5 

1 . 3 

8 . 5 

.05 

0 

0 

54 . 1 

-.2 

110.4 

31.8 

1 . 3 

8 . 9 

.  05 

0 

0 

55 . 2 

-  .  2 

112.6 

32.2 

1 . 3 

9.2 

.05 

0 

0 

56.3 

-.2 

114.9 

32.5 

1 . 3 

9 . 6 

.  05 

0 

0 

57 . 4 

-.2 

117.2 

32.8 

1 . 3 

10.0 

.05 

0 

0 

58.6 

-.2 

119.5 

33 . 1 

1  .  3 

10.4 

.05 

0 

0 

59.7 

-.2 

121.9 

33.5 

1 . 3 

10.8 

.05 

0 

0 

60.9 

-.2 

124 . 3 

33.8 

1.3 

11.2 

.05 

0 

0 

62.1 

-.2 

126.8 

34 . 1 

1 . 3 

11.7 

.05 

0 

0 

63.4 

-.2 

129.4 

34 . 5 

1.3 

12 . 1 

.  05 

0 

0 

64 . 7 

-.2 

131.9 

34.8 

1 . 3 

12.6 

.05 

0 

0 

65.9 

-.2 

134.6 

35.2 

1.3 

13 . 1 

.  05 

0 

0 

67 . 3 

-.2 

137.3 

35.5 

1 .  3 

13 . 6 

.05 

0 

0 

68 . 6 

-.2 

140.0 

35.9 

1.3 

14.2 

.05 

0 

0 

70.0 

-.2 

142.8 

36.2 

1.3 

14.7 

.05 

0 

0 

71.4 

-.2 

145.7 

36.6 

1.3 

15.3 

.05 

0 

0 

72.8 

-.2 

148.6 

37.0 

1.3 

16.0 

.05 

0 

0 

74.3 

-.2 

151.6 

37.3 

1.3 

16.6 

.05 

0 

0 

75.8 

-.2 

154.6 

37.7 

1.3 

17 . 2 

.05 

0 

0 

77 . 3 

-.2 

157.7 

38.1 

1.3 

17.9 

.05 

0 

0 

78.8 

-.2 

160.8 

38.5 

1 . 3 

18.7 

.05 

0 

0 

80.4 

-.2 

164 . 1 

38.9 

1.3 

19.4 

.05 

0 

0 

82 . 0 

-.2 

300.0 

15.6 

1  .  3 

4 . 2 

.  05 

0 

0 

90 . 0 

-  .  2 

303.0 

16.2 

1  .  3 

4  .  3 

.  05 

0 

0 

90.9 

-.2 

306 . 0 

16.9 

1 . 3 

4 . 5 

.05 

0 

0 

91.8 

-.  2 

309 . 1 

17.5 

1 . 3 

4.7 

.  05 

0 

0 

92 . 7 

-.2 

312.2 

18 . 2 

1.3 

4.9 

.  05 

0 

0 

93 . 7 

-.2 

315.3 

19.0 

1.3 

5.1 

.05 

0 

0 

94.6 

-.2 

318.5 

19.7 

1 . 3 

5.3 

.  05 

0 

0 

95.5 

-  .  2 

321.6 

20.5 

1.3 

5.5 

.05 

0 

0 

96.5 

-.2 

324.9 

21.3 

1 . 3 

5.7 

.05 

0 

0 

97.5 

-.2 

328  .  1 

22 . 2 

1.3 

5.9 

.  05 

0 

0 

98 . 4 

-  .  2 

331.4 

23.1 

1  .  3 

6.2 

.  05 

0 

0 

99 . 4 

-  .  2 

334.7 

24 .0 

1 . 3 

6.4 

.05 

0 

0 

100.4 

-.2 

338.0 

25.0 

1  .  3 

6.7 

.  05 

0 

0 

101.4 

-  .  2 

341.4 

26.0 

1 . 3 

6.9 

.  05 

0 

0 

102 . 4 

-.2 

344.8 

27.0 

1 . 3 

7.2 

.  05 

0 

0 

103 . 5 

-.2 

123 


348.3  28.1  1.3  7.5  .05  0  0  104.5  -.2 

351.8  29.2  1.3  7.8  .05  0  0  105.5  -.2 

355.3  30.4  1.3  8.1  .05  0  0  106.6  -.2 

358.8  31.6  1.3  8.4  .05  0  0  107.7  -.2 

362.4  32.9  1.3  8.8  .05  0  0  108.7  -.2 

366.1  34.2  1.3  9.1  .05  0  0  109.8  -.2 

369.7  35.5  1.3  9.5  .05  0  0  110.9  -.2 

373.4  37.0  1.3  9.9  .05  0  0  112.0  -.2 

377.1  38.4  1.3  10.3  .05  0  0  113.1  -.2 

380.9  40.0  1.3  10.7  .05  0  0  114.3  -.2 

384.7  41.6  1.3  11.1  .05  0  0  115.4  -.2 

Range(nm)  0.4  104.6  -.03 

Accuracy(%)  0.3  87.2  -.02 

Availability (%)  0.3  93.5  -.03 

95  95  95  100  100  100  110  110  110  110  110  110  110 

110  110  110  110  110  110  110  110  110  110  110  110  110 


80  80 

80 

80  80 

80  80 

80  80 

80 

80  80 

80 

80  80 

80 

80  80 

80  80 

80  80 

80 

80  80 

80 

70  70 

70 

72  72 

72  72 

75  75 

75 

75  75 

75 

75  75 

75 

75  75 

75  75 

75  75 

75 

75  75 

75 

120  - 

.  03 

123  -. 

03  128  -.03 

134 

-.03 

140 

-.03 

145 

-.03 

150  - 

.  03 

154  -. 

03  158  -.03 

162 

-.03 

166 

-.03 

169 

-.03 

172  - 

.  03 

175  -. 

03  178  -.03 

179 

-.03 

180 

-.03 

180 

-.03 

180  - 

.  03 

180  -. 

03  180  -.03 

180 

-.03 

180 

-.03 

180 

-.03 

180  - 

.  03 

180  -. 

03 

95.0 

-.02 

95.1 

-.02 

95.2  - 

.02 

95.3 

-.02 

95.4 

-.02 

95.5 

-.02 

95 . 6 

-.02 

95.7 

-.02 

95.8  - 

.  02 

95.9 

-.02 

96 . 0 

-.02 

96. 1 

-.02 

96.2 

-.02 

96.3 

-.02 

96.4  - 

.  02 

96.4 

-.02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02 

96.8  - 

.  02 

96.8 

-.02 

96.9 

-.02 

96.9 

-.02 

96.9 

-.02 

97 . 0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  - 

.03 

95.3 

-.03 

95.4 

-.03 

95.5 

-.03 

95.6 

-.03 

95.7 

-.03 

95.8  - 

.  03 

95 . 9 

-.03 

96.0 

-.03 

96 . 1 

-.03 

96.2 

-.03 

96.3 

-.03 

96.4  - 

.03 

96.4 

-.03 

96.5 

-.03 

96.6 

-.03 

96.7 

-.03 

96.7 

-.03 

96.8  - 

.03 

96.8 

-.03 

96.9 

-.03 

96.9 

-.03 

96.9 

-.03 

97.0 

-.03 

140  - 

.03 

144  -. 

03  147  -.03 

150 

-.03 

153 

-.03 

155 

-.03 

158  - 

.  03 

160  -. 

03  162  -.03 

164 

-.03 

165 

-.03 

167 

-.03 

168  - 

.  03 

169  -. 

03  170  -.03 

171 

-.03 

172 

-.03 

173 

-.03 

173  - 

.  03 

174  -. 

03  175  -.03 

175 

-.03 

176 

-.03 

176 

-.03 

176  - 

.  03 

177  -. 

03 

95.0 

-.02 

95.1 

-.02 

95.2  - 

.  02 

95.3 

-.02 

95 . 4 

-.02 

95.5 

-.02 

95.6 

-.02 

95.7 

-.02 

95.8  - 

.02 

95.9 

-.02 

96.0 

-.02 

96. 1 

-.02 

96.2 

-.02 

96.3 

-.02 

96.4  - 

.02 

96.4 

-.02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02 

96.8  - 

.  02 

96.8 

-.02 

96.9 

-.02 

96.9 

-.02 

96.9 

-.02 

97 . 0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  - 

.  03 

95.3 

-.03 

95.4 

-.03 

95.5 

-.03 

95.6 

-.03 

95.7 

-.03 

95.8  - 

.03 

95.9 

-.03 

96.0 

-.03 

96. 1 

-.03 

96.2 

-.03 

96.3 

-.03 

96.4  - 

.  03 

96.4 

-.03 

96.5 

-.03 

96.6 

-.03 

96.7 

-.03 

96.7 

-.03 

96.8  - 

.  03 

96.8 

-.03 

96.9 

-.03 

96.9 

-.03 

96.9 

-.03 

97 . 0 

-.03 
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Upgrading . Scheduling 
323  1989  5 
Off 
0  5 
5  30 
5  10 

Alt. A  0 
Alt.B  0 


.95  . 

95  .95 

.95 

.95 

.  95 

.  95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95 

Cost- 

Eff .Ratio 

780.08 

30.0 

1.3 

7.0 

.05 

0 

0 

20 

-.2 

100.0 

30. 3 

1 . 3 

7 . 3 

.  05 

0 

0 

50.0 

-.2 

82 . 0 

30.6 

1 . 3 

7 . 6 

.  05 

0 

0 

51.0 

-.2 

104.4 

30.9 

1 . 3 

7 . 9 

.  05 

0 

0 

52 . 2 

-  .  2 

106.1 

31.2 

1 . 3 

8 . 2 

.05 

0 

0 

53 . 0 

-.2 

86. 2 

31.5 

1.3 

8.5 

.05 

0 

0 

54 . 1 

-.2 

110.4 

31.8 

1 . 3 

8 . 9 

.  05 

0 

0 

55 . 2 

-.2 

112.6 

32 . 2 

1.3 

9.2 

.05 

0 

0 

56 . 3 

-.2 

90.9 

32 . 5 

1  .  3 

9.6 

.  05 

0 

0 

57.4 

-.2 

117.2 

32.8 

1 . 3 

10.0 

.  05 

0 

0 

58 . 6 

-.2 

119.5 

33.1 

1.3 

10.4 

.05 

0 

0 

59.7 

-.2 

95.9 

33.5 

1 . 3 

10.8 

.05 

0 

0 

60.9 

-.2 

124.3 

33 . 8 

1.3 

11.2 

.05 

0 

0 

62 . 1 

-.2 

126.8 

34 . 1 

1.3 

11.7 

.  05 

0 

0 

63.4 

-.2 

101.4 

34 . 5 

1  .  3 

12 . 1 

.  05 

0 

0 

64 . 7 

-.2 

131.9 

34 . 8 

1 . 3 

12.6 

.05 

0 

0 

65.9 

-.2 

134.6 

35.2 

1 . 3 

13 . 1 

.05 

0 

0 

67.3 

-.2 

107 . 3 

35.5 

1.3 

13.6 

.05 

0 

0 

68.6 

-.2 

140.0 

35.9 

1 . 3 

14 . 2 

.05 

0 

0 

70.0 

-.2 

142.8 

36.2 

1.3 

14.7 

.  05 

0 

0 

71.4 

-.2 

313.7 

36.6 

1.3 

15.3 

.  05 

0 

0 

72.8 

-.2 

148.6 

37.0 

1 . 3 

16.0 

.  05 

0 

0 

74 . 3 

-.2 

151.6 

37.3 

1.3 

16.6 

.  05 

0 

0 

75.8 

-.2 

120.6 

37 . 7 

1 . 3 

17 . 2 

.05 

0 

0 

77 . 3 

-.2 

157.7 

38 . 1 

1 . 3 

17.9 

.  05 

0 

0 

78 . 8 

-.2 

160.8 

38 . 5 

1 . 3 

18.7 

.  05 

0 

0 

80.4 

-.2 

128 . 1 

38.9 

1.3 

19 . 4 

.  05 

0 

0 

82 . 0 

-.2 

300.0 

15.6 

1  .  3 

4 . 2 

.  05 

0 

0 

90 . 0 

-  .  2 

303.0 

16.2 

1  .  3 

4  .  3 

.  05 

0 

0 

90.9 

-  .  2 

306.0 

16.9 

1 . 3 

4 . 5 

.  05 

0 

0 

91.8 

-  .  2 

309.1 

17.5 

1  .  3 

4 . 7 

.  05 

0 

0 

92 . 7 

-.2 

312.2 

18.2 

1  .  3 

4 . 9 

.05 

0 

0 

93.7 

-.  2 

315.3 

19.0 

1  .  3 

5 . 1 

.  05 

0 

0 

94 . 6 

-  .  2 

318.5 

19.7 

1  .  3 

5.3 

.  05 

0 

0 

95.5 

-.2 

321.6 

20.5 

1 . 3 

5.5 

.  05 

0 

0 

96 . 5 

-.2 

324.9 

21.3 

1  .  3 

5.7 

.  05 

0 

0 

97.5 

-.2 

328.1 

22 . 2 

1  .  3 

5.9 

.05 

0 

0 

98 . 4 

-  .  2 

331.4 

23.1 

1  .  3 

6.2 

.05 

0 

0 

99 . 4 

-.2 

3  34.7 

24.0 

1  .  3 

6.4 

.05 

0 

0 

100.4 

-.2 

338.0 

25.0 

1  .  3 

6 . 7 

.  05 

0 

0 

101.4 

-.2 

341.4 

26.0 

1.3 

6.9 

.05 

0 

0 

102.4 

-.2 

34  4.8 

27.0 

1 .  3 

7 . 2 

.  05 

0 

0 

103 . 5 

-  .  2 

125 


348.3 

28.1 

1.3 

7.5  .05 

0 

0 

104.5 

-.2 

351.8 

29 . 2 

1 . 3 

7.8  .05 

0 

0 

105.5 

-.2 

355.3 

30.4 

1.3 

8.1  .05 

0 

0 

106.6 

-.2 

358.8 

31.6 

1.3 

8.4  .05 

0 

0 

107 . 7 

-.2 

362 . 4 

32.9 

1.3 

8.8  .05 

0 

0 

108.7 

-.2 

366.1 

34 . 2 

1.3 

9.1  .05 

0 

0 

109 . 8 

-.2 

369.7 

35.5 

1.3 

9.5  .05 

0 

0 

110.9 

-.2 

373.4 

37 . 0 

1.3 

9.9  .05 

0 

0 

112.0 

-.2 

377 . 1 

38 . 4 

1.3 

10.3  .05 

0 

0 

113 . 1 

-.2 

380.9 

40.0 

1.3 

10.7  .05 

0 

0 

114 . 3 

-.2 

384.7 

41.6 

1 . 3 

11.1  .05 

0 

0 

115.4 

-.2 

Range (nm) 

0.4  104 

.6 

-.03 

Accuracy (%) 

0.3  87 

.  2 

-.02 

Availability (%) 

0.3  93 

.5 

-.03 

95 

95 

95  100 

100 

100  110 

110  110 

110 

110  : 

110  110 

110  110  110  110 

110 

110  110 

110  110 

110 

no  : 

110  110 

80  80 

80 

80  80 

80  80 

80  80 

80 

80  80 

80 

80  80 

80 

80  80 

80  80 

80  80 

80 

80  80 

80 

70  70 

70 

72  72 

72  72 

75  75 

75 

75  75 

75 

75  75 

75 

75  75 

75  75 

75  75 

75 

75  75 

75 

120  - 

.03 

123  -. 

03  128  -.03 

134 

-.03 

140 

-.03 

145 

-.03 

150  - 

.  03 

154  -. 

03  158  -.03 

162 

-.03 

166 

-.03 

169 

-.03 

172  - 

.  03 

175  -. 

03  178  -.03 

179 

-.03 

180 

-.03 

180 

-.03 

180  - 

.03 

180  -. 

03  180  -.03 

180 

i  -.03 

180 

-.03 

180 

-.03 

180  - 

.03 

180  -. 

03 

95.0 

-.02 

95.1 

-.02 

95.2  -. 

02 

95.3 

-.02 

95.4 

-.02 

95.5 

- 

.  02 

95.6 

-.02 

95.7 

-.02 

95.8  -. 

02 

95.9 

-.02 

96.0 

-.02 

96.1 

- 

.  02 

96.2 

-.02 

96.3 

-.02 

96.4  -. 

02 

96.4 

-.02 

96.5 

-.02 

96.6 

- 

.  02 

96.7 

-.02 

96.7 

-.02 

96.8  -. 

02 

96.8 

-.02 

96.9 

-.02 

96.9 

- 

.  02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  -. 

03 

95.3 

-.03 

95.4 

-.03 

95.5 

- 

.  03 

95.6 

-.03 

95.7 

-.03 

95.8  -. 

03 

95.9 

-.03 

96 . 0 

-.03 

96.1 

- 

.03 

96.2 

-.03 

96.3 

-.03 

96.4  -. 

03 

96.4 

-.03 

96 . 5 

-.03 

96.6 

- 

.  03 

96.7 

-.03 

96.7 

-.03 

96.8  -. 

03 

96.8 

-.03 

96 . 9 

-.03 

96.9 

- 

.  03 

96.9 

-.03 

97.0 

-.03 

140  - 

.  03 

144  -. 

03  147  -.03 

150 

'  -.03 

153 

-.03 

155 

-.03 

158  - 

.  03 

160  -. 

03  162  -.03 

164 

-.03 

165 

-.03 

167 

-.03 

168  - 

.  03 

169  -. 

03  170  -.03 

171 

-.03 

172 

-.03 

173 

-.03 

173  - 

.  03 

174  -. 

03  175  -.03 

175 

-.03 

176 

-.03 

176 

-.03 

176  - 

.  03 

177  -. 

03 

95.0 

-.02 

95.1 

-.02 

95.2  -. 

02 

95. 3 

-.02 

95.4 

-.02 

95.5 

- 

.  02 

95.6 

-.02 

95.7 

-.02 

95.8  -. 

02 

95.9 

-.02 

96 . 0 

-.02 

96.1 

- 

.  02 

96.2 

-.02 

96.3 

-.02 

96.4  -. 

02 

96.4 

-.02 

96.5 

-.02 

96.6 

- 

.  02 

96.7 

-.02 

96.7 

-.02 

96.8  -. 

02 

96.8 

-.02 

96.9 

-.02 

96.9 

- 

.  02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03 

95.2  -. 

03 

95.3 

-.03 

95.4 

-.03 

95.5 

- 

.  03 

95.6 

-.03 

95.7 

-.03 

95.8  -. 

03 

95.9 

-.03 

96 . 0 

-  .  03 

96. 1 

- 

.  03 

96.2 

-  .  03 

96 . 3 

-.03 

96.4  -. 

03 

96.4 

-  .03 

96.5 

-.03 

96 . 6 

- 

.  03 

96.7 

-  .  03 

96.7 

-.03 

96.8  -. 

03 

96.8 

-.03 

96 . 9 

-.03 

96.9 

- 

.  03 

96.9  -.03  97.0  -.03 
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Development . Lead .Time 
323  1989  5 
On 
0  5 
5  10 
5  10 

Alt. A  1 
Alt . B  0 


.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95  .95 

.95 

.95 

.95 

.95  . 

95 

Lagrangian 

634 . 52 

30.0 

1.3 

7.0 

.05 

0 

0 

20 

-  .  2 

100.0 

30.3 

1 . 3 

7.3 

.05 

0 

0 

50.0 

-.2 

82.0 

30.6 

1.3 

7.6 

.05 

0 

0 

51.0 

-.2 

104.4 

30.9 

1 . 3 

7.9 

.05 

0 

0 

52 . 2 

-.2 

106 . 1 

31.2 

1.3 

8.2 

.05 

0 

0 

53 . 0 

-.2 

86 . 2 

31.5 

1.3 

8.5 

.  05 

0 

0 

54 . 1 

-.2 

110.4 

31.8 

1 . 3 

8.9 

.05 

0 

0 

55.2 

-.2 

112.6 

32 . 2 

1 . 3 

9.2 

.  05 

0 

0 

56.3 

-.2 

90.9 

32.5 

1.3 

9.6 

.05 

0 

0 

57.4 

-.2 

117.2 

32.8 

1.3 

10.0 

.  05 

0 

0 

58.6 

-.2 

119.5 

33.1 

1.3 

10.4 

.05 

0 

0 

59.7 

-.2 

95 . 9 

33 . 5 

1.3 

10.8 

.  05 

0 

0 

60.9 

-.2 

124 . 3 

33.8 

1.3 

11.2 

.05 

0 

0 

62.1 

-.2 

126.8 

34.1 

1.3 

11.7 

.  05 

0 

0 

63 . 4 

-.2 

101.4 

34.5 

1.3 

12.1 

.05 

0 

0 

64.7 

-.2 

131.9 

34.8 

1.3 

12.6 

.05 

0 

0 

65.9 

-.2 

134.6 

35.2 

1.3 

13.1 

.05 

0 

0 

67.3 

-.2 

107.3 

35.5 

1.3 

13.6 

.05 

0 

0 

68.6 

-.2 

140.0 

35.9 

1.3 

14.2 

.05 

0 

0 

70.0 

-.2 

142.8 

36.2 

1.3 

14.7 

.05 

0 

0 

71.4 

-.2 

113.7 

36.6 

1.3 

15.3 

.05 

0 

0 

72.8 

-.2 

148.6 

37.0 

1.3 

16.0 

.05 

0 

0 

74 . 3 

-.2 

151 . 6 

37 . 3 

1 . 3 

16.6 

.  05 

0 

0 

75.8 

-.2 

120.6 

37.7 

1.3 

17.2 

.  05 

0 

0 

77.3 

-.2 

157.7 

38.1 

1.3 

17 . 9 

.05 

0 

0 

78.8 

-.2 

160.8 

38.5 

1 . 3 

18.7 

.05 

0 

0 

80.4 

-.2 

128.1 

38.9 

1.3 

19.4 

.05 

0 

0 

82.0 

-.2 

300.0 

15.6 

1  .  3 

4 . 2 

.05 

0 

0 

90.0 

-.2 

303.0 

16.2 

1 . 3 

4 . 3 

.  05 

0 

0 

90.9 

-.2 

306.0 

16.9 

1 . 3 

4 . 5 

.  05 

0 

0 

91.8 

-.2 

309. 1 

17 . 5 

1 . 3 

4 . 7 

.  05 

0 

0 

92 . 7 

-.2 

312.2 

18.2 

1 . 3 

4.9 

.05 

0 

0 

93 . 7 

-.2 

315.3 

19.0 

1 . 3 

5.1 

.05 

0 

0 

94 . 6 

-.2 

318.5 

19.7 

1.3 

5.3 

.05 

0 

0 

95.5 

-.2 

321.6 

20.5 

1.3 

5.5 

.05 

0 

0 

96.5 

-.2 

324 . 9 

21.3 

1.3 

5.7 

.05 

0 

0 

97.5 

-.2 

328 . 1 

22.2 

1.3 

5.9 

.05 

0 

0 

98 . 4 

-.2 

331.4 

23.1 

1 . 3 

6.2 

.05 

0 

0 

99 . 4 

-.2 

334.7 

24 . 0 

1  .  3 

6.4 

.  05 

0 

0 

100. 4 

-.2 

338.0 

25.0 

1 . 3 

6.7 

.  05 

0 

0 

101.4 

-.2 

341.4 

26.0 

1  .  3 

6 . 9 

.05 

0 

0 

102.4 

-.2 

344.8 

27.0 

1.3 

7.2 

.  05 

0 

0 

103 . 5 

-.2 
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348.3 

28.1 

1.3 

7.5 

.05 

0 

0 

104.5 

-.2 

351.8 

29.2 

1.3 

7.8 

.05 

0 

0 

105.5 

-.2 

355.3 

30.4 

1.3 

8.1 

.05 

0 

0 

106.6 

-.2 

358.8 

31.6 

1.3 

8.4 

.05 

0 

0 

107.7 

-.2 

362.4 

32.9 

1.3 

8.8 

.05 

0 

0 

108.7 

-.2 

366.1 

34.2 

1.3 

9.1 

.05 

0 

0 

109.8 

-.2 

369.7 

35.5 

1.3 

9.5 

.05 

0 

0 

110.9 

-.2 

373.4 

37.0 

1.3 

9.9 

.05 

0 

0 

112.0 

-.2 

377 . 1 

38.4 

1.3 

10.3 

.  05 

0 

0 

113.1 

-.2 

380.9 

40.0 

1.3 

10.7 

.05 

0 

0 

114 . 3 

-.2 

384.7 

41.6 

1.3 

11.1 

.05 

0 

0 

115.4 

-.2 

Range(nm) 

0.4 

104 

.  6 

-.03 

Accuracy (%) 

0.3 

87 

.  2 

-.02 

Availability (%) 

0.3 

93 

.5 

-.03 

95  ! 

95  95 

100 

100 

100 

1 

10 

110  110  110 

110  110  110 

110 

110 

110 

110 

110  110  110 

110  110  110 
110  110  110 


80  80 

80 

80  80 

80  80  80  80 

80 

80  80 

80 

80  80 

80 

80  80 

80  80  80  80 

80 

80  80 

80 

70  70 

70 

72  72 

72  72  75  75 

75 

75  75 

75 

75  75 

75 

75  75 

75  75  75  75 

75 

75  75 

75 

120  - 

.  03 

123  - 

.03  128  -.03 

134 

-.03 

140 

-.03 

145 

-.03 

150  - 

.03 

154  - 

.03  158  -.03 

162 

-.03 

166 

-.03 

169 

-.03 

172  - 

.03 

175  - 

.03  178  -.03 

179 

-.03 

180 

-.03 

180 

-.03 

180  - 

.03 

180  - 

.03  180  -.03 

180 

i  -.03 

180 

-.03 

180 

-.03 

180  - 

.03 

180  - 

.03 

95.0 

-.02 

95.1 

-.02  95.2  -. 

02 

95.3 

-.02 

95.4 

-.02 

95.5 

-.02 

95.6 

-.02 

95.7 

-.02  95.8  -. 

02 

95.9 

-.02 

96.0 

-.02 

96.1 

-.02 

96.2 

-.02 

96.3 

-.02  96.4  -. 

02 

96.4 

-.02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02  96.8  -. 

02 

96.8 

-.02 

96.9 

-.02 

96.9 

-.02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03  95.2  -. 

03 

95.3 

-.03 

95.4 

-.03 

95.5 

-.03 

95.6 

-.03 

95.7 

-.03  95.8  -. 

03 

95.9 

-.03 

96.0 

-.03 

96.1 

-.03 

96.2 

-.03 

96.3 

-.03  96.4  -. 

03 

96.4 

-.03 

96.5 

-.03 

96.6 

-.03 

96.7 

-.03 

96.7 

-.03  96.8  -. 

03 

96 . 8 

-.03 

96.9 

-.03 

96.9 

-.03 

96.9 

-.03 

97 . 0 

-.03 

140  - 

.03 

144  - 

.03  147  -.03 

150 

i  -.03 

153 

-.03 

155 

-.03 

158  - 

.03 

160  - 

.03  162  -.03 

164 

-.03 

165 

-.03 

167 

-.  03 

168  - 

.03 

169  - 

.03  170  -.03 

171 

-.03 

172 

-  .  03 

173 

-.03 

173  - 

.  03 

174  - 

.03  175  -.03 

175 

-.03 

176 

-.03 

176 

-.03 

176  - 

.03 

177  - 

.03 

95.0 

-.02 

95.1 

-.02  95.2  -. 

02 

95.3 

-.02 

95.4 

-.02 

95.5 

-.02 

95.6 

-.02 

95.7 

-.02  95.8  -. 

02 

95.9 

-.02 

96.0 

-.02 

96.1 

-.02 

96.2 

-.02 

96.3 

-.02  96.4  -. 

02 

96.4 

-.  02 

96.5 

-.02 

96.6 

-.02 

96.7 

-.02 

96.7 

-.02  96.8  -. 

02 

96.8 

-.02 

96.9 

-.02 

96.9 

-.02 

96.9 

-.02 

97.0 

-.02 

95.0 

-.03 

95.1 

-.03  95.2  -. 

03 

95.3 

-.03 

95.4 

-.03 

95.5 

-.03 

95.6 

-.03 

95.7 

-.03  95.8  -. 

03 

95.9 

-.03 

96.0 

-.03 

96 . 1 

-.03 

96.2 

-.03 

96.3 

-.03  96.4  -. 

03 

96.4 

-.03 

96 . 5 

-.03 

96 . 6 

-.03 

96.7 

-.03 

96 . 7 

-.03  96.8  -. 

03 

96.8 

-.03 

96.9 

-.03 

96 . 9 

-.03 

96.9 

-.03 

97 . 0 

-.03  . 
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APPENDIX  C:  SAMPLE  OUTPUT  FILES 


Dynamic  End-Of-Lif e-Cycle  Planning  Model 
Case  of  Technological  Change 

Approach  -  Lagrangian  Budget  =  $  601.74 

Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 
ignored  in  the  other  approaches. 

Iterations  of  the  Lagrangian  Relaxation  Method 


Lamdal 

Lamda3 

Lamda2  Slopel 

Slope3 

Slope2 

LI 

L3  L2 

0. 

.  025 

2.080  -306.09 

-32.79 

27 . 04 

18. 

8 

11.4  66. 

.025 

.029 

2.080  -32.79 

-.01 

27.04 

11 . 

4 

11.4  66. 

.  029 

.  032 

2.080  -.01 

27.04 

27.04 

11. 

4 

11.4  66. 

.  029 

.032 

.032  -.01 

27.04 

27.04 

11. 

4 

11.4  11. 

Network 

of 

:  Possible  Solutions  in  Nodal  Form 

Note  -  ' 

rhe 

!  criterion  function  of  a 

node  will 

be  very 

low 

if 

all  paths 

leading 

to  it  fail  to  satisfy  lead 

time  req. 

and/or  : 

minimum 

MOE  levels. 

No.  Sys. 

Year 

EP  Crit.Fn.  Best . Pred . Node  Best 

.Pred. Arc  Best. Alt 

1 

0 

1989 

0  0. 

2 

1 

1989 

1  1.411 

1 

1 

0 

3 

1 

1990 

2  .649 

1 

2 

0 

4 

1 

1991 

3  -.311 

1 

3 

0 

5 

1 

1992 

4  -1.374 

1 

4 

0 

6 

1 

1993 

5-1 . 000E+007 

1 

5 

0 

7 

1 

1994 

6-1 . 000E+007 

1 

6 

0 

8 

2 

1994 

7  -4.107 

2 

7 

1 

9 

2 

1995 

9  -4.476 

3 

11 

1 

10 

2 

1996 

13  -5.013 

4 

17 

1 

11 

2 

1997 

19  -5.659 

5 

25 

1 

12 

2 

1998 

27  -6.302 

5 

33 

1 

13 

2 

1999 

37  -6.876 

5 

43 

1 

14 

2 

2000 

49  -7.378 

5 

53 

1 

15 

2 

2001 

59-1 . 000E+007 

4 

60 

2 

16 

2 

2002 

67-1 . OOOE+O07 

5 

67 

1 

17 

2 

2003 

73-1 .000E+008 

0 

0 

0 

18 

2 

2004 

77-1 . 000E+008 

0 

0 

0 

19 

3 

1999 

79  -7.772 

8 

79 

1 

20 

3 

2000 

81  -7.772 

9 

83 

1 

21 

3 

2001 

85  -7.962 

10 

89 

1 

22 

3 

2002 

91  -8.285 

11 

97 

1 

23 

3 

2003 

99  -8.598 

12 

107 

1 

24 

3 

2004 

109  -8.860 

13 

119 

1 

25 

3 

2005 

121  -9.075 

14 

131 

1 

26 

3 

2006 

133  -9.316 

14 

141 

1 

27 

3 

2007 

145  -9.510 

14 

151 

1 

28 

3 

2008 

157  -9.657 

14 

161 

1 

29 

3 

2009 

169  -9.755 

14 

171 

1 

30 

3 

2010 

181-1. 000E+007 

14 

182 

2 

31 

3  2011 

191-1. O0OE+0O8 

15 

191 

1 

32 

3  2012 

199-1. 000E+008 

16 

199 

1 

33 

3  2013 

205-1 . 000E+008 

0 

0 

0 

34 

3  2014 

209-1. 000E+008 

0 

0 

0 

Solution  Path 

Transition 

Alternative 

Life-Cycle 

Life-Cycle 

Criterion 

Node 

Year 

System  Employed 

Cost 

Effectiveness 

Func. Value 

19 

1999 

3  1 

273.366 

5.029 

-7.772 

8 

1994 

2  1 

321.339 

4.702 

-4 . 107 

2 

1989 

1  0 

-20.000 

.775 

1.411 

The  total  life-cycle  cost  and  effectiveness  for  the  optimal  solution 
are  574.70  and  10.506  respectively. 

Optimal  value  of  the  Lagrangian  multiplier  is  .03180 
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Dynamic  End-Of-Life-Cycle  Planning  Model 
Case  of  Upgrading . Schedul ing 


Approach  -  Satisficing  Budget  =  $  780.08 

Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 
ignored  in  the  other  approaches. 


Sys . Code 
1 
2 
3 


Min. Life 
0 
5 
5 


Max . Life 
5 

10 

10 


Base. Yr. 
1989 
1994 
1999 


No. Nodes 
6 

11 

16 

Lead  Time 
0 
0 


No. Arcs 
6 
72 
132 


Alternative  Name 
Alt.  A 
Alt .  B 


Alt. Code 
1 
2 


Measure  of  Effectiveness 
Code  Name  Weightage 

1  Range(nm)  .400 

2  Accuracy (%)  .300 

3  Availability (%)  .300 


Present  System 
Level  Exp. RateOf Change 

104.600  -.030 

87.200  -.020 

93.500  -.030 


Candidate  Arcs  For  Solution  Path 


Note  -  LCC  is  made  to  be  unfavorably  high  if  the  minimum  reguirement 
level  for  any  component  MOE  of  LCE  has  not  been  met.  The  limiting  MOE 
is  given  as  the  highest  index  of  the  M0E(s)  that  failed  to  meet  the 
minimum  reguirement  level (s)  during  the  life-span  covered  by  the  arc. 


Arc 

Alt. 

Tail 

Head 

Li fe-Cycle 

Limiting 

No. 

Code 

Node 

Node 

Effectiveness  Cost 

MOE 

1 

0 

1 

2 

.7745  -20.00 

0 

2 

0 

1 

3 

1.5287  27.67 

0 

3 

0 

1 

4 

2.2632  80.95 

0 

4 

0 

1 

5 

2 . 97839999999 . 00 

1 

5 

0 

1 

6 

3.67489999999.00 

1 

6 

0 

1 

7 

4 .35309999999 .00 

2 

7 

1 

2 

8 

4.7024  302.34 

0 

8 

2 

2 

8 

4.9501  385.58 

0 

9 

1 

2 

9 

5.41609999999.00 

1 

10 

2 

2 

9 

5.7008  415.63 

0 

11 

1 

3 

9 

4 .74309999999.00 

1 

12 

2 

3 

9 

5.0031  374.09 

0 

13 

1 

2 

10 

6.11099999999.00 

1 

14 

2 

2 

10 

6.4318  443.29 

0 

15 

1 

3 

10 

5.46269999999.00 

1 

16 

2 

3 

10 

5.7617  403.70 

0 

17 

1 

4 

10 

4.8083  299.01 

0 

18 

2 

4 

10 

5.0436  362.25 

0 

19 

1 

2 

11 

6.78779999999.00 

3 

20 

2 

2 

1 1 

7 .14369999999.00 

3 

21 

1 

3 

11 

6.16359999999.00 

1 

131 


22 

2 

3 

11 

6.5004 

430.97 

0 

23 

1 

4 

11 

5.53779999999.00 

1 

24 

2 

4 

11 

5.8083 

391.28 

0 

25 

1 

5 

11 

4 . 8860 

270.43 

0 

26 

2 

5 

11 

5.0842 

351.22 

0 

27 

1 

2 

12 

7.44699999999.00 

3 

28 

2 

2 

12 

7 . 83679999999 . 00 

3 

29 

1 

3 

12 

6.84609999999.00 

3 

30 

2 

3 

12 

7 .21969999999 . 00 

3 

31 

1 

4 

12 

6.24819999999.00 

1 

32 

2 

4 

12 

6.5529 

418.01 

0 

33 

1 

5 

12 

5.6271 

313.96 

0 

34 

2 

5 

12 

5.8550 

379.75 

0 

35 

1 

6 

12 

4.9637 

278.45 

0 

36 

2 

6 

12 

5.1247 

340.96 

0 

37 

1 

2 

13 

8.08899999999.00 

3 

38 

2 

2 

13 

8.51169999999.00 

3 

39 

1 

3 

13 

7.51089999999.00 

3 

40 

2 

3 

13 

7.92009999999.00 

3 

41 

1 

4 

13 

6.93999999999.00 

3 

42 

2 

4 

13 

7.27799999999.00 

3 

43 

1 

5 

13 

6.34879999999.00 

1 

44 

2 

5 

13 

6.6054 

406.06 

0 

45 

1 

6 

13 

5.7164 

320.43 

0 

46 

2 

6 

13 

5.9016 

369.09 

0 

47 

1 

7 

13 

5.0290 

269.00 

0 

48 

2 

7 

13 

5.1529 

330.40 

0 

49 

1 

3 

14 

8.15829999999.00 

3 

50 

2 

3 

14 

8.60209999999.00 

3 

51 

1 

4 

14 

7.61369999999.00 

3 

52 

2 

4 

14 

7.98389999999.00 

3 

53 

1 

5 

14 

7.05159999999.00 

3 

54 

2 

5 

14 

7.33629999999.00 

3 

55 

1 

6 

14 

6.4494 

359.70 

0 

56 

2 

6 

14 

6.6580 

395.03 

0 

57 

1 

7 

14 

5.7915 

309.51 

0 

58 

2 

7 

14 

5.9340 

357.99 

0 

59 

1 

4 

15 

8.26979999999.00 

3 

60 

2 

4 

15 

8.67139999999.00 

3 

61 

1 

5 

15 

7.73599999999.00 

3 

62 

2 

5 

15 

8.04779999999.00 

3 

63 

1 

6 

15 

7.1632 

396.23 

0 

64 

2 

6 

15 

7.3944 

418.86 

0 

65 

1 

7 

15 

6.5340 

347.41 

0 

66 

2 

7 

15 

6.6945 

383 . 44 

0 

67 

1 

5 

16 

8.40249999999.00 

3 

68 

2 

5 

16 

8.74059999999.00 

3 

69 

1 

6 

16 

7.85829999999.00 

3 

70 

2 

6 

16 

8. 11159999999.00 

3 

71 

1 

7 

16 

7 .2570 

382 . 67 

0 

72 

2 

7 

16 

7.4349 

406.84 

0 

73 

1 

6 

17 

8.53509999999.00 

3 

74 

2 

6 

17 

8.80989999999.00 

3 

75 

1 

7 

17 

7.96109999999.00 

3 

76 

2 

7 

17 

8. 15599999999.00 

3 
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77 

1 

7 

18 

78 

2 

7 

18 

79 

1 

8 

19 

80 

2 

8 

19 

81 

1 

8 

20 

82 

2 

8 

20 

83 

1 

9 

20 

84 

2 

9 

20 

85 

1 

8 

21 

86 

2 

8 

21 

87 

1 

9 

21 

88 

2 

9 

21 

89 

1 

10 

21 

90 

2 

10 

21 

91 

1 

8 

22 

92 

2 

8 

22 

93 

1 

9 

22 

94 

2 

9 

22 

95 

1 

10 

22 

96 

2 

10 

22 

97 

1 

11 

22 

98 

2 

11 

22 

99 

1 

8 

23 

100 

2 

8 

23 

101 

1 

9 

23 

102 

2 

9 

23 

103 

1 

10 

23 

104 

2 

10 

23 

105 

1 

11 

23 

106 

2 

11 

23 

107 

1 

12 

23 

108 

2 

12 

23 

109 

1 

8 

24 

110 

2 

8 

24 

111 

1 

9 

24 

112 

2 

9 

24 

113 

1 

10 

24 

114 

2 

10 

24 

115 

1 

11 

24 

116 

2 

11 

24 

117 

1 

12 

24 

118 

2 

12 

24 

119 

1 

13 

24 

120 

2 

13 

24 

121 

1 

9 

25 

122 

2 

9 

25 

123 

1 

10 

25 

124 

2 

10 

25 

125 

1 

1 1 

25 

126 

2 

11 

25 

127 

1 

12 

25 

128 

2 

1  2 

25 

129 

1 

1  3 

25 

130 

2 

13 

25 

131 

1 

14 

25 

8.64659999999.00 

3 

8.85799999999.00 

3 

5.0290 

255.72 

0 

5.1529 

337.32 

0 

5.7915 

296.82 

0 

5.9340 

365.87 

0 

5.0943 

263.82 

0 

5.1934 

327.16 

0 

6.5340 

335.28 

0 

6.6945 

392.24 

0 

5.8667 

303.43 

0 

5.9806 

355.21 

0 

5.1472 

254 . 57 

0 

5.2216 

317.61 

0 

7 .2570 

371.05 

0 

7.4349 

416.48 

0 

6.6187 

340.50 

0 

6.7470 

381.12 

0 

5.9275 

292 . 74 

0 

6.0130 

345.21 

0 

5.2002 

229.69 

0 

5.2497 

308.53 

0 

7.96109999999.00 

3 

8.15599999999.00 

3 

7.3509 

374.97 

0 

7.4932 

404.96 

0 

6.6872 

328.45 

0 

6.7835 

370.73 

0 

5.9884 

266.54 

0 

6.0454 

335.72 

0 

5.2531 

237.36 

0 

5.2779 

299.28 

0 

8 . 64659999999 . 00 

3 

8.85799999999.00 

3 

8 . 06389999999.00 

3 

8.21979999999. 00 

3 

7.4269 

361.66 

0 

7.5337 

394.22 

0 

6.7558 

301 . 03 

0 

6 .8200 

360.88 

0 

6.0493 

272 .86 

0 

6.0778 

325.99 

0 

5.3060 

229.22 

0 

5.2936 

290.87 

0 

8 . 75809999999 . 00 

3 

8 . 92719999999 . 00 

3 

8 . 14729999999 . 00 

3 

8 .26419999999 . 00 

3 

7 . 5030 

333.10 

0 

7.5742 

384.04 

0 

6.8244 

306.08 

0 

6.8565 

350.70 

0 

6.1101 

263.45 

0 

6.0959 

317.19 

0 

5.3466 

206.34 

0 

i  :n 


132 

2 

14 

25 

133 

1 

10 

26 

134 

2 

10 

26 

135 

1 

11 

26 

136 

2 

11 

26 

137 

1 

12 

26 

138 

2 

12 

26 

139 

1 

13 

26 

140 

2 

13 

26 

141 

1 

14 

26 

142 

2 

14 

26 

143 

1 

15 

26 

144 

2 

15 

26 

145 

1 

11 

27 

146 

2 

11 

27 

147 

1 

12 

27 

148 

2 

12 

27 

149 

1 

13 

27 

150 

2 

13 

27 

151 

1 

14 

27 

152 

2 

14 

27 

153 

1 

15 

27 

154 

2 

15 

27 

155 

1 

16 

27 

156 

2 

16 

27 

157 

1 

12 

28 

158 

2 

12 

28 

159 

1 

13 

28 

160 

2 

13 

28 

161 

1 

14 

28 

162 

2 

14 

28 

163 

1 

15 

28 

164 

2 

15 

28 

165 

1 

16 

28 

166 

2 

16 

28 

167 

1 

17 

28 

168 

2 

17 

28 

169 

1 

13 

29 

170 

2 

13 

29 

171 

1 

14 

29 

172 

2 

14 

29 

173 

1 

15 

29 

174 

2 

15 

29 

175 

1 

16 

29 

176 

2 

16 

29 

177 

1 

17 

29 

178 

2 

17 

29 

179 

1 

18 

29 

180 

2 

18 

29 

181 

1 

14 

30 

182 

2 

14 

30 

183 

1 

15 

30 

184 

2 

15 

30 

185 

1 

16 

30 

186 

2 

16 

30 

5.3218  282.34  0 

8.84849999999.00  3 

8.97529999999.00  3 

8.23059999999.00  3 

8.30849999999.00  3 

7.5791  336.97  0 

7.6146  373.47  0 

6.8930  295.48  0 

6.8769  341.56  0 

6.1568  239.38  0 

6.1283  308.21  0 

5.3871  213.78  0 

5.3376  274.17  0 

8.93889999999.00  3 

9.02339999999.00  3 

8.31389999999.00  3 

8.35289999999.00  3 

7.6551  325.27  0 

7.6373  364.02  0 

6.9455  270.30  0 

6.9134  332.17  0 

6.2034  245.62  0 

6.1464  299.60  0 

5.4277  206.71  0 

5.3533  266.52  0 

9.02929999999.00  3 

9.07159999999.00  3 

8.39719999999.00  3 

8.37779999999.00  3 

7.7134  299.05  0 

7.6778  354.26  0 

6.9981  275.42  0 

6.9339  323.17  0 

6.2500  237.47  0 

6.1646  291.56  0 

5.4682  185.59  0 

5.3691  258.95  0 

9.11969999999.00  3 

9.09849999999.00  3 

8.46109999999.00  3 

8.42209999999.00  3 

7.7717  303.13  0 

7.7005  344.91  0 

7.0506  266.25  0 

6.9543  314.78  0 

6.2967  215.24  0 

6.1827  283.58  0 

5.4806  192.99  0 

5.3815  251.77  0 

9.18889999999.00  3 

9.14669999999.00  3 

8.52489999999.00  3 

8.44699999999.00  3 

7.8299  293.02 

7.7232  336.21 
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o  o 


187 

1 

17 

30 

7.1031 

242.97 

0 

188 

2 

17 

30 

6.9748 

306.43 

0 

189 

1 

18 

30 

6.3109 

221.65 

0 

190 

2 

18 

30 

6.1970 

276.03 

0 

191 

1 

15 

31 

9.25819999999.00 

3 

192 

2 

15 

31 

9.17359999999.00 

3 

193 

1 

16 

31 

8.58879999999.00 

3 

194 

2 

16 

31 

8.47179999999.00 

3 

195 

1 

17 

31 

7.8882 

268.76 

0 

196 

2 

17 

31 

7.7459 

327.52 

0 

197 

1 

18 

31 

7 . 1192 

248.45 

0 

198 

2 

18 

31 

6 . 9908 

298 . 54 

0 

199 

1 

16 

32 

9.32739999999.00 

3 

200 

2 

16 

32 

9 .20059999999 . 00 

3 

201 

1 

17 

32 

8.65259999999.00 

3 

202 

2 

17 

32 

8.49669999999.00 

3 

203 

1 

18 

32 

7.9060 

273 . 39 

0 

204 

2 

18 

32 

7.7637 

319 . 33 

0 

205 

1 

17 

33 

9.39669999999.00 

3 

206 

2 

17 

33 

9.22759999999.00 

3 

207 

1 

18 

33 

8 . 67209999999 . 00 

3 

208 

2 

18 

33 

8.51619999999.00 

3 

209 

1 

18 

34 

9.41779999999.00 

3 

210 

2 

18 

34 

9.24869999999.00 

3 

Network  of  Possible  Solutions  in  Nodal  Form 

Note  -  The  criterion  function  of  a  node  will  be  very  low  if  all  paths 
leading  to  it  fail  to  satisfy  lead  time  req.  and/or  minimum  HOE  levels. 

No.  Sys.  Year  EP  Crit.Fn.  Best . Pred . Node  Best . Pred . Arc  Best. Alt. 


I 

0 

1989 

0  0. 

2 

1 

1989 

1  20.000 

1 

1 

0 

3 

1 

1990 

2  -27.666 

1 

2 

0 

4 

1 

1991 

3  -80.950 

1 

3 

0 

5 

1 

1992 

4-1 . 000E+007 

1 

4 

0 

6 

1 

1993 

5-1 . OOOE+0O7 

1 

5 

0 

7 

1 

1994 

6-1 . 000E+007 

1 

6 

0 

8 

2 

1994 

7  -282.339 

2 

7 

1 

9 

2 

1995 

9  -395.635 

2 

10 

2 

10 

2 

1996 

13  -379.964 

4 

17 

1 

11 

2 

1997 

19  -458.632 

3 

22 

2 

12 

2 

1998 

27  -498.963 

4 

32 

2 

13 

2 

1999 

37-1 . 000E+007 

2 

37 

1 

14 

2 

2000 

49-1 . OOOE+007 

3 

49 

1 

15 

2 

2001 

59-1 . OOOEf 007 

4 

59 

1 

16 

2 

2002 

67-1 . 000E+008 

0 

0 

0 

17 

2 

2003 

73-1 . 000E+008 

0 

0 

0 

18 

2 

2004 

77-1 . 000E+008 

0 

0 

0 

19 

3 

1999 

79  -538.063 

8 

79 

1 

20 

3 

2000 

81  -579.163 

8 

81 

1 

21 

3 

2001 

85  -617.622 

8 

85 

1 

22 

3 

2002 

91  -653.389 

8 

91 

1 

23 

3 

2003 

99  -708.410 

10 

103 

1 
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24 

3  2004 

109  -741.623 

10 

113 

1 

25 

3  2005 

121  -791.733 

11 

125 

1 

26 

3  2006 

133  -835.932 

12 

137 

1 

27 

3  2007 

145-1. 000E+007 

13 

149 

1 

28 

3  2008 

157-1. 000E+007 

12 

157 

1 

29 

3  2009 

169-2. 000E+007 

13 

169 

1 

30 

3  2010 

181-1. OO0E+OO8 

0 

0 

0 

31 

3  2011 

191-1 . OOOE+008 

0 

0 

0 

32 

3  2012 

199-1. 000E+008 

0 

0 

0 

33 

3  2013 

205-1 . 000E+008 

0 

0 

0 

34 

3  2014 

209-1 . 000E+008 

0 

0 

0 

Minimum  Average 

Annual  Cost  Solution 

Path 

(  Time 

value  of 

money  is  ignored,  ie 

.  interest 

rate  is  assumed 

to  be  0  ) 

Transition 

Alternative 

Life-Cycle 

Life-Cycle 

Criterion 

Node 

Year 

System  Employed 

Cost 

Effectiveness 

Func . Value 

24 

2004 

3  1 

361.659 

7.427 

-741.623 

10 

1996 

2  1 

299.014 

4.808 

-379.964 

4 

1991 

1  0 

80.950 

2.263 

-80.950 

Minimum  Annualised  Cost  Solution  Path 

(  Monetary  interest  assumed  constant  at  5.000  %  per  annum  ) 


Transition 

Alternative 

Life-Cycle 

Life-Cycle 

Criterion 

Node 

Year 

System 

Employed 

Cost 

Effectiveness 

Func. Value 

22 

2002 

3 

1 

371.050 

7.257 

-653 .389 

8 

1994 

2 

1 

302 .339 

4 . 702 

-282 .339 

2 

1989 

1 

0 

-20.000 

.775 

20 . 000 
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